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ABSTRACT 


Cracking of asphaltic concrete pavements due to thermal effects 
is recognized as a serious form of pavement distress in areas of low 
temperature climatic environment. In this study pavement temperature 
data obtained from field installations are summarized, an analytical 
model for predicting pavement temperatures is developed and using such 
temperature data together with appropriate materials definition the 
results of various thermal stress prediction methods are assessed in 
relation to observed low temperature behavior of pavement structures 


within two field test projects located in Western Canada. 


The study firstly presents a review of various factors found 
from previous investigations to be related to the low temperature 
cracking problem. Following this review recorded temperatures in 
several pavement structures subjected to temperatures approaching -40 F 
are summarized. The summary includes specific information on daily air 
and pavement temperature variations, time rates of pavement temperature 
change and thermal regimes within subgrades of different pavement 
structures subjected to near identical climatic conditions. To extend 
this information to other climatic areas an analytical model for pre- 
dicting pavement temperatures is developed and necessary input para- 
meters for the model are described. Comparisons are given between pre- 
dicted and recorded temperatures in various pavement structures. Such 


comparisons have shown excellent agreement which suggests that the model 
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can be used to obtain reliable predictions of pavement temperatures 


under actual, climatyve conditions. 


Concepts used to characterize the time and temperature de- 
pendence of asphaltic concretes incorporated in the test projects are de- 
scribed and resulting stiffness values are presented. Tensile properties 
of pavement cores obtained from the projects and determined by means of 
the tensile splitting test are summarized. Using the stiffness values 
together with recorded and predicted temperature data thermally induced 
stresses within the pavements are computed by various stress predictive 
methods. Theoretical considerations for the formulation of the predic- 
tive methods are reviewed and numerical solutions of the stress equations 
are presented. Employing a postulated fracture criterion, computed 
stresses are compared with laboratory evaluated fracture strength- 
temperature "relationships of “the pavement “cores: to’ predict times and 
temperatures of fracture.” These predicted ‘fracture conditions are then 
compared with times and temperatures of observed cracking of the corre- 
sponding pavements. Based on these comparisons a simplified guideline 


for reducing the occurrence of transverse cracking is proposed. 


The principal conclusion of this study is that stress pre- 
diction methods coupled with the use of recorded or predicted temperature 
data and appropriate materials characterization provide a means of 
reliably assessing the low temperature fracture susceptibility of 
asphaltic concrete mixes. For field data currently available such 
assessments can be made using thermal stresses computed by means of a 


suitable pseudo-elastic beam analysis. 
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CHAPTER I 


INTRODUCTION 


Pry Introductory Remarks 


The performance of asphaltic concrete pavements is dependent on 
many individual and interacting load, environment, construction and 
maintenance variables. In areas subjected to low temperature climatic 
environments, significant deterioration in riding quality and increase 
in maintenance cost of these pavements, commonly referred to as flexible 
pavements, have been attributed, in part, to non-traffic load associated 
distress mechanisms. One manifestation of such distress is low tempera- 
ture transverse crackine. During, the past.decade: many agencies ,and 
individuals have expended considerable effort towards developing 
solutions to the transverse cracking problem. The majority of investi- 
gations related to the problem have involved the collection of detailed 
field inventory data, observation of controlled field test projects and 
the testing of pavement component materials under various laboratory 


conditions. 


Documented field inventories and observations of controlled 
field test projects have enabled many significant factors associated 
with low temperature cracking of asphaltic concrete pavements to be 


defined. ‘These test projects have also yielded specific information on 
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the temperature distributions within pavement structures during periods 
of observed fracture. The low temperature properties of constituent 
materials of pavements, determined by various material characterization 
techniques, have been compared with observed field behavior. Such com- 
parisons have permitted the influence of material properties on the low 
temperature fracture susceptibility of asphaltic concrete mixtures to be 
assessed and have enabled many factors contributing to subsequent 


pavement performance losses to be identified. 


The form of pavement distress associated with low temperature 
cracking is but one of a series of distress mechanisms that must be 
considered in an overall pavement design process. Since a pavement is a 
complex structure, subjected to a variety of loading and environmental 
conditions, the solution of the overall pavement design problem requires 
the development of a comprehensive and coordinated framework in which 
various components of the problem can be organized. In recent years, 
systems engineering principles, such as those presented by Hudson et al. 
(1969) and Hutchinson and Haas (1968), have been employed to structure 


pavement design processes. 


The significant elements of the pavement design system proposed 
by Hudson et al. (1969) are shown in FIGURE I-1 and include: 1) indi- 
vidual and interacting load, environment, construction and maintenance 
input variables; 2) physical characteristics of pavements, which include 
the basic properties which characterize material behavior; 3) primary 
and limiting output functions of the system; 4) a decision criteria 


based on economic and performance considerations; and 5) feedback and 
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interaction between the various units of the system. 


The pavement design process proposed by Hutchinson and Haas 
(1968) is subdivided into six major phases which include: 1) problem 
definition; 2) solution generation; 3) solution analysis; 4) evaluation 
and optimization; 5) implementation; and 6) performance assessment. 
These authors reviewed existing knowledge within each of the six phases 
of the design system and concluded that one of the most important de- 
ficiencies in current knowledge is a quantitative understanding of non- 
load associated factors which cause progressive deterioration in 


pavement serviceability. 


Within the proposed pavement design systems, component design 
subsystems may be structured in which one or a combination of forms of 
pavement distress are considered. One such subsystem which considers 
the response of flexible pavements at low temperatures has been pro- 
posed by Haas and Anderson (1969). These authors organized existing 
knowledge of pavement distress associated with low temperature cracking 
within the design framework proposed by Hutchinson and Haas (1968), and 
concluded that the primary aspect of low temperature response of 
flexible pavements that must be considered in such a design subsystem 


is their susceptibility to fracture. 


With reference to the subsystem presented by Haas and Anderson 
(1969), possible causes of low temperature cracking were grouped into 


the following three categories: 


1) Exceeding the tensile strength (or tolerable strain) of the 


bituminous surface as a result of thermally induced stresses 
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2 ) 


3) 


(and strains), 
a) without considering traffic loads, 


b) when considering traffic loads. 


Freezing shrinkage and cracking of the subgrade, followed by 
crack propagation through the base and asphaltic concrete 


component layers. 


Freezing shrinkage and cracking of the base or subbase 
materials, followed by crack propagation through the 


asphaltic concrete. 


Considering category l(a), the major elements of the solution 


analysis and evaluation phases of the aforementioned subsystem are shown 


in FIGURE I-2. As illustrated in this figure, an assessment of the low 


temperature fracture susceptibility of asphaltic concrete paving 


mixtures requires a knowledge of the thermal regimes to which asphaltic 


concrete are subjected, definition of the low temperature stiffness and 


strength properties of asphaltic concretes, and a model by which the low 


temperature response of paving mixtures can be realistically evaluated. 


1.2. Purpose and Scope 


The objectives of this investigation are: 


1) 


2) 


to summarize, in a statistical manner, recorded temperatures 
in asphaltic concrete pavement structures subjected to low 


temperature climatic environments; 


to develop an analytical model for predicting temperatures 


in asphaltic concrete pavement structures and to compare 
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2) 


3) 


(and strains), 
a) without considering traffic loads, 


b) when considering traffic loads. 


Freezing shrinkage and cracking of the subgrade, followed by 
crack propagation through the base and asphaltic concrete 


component layers. 


Freezing shrinkage and cracking of the base or subbase 
materials, followed by crack propagation through the 


asphaltic concrete. 


Considering category l(a), the major elements of the solution 


analysis and evaluation phases of the aforementioned subsystem are shown 


in FIGURE I-2. As illustrated in this figure, an assessment of the low 


temperature fracture susceptibility of asphaltic concrete paving 


mixtures requires a knowledge of the thermal regimes to which asphaltic 


concrete are subjected, definition of the low temperature stiffness and 


strength properties of asphaltic concretes, and a model by which the low 


temperature response of paving mixtures can be realistically evaluated. 


i532, Purpose and Scope 


The objectives of this investigation are: 


1) 


2) 


to summarize, in a statistical manner, recorded temperatures 
in asphaltic concrete pavement structures subjected to low 


temperature climatic environments; 


to develop an analytical model for predicting temperatures 


in asphaltic concrete pavement structures and to compare 
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predicted and recorded temperatures; 


3) to characterize the stiffness and low temperature tensile 
strengths of various asphaltic concrete paving mixtures 
which are known to exhibit different low temperature 


behavior; 


4) to assess various stress predictive methods by comparing 
predicted and observed field fracture times and tempera- 


tures of various asphaltic concrete pavement surfaces; and 


5) to develop quantitative guidelines for the design of 
asphaltic concrete paving mixtures subjected to climatic 
conditions, where, neglecting all other loading mechanisms, 
the thermally induced tensile stress may exceed the tensile 


strength of the asphaltic concrete. 


In the present study data obtained from two extensive field test 
projects in Western Canada, designed to study the low temperature re- 
sponse of asphaltic concrete pavements, have been analyzed to meet these 


objectives. 


i.3, Outline of Contents 
The results of previous laboratory and ficld investigations re- 
lated to the low temperature cracking problem of asphaltic concrete 


pavements are briefly reviewed in CHAPTER II. 


In CHAPTER III recorded temperatures in asphaltic concrete 
pavements, during selected time periods of low temperature climatic 


environments, are summarized. The summary includes daily temperature 
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6 
variations, temperature gradients and time rates of temperature change 


in pavement surfaces and subgrade thermal regimes. 


CHAPTER IV contains a review of heat transfer concepts as 
applied to pavement structures. Methods employed to quantify the in- 
fluence of various heat transfer modes on temperature regimes of 
pavement structures are discussed and thermal properties of materials 


comprising such structures are presented. 


Assuming a one-dimensional heat conduction problem, a finite 
difference model for predicting temperatures in layered systems, in 
response to imposed climatic conditions, is developed in CHAPTER V. 
Details of the temperature prediction computer program are contained in 
APPENDIX A. Computed temperatures are compared with measured tempera- 


tures in various pavement structures. 


Tensile strength and stiffness properties of various asphaltic 
concrete paving mixtures are presented in CHAPTER VI. These properties 
are discussed in relation to observed low temperature behavior of 


pavements in which the mixtures are incorporated. 


In CHAPTER VII, numerical methods used to predict thermal 
stresses in asphaltic concrete pavements are discussed and theoretical 
considerations for the formulation of the stress prediction methods are 
reviewed. The various stress analyses are then used to compute 
thermally induced stresses within the asphaltic concretes whose strength 
and stiffness properties are given in CHAPTER VI. Details of the com- 
puter programs developed for these stress computations are contained in 


APPENDICES B, C, D and E. Employing a postulated fracture criterion, 
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7 
predicted times and temperatures of fracture of the paving mixtures are 
compared with times and temperatures of observed cracking of the corre- 
sponding pavements in the field projects. From the results of such 
comparisons, the stress predictive methods are assessed and a corre- 
lation between observed and predicted low temperature fracture suscepti- 


bidity of the “asphaltic concretes is derived. 


Utilizing) field data currently available, together with the 
results of the comparisons between observed and predicted pavement 
fracture, presented in CHAPTER VII, tentative guidelines for the design 
of asphaltic concrete paving mixtures with respect to low temperature 
cracking are developed in CHAPTER VIII. A short summary, the con- 
clusions of this investigation, and recommendations for further study 


are contained in CHAPTER IX. 
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FIGURE I-2 MAJOR ELEMENTS OF THE SOLUTION ANALYSIS 
AND EVALUATION PHASES FOR THE DESIGN SUBSYSTEM. 
(after Haas and Anderson, I969 ) 


CHAPTERed1 


A SELECTED REVIEW OF VARIOUS FACTORS 
RELATED TO LOW TEMPERATURE CRACKING OF 


FLEXIBLE PAVEMENTS 


2.1. Introduction 

Non-load associated fracture of flexible pavements subjected 
to low temperature climatic environments, such as experienced in the 
Northern United States, was recognized as a serious problem by Rader 
(1935). The problem of low temperature cracking became acute in 
Canada during the mid-1950's and early 1960's when the highway network 
was rapidly expanding and performance requirements for pavements in- 
creased. During this period of time both field and laboratory investi- 
gations were initiated in an attempt to define the extent and the causes 
of transverse cracking. Among the first results of these early studies 
were those reported by Shields (1964), Shields and Anderson (1964), and 


Domaschuk et al. (1964). 


In 1965, the Canadian Good Roads Association (CGRA) recognized 
the severity of the transverse cracking problem and placed it within the 
top priority of research needs in Canada. The paper presented by 
Monismith et al. (1965) to the Association of Asphalt Paving Technolo- 


gists concerning thermally induced stresses in asphaltic concrete paving 
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mixtures and the Association's 1966 Symposium on Non-Traffic Load 
Associated Cracking gave further recognition to the problem. Since the 
mid 1960's field test projects designed to study the low temperature 
response of asphaltic concrete pavements have been constructed and 
various laboratory studies related to the problem have been either 
initiated or continued. An ad hoc committee on Low Temperature Behavior 
of Flexible Pavements was established within CGRA in 1969, the purpose 
of which was to dimension the low temperature pavement cracking problem 
in Canada and to document used or suggested design approaches to the 
problem. The results of this comprehensive review of the subject area 
has been presented by Haas et al. (1970). In view of this comprehensive 
review, the content of this chapter is limited to a brief description of 
low temperature transverse cracks, a review of performance and economic 
considerations pertaining to the problem, and a summary of some vari- 
ables which have been found to be associated with transverse crack 


development. 


Sia Description of Low Temperature Transverse Cracks 


In areas of low temperature climatic environment, the fracture 
pattern of asphaltic concrete surfaces placed on granular base courses 
or directly on prepared subgrades is commonly one of regular spaced 
transverse cracks in association with longitudinal cracks along the 
center line and shoulders of paved surfaces. Such a fracture pattern 
differs from the more irregular pattern generally exhibited by pavement 
structures that incorporate portland cement treated or stabilized bases. 


The irregular spacing of transverse and longitudinal cracks within the 
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surface of such structures may partially be the result of hydration 
cracks within the treated base being reflected through the asphaltic 
concrete. The following discussion is confined to transverse cracks 


associated with the former pavement structures. 


Four types of transverse cracks, as identified by McLeod (1968), 
are shown in FIGURE II-1. Each type is commonly found in pavements 
where low temperature cracking has become a serious problem. From con- 
tinuous crack mapping of different pavement structures included in an 
extensive field test project designed to study the low temperature 
response of asphaltic concrete pavements, Young et al. (1969) found that 
the pattern of transverse crack development often follows a systematic 
process. This process involves the formation of initial cracks spaced 
at relatively large intervals followed by the development of subsequent 
transverse cracks spaced at approximately regular intervals between the 
initially formed cracks. As transverse crack frequency increases longi- 
tudinal cracks may develop in the wheel paths causing a block fracture 
pattern to develop. Such a fracture pattern of pavements located in 


Northern Ontario was reported by Kingham (1966). 


Results of transverse crack surveys, such as presented by 
Culley (1966), Shields (1964), Shields et al. (1964) and (1969), 
Anderson et al. (1966), Fromm (1966), and Young et al. (1969), have 
shown that within a given locale asphaltic concrete pavements commonly 
exhibit a wide range in transverse crack frequencies. Typical of the 
range experienced is that reported by Shields (1964). This particular 


survey included all main paved highways within a fifty mile radius of 


iat 
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13 
the City of Edmonton, Alberta, in which crack frequencies were found 
to range from pavements exhibiting no transverse cracks to others 
exhibiting 400 cracks per mile. Factors believed contributing to such 
a range in crack frequency are discussed in Section 2.4 of this 


chapter. 


net Gack, An Maat hon 

Continuously monitoring temperature and crack detection 
systems, employed in field test projects located in Southern Manitoba 
and Central Alberta, have enabled temperature conditions conducive to 
transverse cracking of various pavement structures to be defined. 
Young et al. (1969) reported that transverse cracking of pavement 
structures located at the Manitoba test site resulted from the effects 
of prolonged low temperatures with the majority of cracks being 
initiated at a time approximately equivalent to the time at which 
minimum daily air temperatures were recorded. Shields et al. (1969) 
reported that initial transverse cracking at the Alberta test site was 
confined to periods of seasonal minimum pavement temperatures. Further 
evidence of an apparent relationship between transverse crack develop- 
ment and low temperature climatic environment has been presented by 
Fromm (1966) and Kingham (1966). These investigators reported that 
pavements located in the northern regions of the Province of Ontario 
commonly exhibit higher crack frequencies than pavements located in 
southern areas of the province. Although it is generally agreed that 
the majority of transverse pavement cracks in Canada initiate during 


severe winter climatic conditions Culley (1966), Domaschuk et al. (1964) 
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and Iluculak (1964), reported that transverse cracks become evident 
after prolonged cold temperatures followed by a sudden rise in 


temperature. 


{ 
Pegaso) crack Width Development 


Among the most extensive studies related to transverse crack 
development have been those reported by Rix (1969), and Shields 
et al. (1969). From measurements between surface gauge points, in- 
Stalled at transverse crack locations within various mature pavement 
structures in Alberta and made using a mechanical extensometer, Rix 
(1969) found that seasonal crack width changes could not be accounted 
for by simple contraction and expansion of the asphaltic concrete 
surfaces. Maximum crack widths tended to increase with increased depth 
of frost penetration within the subgrades of the pavement structures 
and minimum crack widths, measured during summer months, generally 
exceeded minimum seasonal widths of the previous year. Permanent 
crack width increases, for pavements located in Northern Ontario, have 


also been reported by Kingham (1966). 


Employing measurement techniques as used by Rix (1969), Shields 
et al. (1969) reported that initial surface crack widths of a pavement 
in Central Alberta were in the order of 0.4 millimeters. Following a 
period of low temperature, an increase in pavement surface temperature 
resulted in generally larger crack widths and approximately fifty days 
after crack initiation a sudden cold spell resulted in further contrac- 
tion with maximum crack widths in the order of 6.0 millimeters. During 


the second winter crack openings became much more apparent and maximum 
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widths of approximately 12.0 millimeters were observed at a time closely 
related to the time of maximum recorded depth of frost penetration. The 
authors concluded that since a direct relationship between crack opening 
and surface temperature was not apparent some mechanism other than 


Simple surface contraction must be operative. 


With the exception of these investigations, the influence of 
subgrade, base and subbase material properties on crack width develop- 
ment has not been included in field studies directly related to the 
transverse cracking problem. However, Young et al. (1969) reported 
that initial crack widths of 0.19 inches were measured at the surface 
of asphaltic concrete structures incorporating a highly plastic clay 
subgrade, while, surface cracks in structures incorporating a sand 
subgrade were very minute and difficult to detect by visual observation. 
Further evidence of the influence of subgrade soil type on crack width 
development has been presented by Hamilton (1966). Employing uni- 
directional closed system freezing tests, Hamilton (1966) found that 
various soils exhibit shrinkage normal to the direction of freezing and 
that the magnitude of freezing shrinkage is dependent on soil density, 


plasticity and degree of saturation. 


Deore. worack Depth 


Asphaltic concrete coring data supplemented by general field 
observations, such as reported by Huculak (1964), Anderson et al. 
(1966), Kingham (1966), Shields et al. (1969) and Young et al. (1969), 
have indicated that in some cases transverse cracks are limited in depth 


to the asphaltic concrete, while in others, cracks extend through the 
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surface and into subsurface layers. Deep subsurface cracks are gener- 
ally found in pavement structures having clay subgrade soils. In both 
the preceding cases the influence of subsurface layers on promoting 
transverse cracking of the asphaltic concrete is unknown. However, 
from results of detailed field and laboratory studies, Anderson et al. 
(1966) concluded that mechanisms promoting transverse cracking cannot 
be simply and explicitly stated but undoubtedly involves thermal con- 
traction of surface components which may be masked by the behavior of 


the subgrade under certain conditions. 


The possibility of transverse cracking due to thermal con- 
traction of the asphaltic concrete surface has been supported by 
results of experimental studies presented by Monismith et al. (1965), 
Hills and Brien (1966), Haas (1968) and Haas and Topper (1969). 
Employing various material characterization techniques, theoretical 
considerations and stress predictive methods, these investigators have 
all found that at low temperatures predicted thermally induced stresses 
within asphaltic concrete surfaces may be of a magnitude far exceeding 


the tensile strength of asphaltic concrete. 


Dele. VELticau Displacement 


The initial effect of transverse cracking on pavement riding 
quality is generally minor. However, the intrusion of water at crack 
locations into subgrades exhibiting high potential swelling character- 
istics can result in heaving at crack locations. On sand subgrades, 
removal of fine particles by water can result in depressions at crack 


locations. Kingham (1966) presented data showing that the differential 
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heave parallel to transverse cracks can be as much as 0.44 inches, with 
the maximum and minimum vertical movement occurring at the pavement 

edge and pavement center line, respectively. Vertical displacements at 
transverse cracks located in various pavements have been reported by 
Rix (1969). Rix found that such displacements are primarily controlled 
by the moisture content of the subgrade, while the direction of movement 


appears related to soil plasticity as well as moisture content. 


In addition to vertical displacement, the loss of load transfer 
at crack locations can result in additional cracks being formed imme- 
diately adjacent to the original cracks. As cracking progresses, 
spalling may develop with progressive removal of surface material by 
traffic and the eventual formation of potholes. Such events can 
readily lead to rapid pavement performance losses, increased maintenance 


costs and decreased pavement surface life. 


2.3. Performance and Economic Considerations 

Anderson et al. (1966) presented data showing the severe 
performance losses that can occur in pavements exhibiting ineffectively 
maintained transverse cracks. The performance-age history of one such 
pavement is shown in FIGURE II-2. In this figure the predicted per- 
formance-age history was determined by means of a mathematical model 
relating performance to age, strength as measured by deflection and to 
traffic conditions. The model was developed from data collected on 
approximately 469 inventory pavement sections on clay subsoils in 
Alberta. Assuming a terminal serviceability of 4.5, the pavement 


reached this level in approximately one fourth of its expected life 
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span. Resurfacing of such a pavement will raise the level of pavement 
serviceability, however, existing cracks generally reflect through the 
new surface causing increased maintenance problems and continued 


performance losses. 


Haas (1968) reviewed the economic importance of the low 
temperature cracking problem in Canada. Based on approximate 1967 
costs of crack filling operations in the prairie provinces of Canada, 
Haas concluded that a conservative estimate of annual maintenance costs 


of these operations would be at least one million dollars. 


In addition to being a source of maintenance expenditures, the 
reduction in pavement service life caused by transverse cracking may be 
of much more importance in an economic sense. To the author's knowl- 
edge an objective analysis of such costs has not been performed. 
However, based on a simplified economic evaluation of 1966-67 annual 
costs of pavement construction in Western Canada, Haas (1968) con- 
pinded that it is not unreasonable to suggest that low temperature 
cracking of pavements in Canada is an annual multi-million dollar 


problem. 


2.4. Factors Related to Transverse Crack Development 


From the results of extensive field inventory data and field 
observations Shields (1964) postulated possible mechanisms promoting 
transverse cracking which included simple contraction of the asphaltic 
concrete and/or shrinkage of the subgrade. Considering these possible 
mechanisms, together with the wide range of materials and climatic 


conditions encountered, Shields categorized the many factors believed 
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pertinent to transverse crack development by means of FIGURE II-3. 
Factors of possible significance to transverse crack development have 
also been summarized by Haas (1968). In this latter summary, external 
and component variables of possible significance were summarized with 
respect to the pavement component layer in which cracking may initiate. 
Some factors which have been found to be of recognized significance to 
the low temperature cracking problem, and which are reviewed in 


following subsections, include; 


1) asphalt source, 

2) asphalt grade, 

3) asphalt stiffness, 

4) asphaltic concrete fracture characteristics, 
5) asphaltic concrete age, 

6) asphaltic concrete thickness, and 


7) subsurface materials. 


Ls Asphalt Source 


Results of transverse crack frequency surveys have shown that 
one of the most significant variables associated with the occurrence of 
low temperature cracking is the crude source of the asphalt cement, as 
indicated by differences in low temperature field behavior of asphaltic 
concretes composed of similar asphalts but obtained from different 
suppliers. Shields (1964), Shields and Anderson (1964) and Anderson 
et al. (1966) presented crack count data for many miles of the Alberta 
highway system. The data indicated a definite relationship between 


frequency of cracking and asphalt supplier, although exceptions were 
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noted. Fromm (1966) reported a similar relationship for pavements 
within the Province of Ontario. The results of field test projects, 
such as reported by Culley (1966), McLeod (1969) and Shields et al. 
(1969) , have also shown that pavement structures of similar properties 
but composed of asphalts from different crude sources exhibit different 
crack frequencies. Although field observations suggest relationships 
between crack frequencies and asphalt sources, results of extensive 
laboratory investigations such as reported by Anderson et als ((1966);, 
have revealed that individual asphalt properties obtained by conven- 
tional laboratory testing methods cannot be definitely correlated with 


this observed behavior. 


2.4.2. Asphalt Grade 

Incorporated within field test projects in the Provinces of 
Ontario, Manitoba, Saskatchewan and Alberta were asphaltic concrete 
pavement surfaces composed of asphalts of different penetrations as 
measured at 77 F and different absolute viscosity values as determined 
at 140 F. The results of these test projects, presented by McLeod 
(1969), Young et al. (1969), Culley (1966) and Shields et al. (1969), 
respectively, have shown that those asphaltic concretes composed of 
similar penetration grade asphalts but having different viscosities 
exhibit significantly different crack frequencies. Asphaltic concretes 
which incorporated asphalts of the same relative viscosity but of 
different penetrations were also found to exhibit different crack 
frequencies. In the former case, cracking was found to be most 


predominate in the asphaltic concretes composed of the lowest viscosity 


ig || ieee or 
_Rigemevagirsh mienoi seta} we linte « booger REET mMat aie 
th9etor4 2e0F (81889 29 aXiyson 4177 0x00 To sumiverineda alll 
sf 23 ebioade Ban (C001) boesiol(abeT) Xai Lad Yeh bem 
eeitreqorq 1slini2, Yo deer bnss inten sg tad mouidconit « 21 
SMTA Lh sidirixe eearwes Shinn Inqiet¥hh port et lenges to hoeogan 
eqitenoiiniey Idegnie cihitavibeds bioit, dyvatitA © ee. 
owhetorxs to ehizay .2scgo ¢ hneetns bas steel ae 0 
, (O0OL). tg. 39 toershah Yi hetroqes: en. iw stoi tay ha aovae e 
tovies vd hentsydo eaririqory ifeigen louliv iba? 187, bale ms 


mtiw betsiorroo \Ylesintieh sd sonne> ubodtem hissed croimgbdat f 


ToL vated bornsede-e 


J 


Sieiones atpishges s1ew syredts bas Aewadodsten? Shoe tniet 0 
es paoiietiensd tasxsttib 2; farigas 20 boeognas, cobeteaue 

bontineseh 2k eovlev vttzosziy otwloed;; Siretetanh bas. 4) Wh re 

besioM ud Sormeosy <2999( 074 Ja9d d£0d7 Yo ertuege oat FON ay 


. aise ehioide hyve (eaRt) OL Sano e (CO@Ts, ks —— 


; He fetsese snd ode do a tsdqea bedamagere 


es! fib Aidinxs ot bainet pelg ST 9 ie 


ag or ; bro? enw anidero <9AgD. “Tonto od3, 
= “218 20, hezoqmen re teronos abspeiges sie | 


ou 
asphalt, while in the latter, highest crack frequencies were associated 


with pavements incorporating the lowest penetration asphalt. 


Relationships between absolute viscosity and penetration values 
of asphalts, measured at 140 and 77 F, respectively and as presented by 
McLeod (1967), are shown in FIGURE II-4. McLeod denoted lines A, B and 
C as being representative of high, low and intermediate viscosity 
asphalts, respectively. From this figure, the wide range in absolute 
viscosity which may be experienced when grading asphalt by penetration 
and the range in asphalt penetration when grading by viscosity is 
apparent. On the basis of the results of the previously mentioned 
field studies and in particular those presented by Young et al. (1969) 
the low temperature fracture susceptibility of an asphaltic concrete 
paving mixture is dependent on both the penetration and the viscosity 


of the asphalt comprising the mixture. 


Ie: ae Asphalt Stiffness 


In recent years the stiffness concept, introduced by Van der 
Poel (1954), has received fairly extensive attention in relation to the 
low temperature cracking problem. Stiffness, as defined by Van der 
Poel, is a time and temperature dependent ratio of tensile stress to 
strain of asphalt materials subjected to various loading conditions. 
Further details of the stiffness concept are contained in CHAPTER VI. 
Various methods used to determine stiffness moduli of asphalts at low 
temperatures have been summarized by Anderson and Haas (1970) and 
Hajek (1971). These methods may be classified as direct methods, 


based on direct testing of the asphalt, and indirect methods which make 
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use of nomographs for estimating low temperature stiffness values from 


a knowledge of asphalt properties determined at higher temperatures. 


Employing a direct tension test on thin films of two typical 
asphalts used in Canada, Haas (1968) found that the stiffness moduli 
thus determined could be subjectively related to observed differences 
in low temperature behavior of asphalt concrete paving mixtures incor- 
porating the asphalts. Using an indirect method of stiffness determi- 
nation Young et al. (1969) found that the crack frequency of particular 
asphaltic concrete pavements could be related to the low temperature 
stiffness moduli of the asphalts comprising the mixtures. In both of 
these studies, asphalts which exhibited highest stiffness moduli at 
low temperatures were associated with pavements which exhibit highest 


transverse crack frequencies. 


At the present time, direct means of measuring the low tempera- 
ture stiffness of asphalt in a simple and repeatable manner are not 
available. Therefore, low temperature asphalt specification require- 


ments, based on a stiffness approach, have not yet been developed. 


a eee Asphaltic Concrete Fracture Characteristics 


Visual observations of asphaltic concrete cores taken at 
transverse crack locations have indicated that the majority of 
fracture takes place through the asphalt phase. Shields (1966) re- 
ported that within three separate pavements approximately 95 per cent 
of the fracture was through the asphalt phase. He further suggested 
that, because of the soft stone used in the asphaltic concretes, the 


observed stone fracture may have occurred during compaction procedures. 
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Puzinauskas (1966) found that asphaltic concrete cores taken from some 
Ontario pavements exhibited fractured aggregate particles, although to 
a variable degree. Burgess (1966) reported having observed varying 


degrees of large stone fracture where transverse cracks had occurred. 


From observations of asphaltic concrete specimens subjected to 
direct tension tests at temperatures ranging from 30 F to -30 F, Hajek 
(1971) reported that the percentage of aggregate fracture decreased 
with increasing test temperature. He stated that these observations, 
together with those reported by Puzinauskas (1966), lend support to the 
initiation of transverse cracks at low temperatures because it would be 
unlikely that such cracking would occur through the aggregate at high 


temperatures. 


From. the résults of detailed field inventories and test 
projects various investigators, such as Anderson et al. (1966), Shields 
et al. (1969) and Young et al. (1969), have shown that no apparent 
relationship exists between individual asphaltic concrete mix variables, 
such as asphalt content and aggregate gradation, and crack frequency. 
However, Christianson (1970) and Anderson and Shields (1971) have 
presented data showing that tensile fracture properties of field cores 
and laboratory prepared asphaltic concrete specimens, which are 
dependent on such variables, can be correlated with low temperature 


cracking. 


From results of tensile splitting tests performed at 0 F on 
various asphaltic concrete specimens Anderson and Shields (1971) found 


that the tensile failure strain capability of the specimens increased 
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with increasing viscosity at 140 F for a given penetration grade 
asphalt, and increased with increasing asphalt penetration value at 

77 F for a given viscosity. The authors suggested that these findings 
reflect the experience gained with high and low viscosity asphalt 
cements, but is not readily deduced from conventional asphalt tests at 
higher test temperatures. The authors also indicated that although a 
definite limiting failure strain cannot be stated, previous correlations 
with field behavior suggest that values less than 10 x 10°* in per in 
are generally associated with paving mixtures having high potential for 


cracking under climatic conditions of Western Canada. 


Based on results of various test projects, McLeod (1969) 
prepared a tentative design guide for asphaltic concretes with re- 
spect to low temperature cracking. Within this guide, levels of 
peraitic concrete stiffness at which cracking is expected and levels 
at which cracking should be eliminated are given for various minimum 
expected temperatures. The guide, as more fully described in CHAPTER 
VIII, is based on asphalt cement stiffness moduli determined by indirect 
means and dense, well-graded asphaltic concrete mixtures. Such a design 
approach to the problem reflects the findings of field inventories, 
full-scale test sections and laboratory investigations that have shown 
that the major variable contributing to the low temperature fracture 


susceptibility of an asphaltic concrete paving mixture is the asphalt 


used. 
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2.4.5. Asphaltic Concrete Age 


Results of transverse crack frequency surveys involving more 
than 1900 miles of asphaltic concrete pavements in Alberta, summarized 
by Kathol (1968), have indicated a general increase in crack frequency 
with increasing pavement age. Kingham (1966) suggested a similar 
relationship for pavements in the Province of Ontario. Although 
regional transverse crack surveys have suggested a relationship between 
crack frequency and surface age, surveys of specific field projects, 
such as those described by Shields et al. (1969) and Young et al. 
(1969), suggest that climatic conditions, asphalt source, asphalt grade 
and subgrade soil type are more important factors influencing initial 


as well as subsequent crack frequencies. 


Studies of changes in asphalt properties with time as related 
to transverse crack development have been reported by Anderson and 
Shields (1971) and Culley (1969). From a comparison of the in-service 
performance of three sources of similar penetration grade asphalts, but 
ranging from low to high viscosity materials at 140 F, the former 
authors found that the low viscosity source exhibited the largest 
variation in conventional physical properties, most rapid densification 
under traffic and the earliest and largest amount of low temperature 
cracking. Culley (1969) found that changes in asphalt properties, 
within five different projects and after twelve months of service, 


showed no definite relationship with the extent of transverse cracking 


within the projects. 
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From the results of these investigations it would appear that 
a general increase in crack frequency with age may partially be due to 
an increase in asphaltic concrete stiffness caused by densification 
under traffic and by an increase in the probability of occurrence of 
low temperature climatic condition conducive to transverse San 


which increases with pavement age. 


PHL WE Oe Asphaltic Concrete Thickness 


Young et al. (1969) have presented data showing that, when a 
paving mixture susceptible to low temperature cracking is used, crack 
frequency is dependent on the structural thickness of the asphaltic 
concrete. At the Manitoba test site, described by these authors, a 10 
inch full depth pavement section was found to exhibit approximately 
one-half the crack frequency as that associated with a 4 inch asphaltic 


concrete surface of identical material properties. 


Hajek (1971) suggested that such differences in pavement be- 
havior may be explained by the relatively good insulating properties of 
asphaltic concrete and the resulting influence of these properties on 
thermally induced stresses within asphaltic concrete surfaces of 
various thicknesses. In a theoretical study of depth and spacing of 
tension cracks in semi-infinite clastic solids, Lachenbruch (1961) has 
presented data showing that the zone of stress relief in such a material 
increases as crack depth increases. Such data may also aid in ex- 


plaining the dependency of crack frequency on asphaltic concrete 


thickness. 
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2.4.7. Subsurface Materials 

In studies related to low temperature cracking of asphaltic 
concrete pavements the role of the subsurface component layers on 
promoting or retarding cracking has largely been ignored. However, 
from observations and data obtained from field test projects and 
inventory sections various investigators, such as Fromm (1966), 
Anderson et al. (1966), McLeod (1968) and Young et al? (1969)3ahave 
recognized that the subgrade can have a significant effect in some 
cracking situations. The most apparent of these situations is 
evidenced by pavement structures having sand subgrades generally ex- 
hibiting considerably higher crack frequencies than pavements placed 
over clay subgrades, all other variables being the same. Compari- 
sons between such pavement structures have been reported by Young 
et al. (1969), where, the transverse crack frequencies of pavement 
structures having sand subgrades were found to be approximately three 


times the crack frequencies of pavements with clay subgrades. 


The physical properties of granular bases and subgrades of 
pavements, representing twenty individual inventory sections in the 
Province of Alberta, have been included in the study presented by 
Anderson et al. (1966). The authors reported that no explicit rela- 
tionship could be detected between the properties of the subsurface 
materials and the crack frequencies within the control sections. 

Haas et al. (1970) indicated the need for the development of techniques 
that are quantitatively able to relate soil parameters to cracking 


susceptibility and crack frequency. 
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2.9. Summary 

In this chapter a description of low temperature transverse 
cracks, a brief review of performance and economic considerations re- 
lated to the transverse cracking problem, and a summary of some of the 
factors which have been found by various agencies and individuals to be 
associated with the occurrence of low temperature cracking have been 
presented. It should be emphasized that reference has been made to 
only a limited number of published references and a more comprehensive 


review of the subject area has been presented by Haas et al. (1970). 


From: the studies reported, the low temperature fracture 
susceptibility of an asphaltic concrete paving mixture appears 
primarily dependent on the asphalt used. This dependency is reflected 
by differences in crack frequencies of pavements composed of asphalts 
obtained from different crude sources, incorporating asphalts of 
various grades as specified in terms of viscosity and/or penetration 
and including asphalts which exhibit different low temperature 
stiffness properties. Once transverse cracking of an asphaltic concrete 
has initiated the extent of cracking appears to be not only dependent on 


the asphalt but also on such variables as asphaltic concrete thickness 


and subsurface materials. 
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FIGURE II-4 
CORRELATION BETWEEN VISCOSITY AT I40F. AND PENETRATION 
AT 77F FOR CURRENTLY USED ASPHALT CEMENTS (after McLeod, 1967) 


CHAPTER III 


TEMPERATURES IN ASPHALTIC CONCRETE 
PAVEMENT STRUCTURES SUBJECTED TO 


LOW TEMPERATURE CLIMATIC ENVIRONMENTS 


St tepelntroduction 

The importance of temperature as related to the behavior of 
asphaltic concrete paving mixtures has long been recognized. At high 
temperatures asphaltic concrete paving mixtures require sufficient 
Stability to resist deformations resulting from various rates and modes 
of traffic loading, while pavement distress resulting from fatigue due 
to repetitive loadings may be encountered over a wide range of tempera- 
tures. In areas of low temperature climatic environment the fracture 
strength of asphaltic concrete must be sufficient to support traffic 
during spring thaw periods. Stability, fatigue and fracture strength 
considerations have become of increased importance as highway and air- 
field facilities are being subjected to greater traffic volumes coupled 
with increased allowable loadings. In addition to pavement distress as- 
sociated with traffic loadings, pavement fracture cue to thermal effects 
is a serious form of distress in areas subjected to low temperatures. 
The response of a pavement system to these traffic and environmental 


conditions, as measured in terms of deflection, stress and strain, can 
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34 
only be expressed in a quantifiable manner if the temperature depen- 


dence of the asphaltic concrete pavement component is defined. 


Two field studies in Western Canada have yielded continuous 
surface and subsurface temperature records in several asphaltic concrete 
pavement structures subjected to climatic temperatures approaching a low 
of -40 F. In this chapter, the results of statistical analyses of 
temperatures within the pavement structures, during periods of low 


temperature climatic conditions, are presented. 


Se ae Description of Test Sites 


Detailed descriptions and preliminary results of the two field 
test projects, located in Alberta and Manitoba, have been presented by 
Shields et al. (1969) and Young et al. (1969), respectively. Tempera- 
ture as well as other data made available by the agencies involved in 
hese projects has provided the necessary information to carry out the 


present study. 


Sagult rl AlbertasTesterrojecte-"Shields et al. (1969) 

The Alberta test site is located approximately 60 miles south 
of Edmonton and forms a portion of a four-lane divided freeway in 
Central Alberta. In order to evaluate more precisely the possible 
environmental conditions contributing to low temperature transverse 
cracking, and the influence of asphalt source, the Provincial Department 
of Highways in cooperation with the Research Council of Alberta decided 
to incorporate within one 13 mile contract three different sources of 


200-300 penetration grade asphalt cement. These three asphalt sources, 
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obtained from the major refiners in the province, represented high, 
intermediate and low viscosity materials, as measured at 140 F. It was 
felt that within a single paving contract, with uniform subgrade and 
base materials, it would be possible to obtain reasonably homogeneous, 
contiguous surfaces, wherein the relative behavior BEeie asphalt 
cements could be evaluated. By careful planning it was possible to 
change the asphalt supply so that one source was present on one roadway 
and the other two sources abutted at a common point on the other 


roadway. 


The geometrics of the pavement system at the site include two 
40 foot wide paved surfaces separated by a 64 foot depressed median. 
The pavement structure, shown in FIGURE III-1(A), consists of a 12 
inch granular base course surfaced with a two inch asphalt bound base 
and a four inch asphaltic concrete surface. The asphalt bound base and 
granular base was placed in the summer of 1965 with the surfacing the 
following year. The subsoils are low plasticity tills and alluvial- 


lacustrine deposits overlying till. 


During the first two winters of service, measurements of air, 
surface and subsurface temperatures were obtained by means of a 
continuously recording multi-position thermograph used in conjunction 
with thermistors placed at various depths in the pavement structure. 
Also, during this period, continuity circuits, consisting of a high 
conductivity silver base paint, placed in a longitudinal sawn cut on 
the paved shoulder, enabled temperature conditions conducive to thermal 


fracture of the pavement sections to be defined. These conditions have 
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been briefly reviewed in Section 2.2.1 of the previous chapter. 


3.2.2. Manitoba Test Project - Young Clr dine (907) 


The Manitoba test project, more commonly referred to as the 
Ste. Anne Test Road, is a joint research project of the Manitoba 
Department of Transportation and Shell Canada Limited. The test road 
was constructed during the summer of 1967 and constitutes a portion of 
the westbound two lanes of the four-lane divided Trans-Canada Highway 


approximately 25 miles east of Winnipeg. 


The project consists of twenty-nine 400 foot pavement sections, 
each having a 24 foot paved surface width, constructed on highly 
plastic clay and sand subgrades. The primary variables included in the 
test project are; three pavement structures, four different asphalt 
binders, two asphalt contents and two aggregate gradations. The three 
structures are shown in FIGURE III-1 (B), (C) and (D). The asphalts 
used were a low viscosity (LV) 150-200 penetration grade asphalt, a 
high viscosity (HV) 150-200 penetration grade asphalt, a LV 300-400 
penetration grade asphalt and a slow curing, SC-5, liquid asphalt. 
Optimum asphalt contents of approximately five per cent by weight of 
dry aggregate, determined by the Marshall design method, and one per 
cent below optimum, were the main asphalt contents included in the 
test project. For some test sections the gradation of the asphaltic 
concrete aggregate, which was used in the majority of the sections, 


was modified by the addition of portland cement. 
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Continuous measurements of air and pavement temperatures were 
obtained by means of electronic temperature recorders monitoring 
thermocouples placed at various depths through each of the three pave- 
ment structures. In seven of the twenty-nine pavement sections, crack 
detection circuits were used in an attempt to define initial times of 
cracking. Details of the temperature and crack detection systems have 


been described by Deme and Fisher (1968). 


3.2.3. Climatic Environment at Test Sites 

Shields et al. (1969) indicated that air temperatures recorded 
at the Alberta test site, during the winters of 1966-67 and 1967-68, 
could be considered equivalent to long-term average conditions for 
Central Alberta. Long-term daily mean minimum and maximum air temper- 
atures, determined from meteorological records for the Edmonton area, 
are shown in FIGURE III-2. As shown in this figure, daily mean minimum 
and maximum air temperatures are less than 30 F for approximately five 
and three months of the year, respectively. Such temperature data 
emphasizes the need of establishing procedures for the design, testing 
and evaluation of asphaltic concrete pavements at low temperatures, as 


experienced in the Edmonton area. 


Various proposed asphaltic concrete design procedures attempt 
to incorporate climatic environment by defining a critical temperature 
with reference to a particular distress criterion. In the design 
of asphaltic concrete paving mixtures for stability, the critical 
temperature is generally considered an extreme high temperature to 


which the asphaltic concrete is likely to be subjected. The 
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frequencies of occurrence of extreme maximum and minimum monthly air 
temperatures for the Edmonton area, calculated using records for the 
period 1881 to 1968, are shown in FIGURES III-3 and 4, respectively. 
During winter months in the Edmonton area, extreme maximum and minimum 
monthly air temperatures of 30 F and -50 F, respectively, may be 
expected. Fromm and Phang (1971) have suggested that minimum expected 
temperatures for a particular region be chosen on a probabilistic 
basis. The authors have presented iso-temperature contours, for 
Canada, based on a one per cent recurrence interval, from which a 
winter design temperature of -40 F, for Central Alberta, is suggested. 
This temperature is in close agreement with an air temperature of -50 
F, as asphaltic concrete temperatures are generally somewhat higher 


than air temperatures during periods of extreme cold. 


Extreme maximum and minimum air temperatures in the Winnipeg 
area are Similar to those of Edmonton. Bihourly air and pavement 
surface temperature frequencies for Structure B, identified in Section 
3.2.2, for the time period from November 21, 1967 to March 31, 1968, are 
shown in FIGURE III-5. During this time, air and surface temperatures 
were less than 30 F for approximately 85 per cent of the time and a 


minimum air temperature of -45 F was recorded. 
For ease of presentation, the one Alberta structure together 


with the three Manitoba structures, as shown in FIGURE III-1, will be 


referred to as Structures A, B, C and D, respectively. 
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Ce Asphaltic Concrete Temperatures 


For each of the four structures, temperature data included in 
the analyses were for the time periods given in TABLE III-1. Malfunc- 
tions of the recording apparatus resulted in incomplete temperature 
data at the two inch depth of Structure B and only those days with 
complete asphaltic concrete temperature records were included in the 


following summaries. 


A summary of the daily mean maximum and minimum air and asphal- 
tic concrete temperatures of each structure is given in TABLE III-2. 
The average temperatures of Structure B were less than those of 
Structures C and D. This difference can mainly be attributed to the 
fact that the time period of complete temperature records for Structure 
B was limited to days of seasonal minimum temperature. However, as 
shown in FIGURE III-6, at equivalent depths and times, asphaltic con- 
crete temperatures of Structure B were less than those of Structure C. 
Young et al. (1969) suggested that such temperature differences may 
partially be accounted for by the fact that Structure B is located in 
an area void of shelter, while Structure C is located in a wooded area 


and sheltered from prevailing winter winds. 


As observed from TABLE) III-2, the daily mean maximum and 
minimum temperatures of Structure C approximated those within the 
upper inches of Structure D. This approximation is made by interpo- 
lating the temperatures of Structure D to depths equivalent to depths 
at which temperatures were measured in the asphaltic concrete surface 


of Structure C. Similar results have been reported by Kallas (1966), 
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who found that temperatures at depths of 2, 4 and 6 inches in 6 and 12 
inch asphaltic concrete surfaces were essentially the same when sub- 


jected to similar climatic environments. 


SeSel: Daily Temperature Variations 


Although average temperatures, such as summarized in TABLE 
III-2, are useful in describing the general climatic conditions to 
which pavements in Western Canada are subjected, during periods of low 
ambient temperatures, such an analysis does not yield specific infor- 
mation as to the daily temperature variations within such pavement 
surfaces. Linear regression analyses, using the method of least 
Squares, were employed to develop relationships between daily air and 
surface temperature ranges and daily ranges of temperature at the 
depths of temperature recordings in the four pavement surfaces. The 
results of these analyses are given in TABLE III-3. The relatively 
high standard errors of estimate of those equations relating daily air-: 
surface temperature variations were expected, since, climatic factors 
other than just air temperature influence pavement temperature. A 


review of some of these factors is given in Section 4.4 of CHAPTER IV. 


The regression equations, shown in TABLE III-3, relating daily 
pavement temperature ranges were equated to daily surface temperature 
variations experienced at the test sites. FIGURE III-7 shows the 
resulting daily pavement temperature variations with respect to depth. 
The daily temperature ranges at the two inch depth of Structures A, B 
and C were approximately 80 per cent of any particular daily pavement 


surface temperature variation. With increased depth, the daily 
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temperature variations decreased, and with the exception of the six 
inch depth shown for Structure A, this decrease appears independent 


of pavement structure. 


te a es Temperature Gradients 
As mentioned in Section 2.2.1 of CHAPTER II, Shields eteal. 


(1969) and Young et al. (1969) reported that transverse cracking at 
the Alberta and Manitoba test sites occurred during times of minimum 
recorded air and pavement temperatures. The time of daily minimum 
air temperatures was chosen as a reference time and further linear 
regression analyses were performed, relating daily minimum air tempera- 
tures and asphaltic concrete temperatures at this particular reference 
point in time. The resultant relationships are shown in TABLE III-4. 
The standard errors of estimate of equations relating daily minimum 
air temperature to pavement surface temperature are approximately one- 
half those relating daily temperature ranges, as previously presented 
in TABLE III-3. The improved correlations may be attributed to the 
fact that at the test sites minimum daily temperatures were generally 
recorded during early morning hours at which time shortwave radiation 


does not have a direct influence on pavement temperatures. 


The regression equations given in TABLE III-4 were equated to 
daily minimum air temperatures experienced at the test sites and the 
resulting minimum air-pavement temperature relationships are shown in 
FIGURE III-8. These relationships indicate that the temperature 


gradients increased with a decrease in daily minimum air temperature. 
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Such an increase of temperature gradient is especially noted for 
Structure D, where, at a minimum daily air temperature of 0 F a 
temperature differential between the surface and 10 inch depth of 15 F 
could be expected while, at -35 F the corresponding temperature dif- 
ferential is increased to approximately 25 F. Such changes in 
temperature gradients with respect to daily minimum air temperatures 
can be attributed to the particular climatic conditions to which the 
pavement structures were subjected. Examination of the temperature 
records revealed that the structures were not subjected to prolonged 
periods of extreme low temperatures. Therefore, even though the 
minimum recorded pavement surface temperature of Structure D, during 
the 1967-68 winter, was -37 F, the minimum recorded 10 inch depth 


temperature was -6 F. 


yoyo. Races OF Temperature Change 


The decrease of daily temperature variations with respect to 
depth, as shown in FIGURE III-7, implies a decrease in the time rate 
of temperature change with increased depth. The time rate of tempera- 
ture change corresponds to a strain rate, which is directly proportion- 
al to the temperature change provided that the thermal coefficient of 
expansion of asphaltic concrete is assumed temperature independent. 
Various investigators, Domaschuk et al. (1964), Monismith et al. (1965) 
and Littlefield et al. (1967) have determined the expansion coefficient 
of asphaltic concretes over various temperature ranges. The reported 


coefficients vary from approximately 1.2 x 107° to 1.7 x 10°5 per degree 
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Fahrenheit. Recently, Burgess etrale) (1971) teported. test results 
on asphaltic concrete specimens simulating those mixtures incorporated 
in the Manitoba test project. The average coefficient of expansion of 
the specimens over a temperature range from 68 to -40 F was approxi- 
mately 1.1 x 107° per degree Fahrenheit. The results of these studies 
Suggest that, for practical purposes, the coefficient of thermal 
expansion of a given asphaltic concrete may be assumed temperature 
independent and, therefore, a strain rate may be considered directly 


proportional to a time rate of temperature change. 


The rates of temperature change which occurred in the asphaltic 
concretes of Structures C and D are shown in FIGURES III-9 and III-10, 
respectively. During the 1967-68 winter season, the rate of surface 
temperature change of each structure was equal to or less than +3 
degrees Fahrenheit per hour for approximately 90 per cent of the time, 
and equal to or less than -3 degrees Fahrenheit per hour for 10 per 
cent of the time. At depths of four inches the corresponding rates of 
temperature change were reduced by a factor of three and at the 10 inch 
depth of Structure D the rates of temperature change were greatly 
reduced and were generally within * 0.5 degrees Fahrenheit per hour. 


Maximum time rates of surface temperature change recorded at the 


Manitoba test site were in the order of 8 to 10 degrees Fahrenheit per 
hour. 
A knowledge of such time rates of temperature change is of 


great value in establishing realistic and meaningful laboratory 


testing environments. Investigators, such as Hills and Brien (1966), 
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Haas (1968) and Fromm and Phang (1971), have assumed various time 
rates of temperature, change in studies related to the low temperature 
fracture susceptibility of asphaltic concrete paving mixtures. Hills 
and Brien (1966) subjected restrained asphaltic concrete specimens to 
a cooling rate of 10 degrees Centigrade per hour, while, Haas (1968) 
used a cooling rate of 30 degrees Fahrenheit per hour, for determining 
thermally induced stress-temperature relationships of asphaltic con- 
crete paving mixtures. Fromm and Phang (1971) proposed a method of 
determining the relative low temperature fracture susceptibility of 
different paving mixtures, in which a rate of 10 degrees Fahrenheit per 
hourenshutalizeds This value was determined from an examination of 
hourly temperature records from various locations within the Province 
of Ontario which indicated that the ambient temperature drop during a 
one hour period does not exceed 10 degrees Fahrenheit per hour at 


temperatures below 0 F. 


The results of the present study suggest that, for the climatic 
environment experienced at the Manitoba test site, time rates of 
asphaltic concrete temperature change are generally less than those 
assumed in the preceding laboratory studies, while, a maximum cooling 
rate of approximately 10 degrees Fahrenheit per hour is consistent 


with that utilized by Fromm and Phang (1971). 


eaten Subgrade Temperatures 


The thermal regimes which existed in the subgrades of Struc- 


tures B, C and D, during the 1967-68 winter, are shown in FIGURES 
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Iii-11, IIl1l-12 and III-13, respectively. The isotherms shown in these 
figures were determined by linear interpolation between recorded 
temperatures measured at two foot intervals in each of the subgrades. 
Minimum ambient temperatures at the Manitoba test site were measured 
during late December and early January during which time maximum sub- 
grade temperature gradients were recorded. Uncertainty exists as to 
whether these subgrade temperature gradients had a direct influence on 
initial transverse cracking of the asphaltic concrete pavement 
Surfaces. However, aS mentioned in Section 2.4.5 of CHAPTER IT, 
transverse crack frequencies at the Manitoba test site were found to 
be dependent on subgrade soil type, and since the strength and volume 
change characteristics of soil are temperature dependent, subgrade 
thermal regimes may have contributed to crack initiation as well as 


subsequent crack development. 


The maximum depth of frost penetration, as calculated from the 
position of the 32 F isotherm below the subgrade surface, within the 
clay subgrade of Structures B and D was 70 inches, while, the maximum 
depth of frost penetration in the sand subgrade of Structure C was 110 
inches. Such a difference in frost penetration depth reflects the 
influence of the different thermal properties and latent heat of water 
of the two soil types. The thermal properties referred to’ are dis= 


cussed in Section 4.6 of the following chapter. 


The bihourly 18 inch and 10 inch depth temperature frequencies 
of Structures B and D, respectively, are shown in FIGURE III-14. The 


similarity of these temperature frequencies, together with the equal 
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depths of frost penetration within the two structures, suggests that, 
for the prolonged subfreezing subgrade temperatures experienced at the 
Manitoba test site, the 10 inch asphaltic concrete pavement surface of 
Structure D and the four inch surface plus 16 inch granular base 
course of Structure B may be considered equivalent with respect to 


their influence on subgrade temperatures. 


3.5. Summary 

in this cNapeer recorded temperatures in various asphaltic 
concrete pavement structures subjected to low temperature climatic 
environments experienced in Western Canada have been summarized. 
Relationships between daily air and asphaltic concrete temperature 
variations, and between daily minimum air temperatures and asphaltic 
concrete temperatures have been developed using linear regression 
techniques. Time rates of temperature change within two asphaltic 
concrete pavement surfaces and thermal regimes which existed in sub- 
grades of three pavement structures, subjected to near identical 


climatic environments, have been presented. 
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TABLE III-1 


TIME PERIOD OF TEMPERATURE DATA 


INCLUDED IN THE ANALYSES 


A 
B 
C 
D 


October 1, 1966 - March 31, 1967, 
January 1, 1968 - February 29, 1968 


November 15, 1967 - February 29, 1968 
November 15, 1967 - March 31, 1968 
November 21, 1967 - March 31, 1968 
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TABLE III-2 


SUMMARY OF AIR AND ASPHALTIC 


CONCRETE TEMPERATURES 


Minimum 
Temperature 
) 


Daily Mean 


Maximum 
Temperature 


Daily Mean 
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Figures in brackets indicate one standard deviation. 
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TABLE III-3 


RELATIONSHIPS BETWEEN AIR AND ASPHALTIC 
CONCRETE DAILY TEMPERATURE RANGES 


Linear Regression 


Coefficient Standard 
of ERVonof 


) * 
Equation Correlation Estimate (F) 


1.30 + .79 
0.90 + .8]1 
0.75 + .42 


he 


Daily range in air temperature (F) 


Daily range in surface temperature (F) 


Ra? etc. 


Daily range in 2 and 4 inch asphaltic concrete 
temperatures (F) 
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TABLE III-4 


RELATIONSHIPS BETWEEN AIR AND ASPHALTIC 
CONCRETE TEMPERATURES AT A TIME OF 
MINIMUM DAILY AIR TEMPERATURE 


Coefficient Standard 
of Error of 
Correlation Estimate (F) 


Linear Regression 
Equation* 


Minimum daily air temperature (F) 


Pavement surface temperature (F) 


etc. 


= Asphaltic concrete temperature at 2 and 4 inch 
depths (F) 
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DAILY MEAN 
MAXIMUM AIR 
TEMPERATURE. 


DAILY MEAN 
MINIMUM AIR 
TEMPERATURE 


Temperature (°F) 


From temperature 
records dating from 
I88!l to 1968 


Months 


FIGURE IITI-2 DAILY MEAN MAXIMUM AND MINIMUM 
AIR TEMPERATURE FOR EDMONTON, ALBERTA. 
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FIGURE IIT-7 DAILY RANGE OF TEMPERATURES AT VARIOUS DEPTHS 
IN ASPHALTIC CONCRETE PAVEMENT SURFACES VERSUS 
DAILY RANGE IN PAVEMENT SURFACE TEMPERATURE. 
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FIGURE III-8ASPHALTIC CONCRETE PAVEMENT TEMPERATURES 


AT VARIOUS DEPTHS VERSUS MINIMUM DAILY 
AIR TEMPERATURE. 
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CHAPTER IV 
HEAT TRANSFER IN PAVEMENT STRUCTURES 


4.1. Introduction 

In the absence of recorded pavement temperatures, such as 
summarized in the previous chapter, a realistic assessment of the 
response of asphaltic concrete pavements to various climatic and 
loading conditions is dependent on a method whereby pavement tempera- 
tures can be ‘reliably predicted. The application of principles of 
heat transfer in which the influence of climatic factors and pavement 
material properties on pavement temperatures are considered enable 
such predictive methods to be developed. In a study of asphaltic 
concrete fracture temperatures constitute a primary input variable for 
thermal stress computations and tensile strength determinations. In 
this chapter heat transfer concepts applied to pavement structures and 
thermal properties of pavement component materials are reviewed. These 
concepts and properties have been employed in CHAPTER V where an 


analytical model for predicting pavement temperatures is developed. 


ALi Governing Heat Transfer Equation 


The one-dimensional heat conduction equation governing the 


dissipation of heat in a homogeneous body is expressed as follows: 
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2 4 
a nee (IV-1) 
ax? at 
where: a = diffusivity of the body, 


J 
u 


the temperature of the body at a specific 
depth=x™and time “t 

In a pavement structure the assumption of one-dimensional heat 
transfer is most likely approximated near the pavement surface where 
the depth to width ratio of the structure is small. Other simplifying 
assumptions, concerning negligible moisture and volume changes due to 
temperature variations, are also inherent in the governing heat con- 
duction equation. As mentioned in Sections 2.2.2 and 2.2.5 of CHAPTER 
II, substantial volume changes are known to occur in certain subgrade 
soil types and such changes are often found to be most prevalent in 
subgrades subjected to frost action. When solving Equation (IV-1) for 
the pavement temperature regimes errors introduced by assuming constant 
volume and moisture conditions are difficult to quantify. However, such 
errors are probably small in comparison to the errors introduced when 
attempting to define the thermal properties of the component layers and 


the influence of various climatic variables on pavement surface 


temperature. 


The solution of Equation (IV-1) is dependent on prescribed ini- 
tial and boundary conditions. In heat conduction problems the tempera- 
ture distribution throughout a body is generally stipulated as an 
initial condition, and in relation to a pavement structure, this 


requires that the temperature gradient at a particular time be known or 
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estimated. The boundary conditions describe the influence of the 
surroundings of a body on its surface. At extreme depths in a pavement 
structure changes in temperature with respect to depth and time become 
small as compared to temperature variations near the pavement surface. 
Therefore, the assumption of a constant temperature at a specified 

depth and for a given time period constitutes a possible boundary condi- 
tion when determining the temperature distribution in a pavement 
structure as a function of time. A relationship between various clima- 
tological factors and pavement surface temperature forms a second 
required boundary condition when solving Equation (IV-1) for the thermal 


regime in a pavement system. 


4.3. Ambient and Pavement Surface Temperatures 


The most readily recognized climatic variable influencing the 
thermal regime of a pavement structure is atmospheric temperature. 
Methods of predicting frost penetration depths in pavement structures, 
such as those presented by Aldrich (1956) and Sanger (1963), express 
atmospheric temperature in terms of the parameters degree-days and 
freezing index. These parameters, together with a knowledge of the 
thermal properties of the component layers, permit approximate depths 


of frost penetration to be calculated. 


The degree-days of freezing for any one day equals the difference 
between the average daily air temperature and a reference temperature, 
provided that the average daily air temperature is less than the 
reference temperature. The reference temperature is generally assumed 


equal to 32 It. The freezing index is defined as the difference between 
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the maximum and minimum points on a cumulative degree-day plot and is a 
measure of the combined duration and magnitude of below freezing temper- 
atures. The freezing index can be calculated on the basis of air or 
pavement surface temperature. However, since air temperature data is 
generally more readily available, the air freezing index is commonly 
computed and converted to a surface freezing index by assuming a con- 


stant ratio between average air and pavement surface temperatures. 


Although air and pavement surface temperatures are related, 
various investigators, Carlson and Kersten (1953), Kersten and Johnson 
(1955), Sanger (1963) and Penner et al. (1966) have found that this 
relationship is dependent on such factors as geographic location, 


thermal properties of the pavement materials and surface type. 


Carlson and Kersten (1953) suggested a value of 0.6 as the 
Tatio of surface freezine index to air freezing index, componly termed 
the "n-factor", for asphaltic concrete, portland cement concrete and 
gravel surfaces in Alaska. Kersten and Johnson (1955) concluded that 
a n-factor of 0.80 was applicable for asphaltic concrete pavement sur- 
faces in the Minneapolis, Minnesota area. Sanger (1963) reviewed the 
results of frost penetration studies performed by the United States 
Corps of Engineers and concluded that, for climatic conditions ex- 
perienced in the northern states, the surface freezing index for bare 
pavements and smooth soil surfaces is approximately 90 per cent of the 
air freezing index. Penner et al. (1966) found n-factors of 0.68 and 
0.58 for foam plastic insulated and uninsulated asphaltic concrete 


pavement structures, respectively, at a test site in Sudbury, Ontario. 
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Methods of predicting depths of frost penetration that incorpo- 
rate the degree-day concept have the disadvantage of yielding only total 
frost penetration depths and not temperature variations in a pavement as 
a function of time and position. Coupled with this factor are the 
difficulties involved in converting an air freezing index to a surface 
freezing index. As observed from the relatively wide range of reported 
n-factors, the relationship between air and pavement surface tempera- 
tures is complex. A rational approach to the development of such a 
relationship involves defining heat transfer phenomena occurring at an 


air-pavement surface boundary. 


4.4. Heat Transfer at an Air-Pavement Surface Boundary 

Thermal energy is transferred to and from a pavement surface by 
means of radiation, conduction and convection. Aldrich (1956) suggested 
that the transfer of thermal energy between the air and a pavement 
surface is most affected by direct and diffuse solar radiation, net 
longwave radiation between the pavement and atmosphere, and convection. 
The interacting radiation and convection heat transfer modes, as pre- 


sented by Aldrich (1956), are shown in FIGURE IV-1. 


The influence of individual climatic variables such as solar 
radiation, wind velocity and cloud cover, on the heat cconomy of a 
pavement surface is difficult to quantify because of their variation 
in intensity, time and geographic location. Therefore, predictive 
methods of defining the influence of such climatic variables on pavement 
surface temperatures rely heavily on empirical relationships determined 


from laboratory and field investigations. Also, the availability of 
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meteorological data at a particular location often dictates the choice 
of methods used to predict the radiation, conduction and convection heat 


"quantities". 


4.4.1. Shortwave Radiation 

Shortwave radiation includes the complete spectrum of solar 
radiation and is the principal source of heat energy received by a pave- 
ment surface during daylight hours. Only a portion of the solar 
radiation received by the earth's atmosphere reaches the surface of the 
earth. The intensity of direct and diffuse incident shortwave radiation 
is dependent on such factors as latitude, altitude, cloud cover, time of 


year and the time of day at the particular observation site. 


Major meteorological stations record the intensity of direct and 
diffused solar radiation in terms of Btu per sq ft per hr or langleys 
per sq ft per hr. Where such data is not available the shortwave 
radiation heat flux at a pavement surface can be estimated by computing 
the theoretical shortwave radiation received at the site for a clear 
day and applying a reduction factor based on the per cent of possible 


daily sunshine received at the site. 


Dempsey and Thompson (1970) observed that the intensity of 
shortwave radiation received by the earth's surface varies parabolically 
from the time of sunrise to the time of sunset. Based on this obser- 
vation, the authors developed a computer program from which the inten- 
sity of shortwave radiation on pavement surfaces, during finite time 


intervals, can be readily predicted. 
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Iimploying a system of solar radiation cells, Straub et yal,. 
(1968) measured the incident solar radiation on an asphaltic concrete 
pavement surface located at Potsdam, New York. From a comparison of 
measured solar radiation values and recorded air and pavement tempera- 
tures, Straub et al. (1968) concluded that for cloudless days and 
climatic conditions experienced at the test site, variations in solar 
radiation have a greater effect on asphaltic concrete pavement tempera- 


tures than caused by changes in air temperature. 


A semiempirical method of simulating maximum pavement tempera- 
tures was developed by Barber (1957). This method includes the effect 
of solar radiation on pavement temperatures by defining an effective 


air temperature as follows: 


Lats" ephartaat/h (IV-2) 
where: T, = effective air temperature, F 

T, = air temperature, F 

a = absorptivity of the pavement surface 

I =e solar tadlation ,ABtu-per sq: tte per hr 

h = surface coefficient, Btuvper sq ft per hriper, FF. 


The absorptivity, a, of a pavement surface equals the ratio of 
solar radiation absorbed to that received by the surface and the 
distribution of incident energy is independent of the temperature or 
physical nature of the surface. In previous pavement temperature 


studies absorptivity values ranging from 0.35 to 1.00 have been 


assumed. 
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The shortwave radiation heat flux influencing temperatures at 


an air-pavement surface boundary may be expressed as follows: 


(IV-3) 


where: | Q = shortwave radiation heat flux: contributing 
to the thermal regime of a pavement 


structure;) Btue permsqrftiiper hr 
a = absorptivity of the pavement surface 


oe = solar energy incident on the pavement 


surtace, Btusper sq it per. hr. 


AO 4.2. Longwave Radiation 


Bodies emit energy at a rate proportional to the fourth power of 
their absolute temperature and the radiation emitted a given wavelength 
is found by Planck's spectral distribution. The transfer of the thermal 
energy by means of longwave radiation at an air-pavement surface boundary 
is most evident during cold cloudless nights. As shown in FIGURE IV-1, 
during the day a pavement surface is warmed by shortwave radiation. 
However, at night the shortwave radiation incident on a pavement surface 
is reduced to zero and thermal energy transferred by longwave radiation 


from the pavement results in a rapid cooling of the pavement surface. 


Straub et al. (1968) approximated the longwave radiation heat 


flux at an air-pavement surface boundary by means of the following 


relationship. 
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where: ors = the heat flux resulting thermal energy 

transferred by longwave radiation, 
Btu per sa ft per,hr 

e = emissivity of the pavement surface 

oO = Stefan-Boltzmann constant, 
0.1714. xn Os? Btu per Sdeie peribr pers), 

T; = pavement surface temperature, R 

to -thaltotemperature, R, 


The emissivity of a surface equals the ratio of emitted radiant 
heat flux per unit area to that of a blackbody at the same temperature 
and is dependent on such variables as temperature, wavelength and sur- 
face composition. Previous asphaltic concrete pavement temperature 


Simulation studies have incorporated emissivity values ranging from 


Grlsato ks00: 


Scott (1957) indicated that longwave radiation at the earth's 
surface is not only dependent on air and surface temperatures but also 
on cloud cover and vapour pressure. Meteorological stations generally 
record the duration of bright sunshine during a day which is determined 
by measuring cloud cover duration. Therefore, when applying a cloud 
cover correction to the longwave radiation entering the energy balance 
at an air-pavement surface boundary, it is necessary to assume that the 


measured cloud cover duration is also representative of that which 


occurs during nighttime hours. 


Geiger (1959) developed an empirical equation in which the 


longwave radiation emitted by the atmosphere is expressed as a function 
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of air temperature and vapour pressure. This relationship is as 


follows: 


Oia oo] ots? Tete to Py) (IV-5) 


In Equation (Ves Ope is the heat flux resulting from longwave 
radiation emitted by the atmosphere with no cloud cover correction. The 
terms o and i are as defined in Equation (IV-4) and Geiger assigned the 
following values to the constants G,; J, and 0: G = 0.77, J = 0728 and 

p = 0.074. Geiger (1959) suggested that the vapour pressure, p, varies 
from approximately 1 to 10 millimeters of mercury for climates near the 


earth's surface. 


4.4.3. Conduction 

The existence of a temperature gradient between a pavement 
surface and the layer of air immediately above the surface results in a 
transfer of thermal energy from the high temperature medium to the low 
temperature medium. The heat transfer mode is that of conduction and 
the rate of heat transfer is proportional to the normal temperature 


gradient at the air-pavement surface boundary. The resulting heat flux 


can be expressed as follows: 


dT 
= = = We 
Cond . dx ( 6) 


where: Peer = the heat flux resulting from heatetranster 


by conduction, Btu per sq ft per hr 


K = thermal conductivity, Btu per hr per ft per F 


dt/dx temperature gradient, F/ft. 
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Conduction is the primary mode of heat transfer in a pavement 
Structure and the energy balance of other heat "quantities", inter- 
acting at a pavement surface, determines the quantity of thermal energy 


remaining for conductive heat transfer within a pavement. 


Geiger (1959) and Barry (1968) suggested that since air is a 
poor conductor heat transfer by conduction in the atmosphere can be 
virtually neglected. The authors also suggested that the low viscosity 
of air and its consequent ease of motion makes convection the primary 


method of atmospheric heat transfer. 


4.4.4. Convection 

The temperature gradient at an air-pavement surface boundary is 
influenced by the mixing motion of the air immediately above the surface. 
Convective heat transfer at the boundary is a process of enero trans- 
port by the combined action of heat conduction and the mixing motion of 
the air and, therefore, a pseudo heat conduction process. Newton's 
Law of Cooling states that the heat flux across the boundary of a 
surface, being cooled by forced convection, is proportional to the 
temperature difference between the surface and the surrounding medium. 


This heat flux can be expressed as follows: 


Q = HAT (IV-7) 


where: Q = heat flux resulting from convective heat 
transter, Btu per sq fe pereir 
I] = surface Convection coef ficaent, 
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AT = temperature difference between the surface 
and the surrounding medium, F. 

The surface convection coefficient, H, which depends upon such 
factors as wind velocity, surface type and surface roughness, is 
difficult to evaluate. Carroll et al. (1966) suggested a value of 2.7 
Btu per sq ft per hr per F for granular surfaces and wind velocities 
ranging from 0 to 10 miles per hour. Vehrencamp (1953) developed the 
following empirical formula for estimating the convection coefficient 


at an air-soil boundary. 


Hi= 08001445. 0:2 Us 7) +20cN0097 (Keaoda) ay (IV-8) 


where: H = surface convection coefficient, 


CaleperYsqtcm*per min” per’ C 


Tm = average of the air and surface temperatures, K 
U = average daily wind velocity, m per sec 

Ts = surface temperature, C 

Tg = air temperature, C. 


The data used to formulate Equation (IV-8) was collected from a 
study of the air temperature distribution close to the surface of a 
smooth, level, dry lake bed. Vehrencamp (1953) indicated that Equation 


(IV-8) yielded results within 10 per cent of experimental data. 


The temperature of a pavement surface is also dependent on such 
phenomena as precipitation and evaporation. lIlowever, the time of 
occurrence and the influence of such hydrologic variables on pavement 
temperature cannot be readily defined. Therefore, heat fluxes resulting 


from such variables have not been considered in previously mentioned 
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asphaltic concrete temperature simulation studies. Sanger (1959) con- 
cluded that phase changes of water, at an air-pavement surface boundary, 
could be safely ignored when calculating frost penetration depths be- 
neath highway and airfield pavements. A similar opinion was held by 


nldrich, (1956) . 


4.5. The Energy Balance Approach 


Utilizing an energy balance approach, the algebraic sum of all 
heat "quantities" at an air-pavement surface boundary is equated to 
zero. Neglecting heat transfer resulting from phase changes in water 
and precipitation, and considering heat transfer towards a pavement 
surface as positive and away from the pavement surface as negative, the 
energy balance equation involving the previously defined heat fluxes is 


as follows: 


+ Q = 0 (IV-9) 


+ + 
Cee Fi iN 4 Seon conv 


Equation (IV-9) provides a basis for establishing a realistic 
relationship between many climatic variables and pavement surface 
temperatures. The degree of accuracy of such a relationship for pre- 
dicting pavement surface temperatures is largely dependent on the 
method and parameters used to predict the influence of the various 
climatic variables on the individual heat fluxes. Jurthcermore, 
Liquation (1V-9) furnishes a second possible boundary condition for the 
numerical solution of Equation (IV-1), from which, the thermal regime 
in a pavement structure can be predicted provided that the thermal 


properties of the component layers are defined. 
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4.6. Thermal Properties of Component Layers 


In order to calculate the thermal regime in a pavement system 
a knowledge of two independent thermal properties, namely, thermal 
conductivity and heat capacity of each component layer is required. 
Thermal conductivity, K, expresses the rate of heat flow through a unit 
area under a unit temperature gradient. Most commonly used units of 
thermal conductivity are Btu per hr per ft per F or Cal per sec per cm 
per C. Thermal conductivity is generally evaluated by measuring the 
temperature gradient along a specimen of the material conducting heat at 


a known rate. 


The heat capacity, C, is the amount of thermal energy required 
to cause a one degree change in a unit mass of the material. The units 
of heat capacity depend on whether volumetric or mass heat capacity is 
uséd. The units of volumetric heat capacity are Btu per cu ft per F or 
Cal per cc per C while mass heat capacity units are Btu per 1b per F or 
Cal per g per C. Often specific heat capacity is specified. The specif- 
ic heat of a material is dimensionless and equals the ratio of the heat 
capacity of the material to the heat capacity of water. Generally, the 
heat capacity of a material is computed from a heat balance between the 


heat gained by water in a calorimeter and the heat gained by the 


material. 


The thermal conductivity and heat capacity of a heterogeneous 
material can be approximated by summing the thermal properties of the 
constituents comprising the material. The volume or weight fractions 


of the individual constituents must be considered in such a summation 
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process. The thermal conductivity and heat capacity of pavement 
materials subjected to freezing and thawing conditions are somewhat 
more difficult to define since latent heat associated with these 


conditions must also be evaluated. 


4.6.1. Thermal Conductivity of Asphaltic Concrete 


Kersten (1949) determined the thermal conductivity of an 
asphaltic concrete paving mixture having a density of 138 pcf and an 
asphalt content of 6 per cent by dry weight of aggregate. The thermal 
conductivity was found to range from 0.82 Btu per hr per ft per F at 


#29 § toyOse86,Btuy,perchryper ft pernFiat 40sF. 


Johnson (1952) reported a thermal conductivity of 0.82 Btu per 
hi per ft per FP for an asphaltic concrete mixture Naving” a density or 
150 pef and an asphalt content of 4.5 per cent by dry weight of 


aggregate. 


Aldrich (1956) suggested a value of 0.84 be used for asphaltic 
concrete paving mixtures near freezing temperatures, and Penner et al. 
(1966) assumed this value when calculating frost penetration depths 
within various asphaltic concrete pavement structures. In recent 
asphaltic concrete temperature prediction analyses, such as presented 
by Corlew and Dickson (1968), Kasianchuk (1968) and Straub et al. 


(1968) , a valuc of 0.70 was assumed, a value suggested by Barber (1957). 


., heat Capacity of Asphaltic Concrete 


Due to the relatively large weight fraction of aggregate in 


aS 
= 
— 
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asphaltic concrete mixtures, as compared to the weight fraction of 
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asphalt, the heat capacity of such a mixture is generally assumed to 
approximate the heat capacity of the aggregate comprising the mixture. 
As discussed in Section 4.6.4 in this chapter, the heat capacity of 
dry aggregate has been found to vary from approximately 0.16 to 0.25 
Btu per 1b per F. From an extensive literature review Johnson (1952) 
concluded that a value of 0.20 Btu per 1b per F was representative for 
most dry soils and rocks, and Dempsey and Thompson (1970) suggested 

a heat capacity of between 0.20 and 0.22 Btu per 1b per F as being 


realistic for most asphaltic concrete paving mixtures. 


£503.04 Thermal Conductivity of Granular Bases and Subgrade Soils 

The most widely used method of estimating the thermal conduc- 
tivity and heat capacity of subsurface pavement materials is based on 
the results of an extensive laboratory investigation reported by 
Kersten (1949). Kersten determined the thermal conductivity of 19 
different soils and related the effect of density, moisture content, 
temperature, texture and mineral composition to the resulting thermal 
conductivities. From the test results, Kersten (1949) developed 
equations which express the thermal conductivity of frozen and unfrozen 
soils in terms of soil type, moisture content and dry density. The 


equations are as follows: 


Fine Textured Soils: 


-0.2) 109-01¥d 
Unfrozen K, = ne 5 ) (IV-10) 
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Frozen ee be 0-08 (10) 40 DE C0) (IV=11) 
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Granular Soils: 


0.01Yd 
Unfrozen K, = Silvie ts wre) sd (iv-12) 
0-013Yd 0.146Yd 
Frozen Kr = oO70 EO) + 0.082 C0) (IV-13) 
where: K, = thermal conductivity of unfrozen soil, 
BLusper hr per sft per F 
Ke = thermal conductivity of frozen soil, 
Deupper he per ft. per & 
Ya = “dry density, pcr 
w = moisture content, percent by dry weight of soil. 


Kersten (1949) considered the above equations to be accurate 


* 25 per cent and suggested that soils with 50 per cent or more 


within 
silt and clay sizes, particles less than 0.05 mm, be defined as fine 
textured soils and those with less than 50 per cent silt and clay size 
Parecicres De GeLined as pranular soiis. “Kersten, indicated that 
Equations (IV-10) and (IV-11) are valid for soils having moisture 
contents equal to or greater than 7 per cent, and Equations (IV-12) and 


(IV-13) are valid for soils having moisture contents equal to or 


ereacer thane, per cent. 


As shown by Equations (IV-10) through (IV-13), the thermal 
conductivity of a soil depends upon whether the soil is in a frozen or 
unfrozen condition and difference in thermal conductivity between these 
two conditions is largely dependent on the moisture content of the 
soil. At low moisture contents, the difference between unfrozen and 


frozen thermal conductivity of a given soil is small. However, with 
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increasing moisture contents the thermal conductivity of a frozen soil 
becomes progressively greater than the thermal conductivity of the soil 
in an unfrozen condition. Since the percentage of water frozen in a 
soil-water system is temperature dependent, the thermal conductivity of 
such a system in a freezing condition is also dependent on temperature. 
Sanger (1963) suggested that the thermal conductivity of soil in a 
freezing condition could be approximated by averaging unfrozen and 


frozen thermal conductivity values. 


Kersten (1949) concluded that an increase in dry density of a 
soil, while maintaining a constant moisture content, resulted in an 
increase in thermal conductivity. In a later summary, Kersten (1959), 
suggested that as a general rule a 1 pcf change in dry density will 


result in approximately three per cent change in thermal conductivity. 


The influence mineral composition on the thermal conductivity 
of soils has been studied by such investigators as Kersten (1949) and 
Farauki (1966). Kersten (1949) found that the thermal conductivity of 
a soil increased as quartz content increased. Farouki (1966) determined 
the thermal conductivity of quartz, sand and kaolinite clay mixtures and 
found that the maximum thermal conductivity was obtained at optimum clay 
contents. Kersten (1959) concluded that the lack of experimental re- 
sults makes it difficult to take mineral composition into account when 


estimating the thermal conductivity of a soil. 


4.6.4. Heat Capacity of Granular Bases and Subgrade Soils 


Kersten (1949) determined the heat capacity of various dry soils 


over a temperature range of 140 F. The heat capacity of each soil was 
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and 140 F, respectively. However, base and subgrade materials comprising 


a pavement structure are seldom, if ever, in a moisture free condition. 


Kersten (1949) suggested that t 


estimate the.specific heat of a 


cone 100 Cg + 

ui” 100% 

Cae 100 Cs + 

Ea 100 + 

where? “Cj; =" heat ’cap 
hi heat cap 

Gee i=») heat «cap 

w = moisture 


of soil. 


The values 1.0 and 0.5 
Capacities, expressed in units 


respectively: 


From Equation (IV-14), 
can be estimated provided that 
heat capacity of the soil ina 


capacities of dry soils vary wi 


he following relationships be used to 


soil-water system. 


Lt, OW 
a (IV-14) 
0.5 w 
F (IV-15) 


acity~of uUNTroOzen Sort,’ Btu per™ib per sr 
acity Om frozen soil, Btu per) 1b per’ F 
acity pot .dny soil... Btu, per..1b; per ih 


content of soil, per cent by dry weight 


in the above equations are the heat 


? 
of, (Btu per ibb, per F.,, of wateryand ice, 


the heat capacity of an unfrozen soil 
the dry density, moisture content and the 
dry state are known. Since the heat 


thin narrow limits an approximated value 


of C. would appear meaningful for practical purposes. lIlowever, in the 


derivation of Equation (IV-15) 


the assumption is made that the total 


moisture in a soil-water system is in a frozen state. The results of 


laboratory studies, such as reported by Lovell (1957) and Yong (1965), 
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have shown that practically all water in granular materials freezes at 
approximately 32 F, while, in soils containing substantial quantities of 
silt and clay size particles a large portion of the water may remain 


unfrozen at temperatures below 32 F. 


Unfrozen moisture content-temperature relationships of various 
soil types have been established by Yong (1965), Nerseova and Tsytovich 
(1963), Williams (1968) and Dillon and Andersland (1966). Although 
general agreement exists concerning the influence of silt and clay size 
particles on unfrozen moisture contents, there are differences of 
opinion as to the influence of initial moisture content on the percent- 


age of water frozen at a given temperature. 


Yong (1965) found that as the initial moisture content of a 
Clay having a liquid limit of 67 percent ;and aplastic limit, of 26 ‘per, 
cent, was increased the unfrozen moisture content at temperatures below 
32 F also increased. The unfrozen moisture content-temperature rela- 
eronships for a clay cand isilt soi1* sas obtained by Yonow (1965) ssare 
shown in FIGURE IV-2. Williams (1968) presented similar unfrozen mois- 
ture content-temperature relationships for various soils. However, 
Williams reported that different total moisture contents had no signifi- 


cant effect on unfrozen moisture contents. 


Nerseova and Tsytovich (1963) suggested that the unfrozen 
moisture content in a soil-water system is dependent on temperature, 
specific surface area of the soil particles, the mineral and chemical 
composition of the soil and the soluble compounds present in the water 


phase. Nersecova and Tsytovich suggested that for practical purposes 
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unfrozen moisture contents can be assumed to be primarily a function of 


temperature. 


Unfrozen moisture content-temperature relationships for 


nonsaline soils, as proposed by Nerseova and Tsytovich (1963), are given 


in TABLE IV-1. 


Dillon and Andersland (1966) developed an empirical relationship 


to predict unfrozen moisture contents of frozen clay soils. This rela- 
tionship is as follows: 
Slee 1 
Wo = To K. EL « k * 100 (IV-16) 
where: W, = unfrozen moisture content, per cent by 
dry weight of soil 
S = a measure of the average specific surface 
of the soil, particles; sq m per ¢. (M) 
Ty -=s) temperature of Prozens sol, Kk 
To = temperature of initial freezing of soil 
pore water, K 
A. = activity ratio, plasticity andex/percent. <= 21 
L = a constant, equal to.l and 2 for soils, having 
non-expanding and expanding lattice structures, 
respectively 
k= o2,8 X10 1548.08 water per. sa.M: 


Dillon and Andersland (1966) collected data on various soils, 


for which unfrozen moisture content information was available, and con- 


cluded that Equation (IV-16) yielded unfrozen moisture contents which 


were in close agreement with experimental values. 
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The results of the preceding investigations enable approximate 
temperature limits of a frozen and unfrozen condition of soil-water 
system to be defined. Between these temperature limits the heat 
capacity of such a system is made large by the latent heat of fusion 
released upon freezing. Williams (1968) determined the apparent 
specific heat capacity of various silt, clay and organic soils. 

Apparent specific heat capacity is distinguished from specific heat 
Capacity in that the former takes into account the latent heat of 
fusion, whereas, "true" specific heat capacity is not associated with a 
phase change. A typical apparent specific heat-temperature relationship 
of a soil, as presented by Williams (1968), is shown in FIGURE IV-3. 
With zeference to this figure,.at temperatures above the freezing 
temperature of a soil-water system the apparent specific heat is 
equivalent to the unfrozen heat capacity of the system. Upon freezing, 
the apparent specific heat is made large by latent heat, and as freezing 
progresses, unfrozen moisture content decreases and the apparent 
specific heat capacity approaches the heat capacity of a frozen system. 
The latent heat component of the apparent specific heat capacity can be 
approximated provided that an unfrozen moisture content-temperature 
relationship is known since the difference in unfrozen moisture content 


between any two temperatures gives the latent heat involved in such a 


temperature change, 


4.7. Summary 


In this chapter heat transfer concepts applied to pavement 


structures and thermal properties of pavement component materials have 
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been briefly reviewed. 


The influence of climatic variables on pavement surface temper- 
ature is difficult to quantify. A rational approach to this problem 
involves an energy balance procedure in which surface temperatures are 
related to various meteorological parameters. The accuracy of such a 
relationship is dependent on the method used to describe the influence 
of changing climatic conditions on surface temperatures, which will 
generally be governed by available meteorological data at a given 


locale. 


From the studies reported, thermal conductivity and heat 
Capacity values of asphaltic concrete paving mixtures have been found 
to vary within narrow limits and, therefore, estimated values would 
appear meaningful for practical purposes. Also, considering the 
moisture content of such paving mixtures to be negligible, these thermal 


properties may be assumed temperature independent. 


Thermal conductivities and heat capacities of granular bases and 
subgrade soils are primarily dependent on moisture content, density and 
temperature. In temperature prediction analyses of pavement structures 
subjected to low temperature climatic environments three different sets 
of these thermal properties must be recognized, depending on whether the 


material is in an unfrozen, freezing or frozen condition. 
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FIGURE IV—| HEAT TRANSFER BETWEEN GROUND SURFACE 
AND AIR ON SUNNY DAY. (after Aldrich, 1956 ) 
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INFLUENCE OF INITIAL MOISTURE CONTENT AND 
TEMPERATURE ON UNFROZEN WATER CONTENT 
OF CLAY AND SILT SOULS. (after Yong, 1965) 
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FIGURE TY-3 
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INFLUENCE OF TEMPERATURE ON APPARENT 
SPECIFIC HEAT CAPACITY. (after Williams, 1968 ) 
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CHAPTER V 


PREDICTED TEMPERATURES IN ASPHALTIC 


CONCRETE PAVEMENT STRUCTURES 


5.1. Introduction 

Numerical solution techniques, such as the method of finite 
differences, coupled with the use of digital computers provide a means 
of readily solving heat transfer problems which cannot ordinarily be 
evaluated analytically because of complex boundary conditions or geo- 
metrics. The prediction of temperature variations in a pavement 
structure as a function of time and position is one such problem. In 
this chapter heat transfer concepts and thermal properties presented 
in the previous chapter, together with a finite difference approxi- 
mation of the governing heat transfer equation, have been used in the 
development of an analytical model for predicting pavement temperatures. 
Comparisons of predicted and recorded temperatures in various asphaltic 
concrete pavement structures subjected to low temperature climatic 


environments are presented. 


5.2. The Finite Difference Approach 


The basic approach in a finite difference approximation of the 
governing differential equation of conductive heat transfer, Equation 


(IV-1) of the previous chapter, is the expansion of the temperature 
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functions in terms of a Taylor's series. Various methods, ranging from 
explicit to fully implicit methods, are available to express the 
differential equation in finite difference form. Independent of the 
choice of methods the problem reduces to a set of algebraic equations 


which can be rapidly solved by computer operations. 


The explicit method provides a noniterative process of deter- 
mining the temperature at a given position and time in relation to 
known preceding temperatures and assuming a one-dimensional analysis, 
only one algebraic equation is required to solve for the unknown 
temperature. However, the time increment at which the temperatures are 
evaluated is restricted by the mathematical stability criteria of the 


finite difference equation. 


Implicit procedures generally involve a matrix type of numerical 
solution at each time step since the unknown temperature is expressed in 
terms of present as well as preceding temperatures. The implicit method 
has no restrictions on the time increment at which the temperatures are 
evaluated. However, the resulting finite difference equations become 
difficult to formulate and program when two- or three-dimensional or 
complex one-dimensional heat transfer problems are to be solved. Such 
is the case in pavement systems where variations of thermal properties 
of the pavement materials with respect to temperature and location are 


required during program operations. 


Various investigators, Carroll et al. (1966), Straub et al. 
(1968), Corlew and Dickson (1968), Dempsey and Thompson (1970) and Ilo 


ect al. (1970) have successfully used the explicit method of the finite 
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difference approximation to predict temperatures in asphaltic concrete 
pavement structures. Many of the concepts and findings of these in- 
vestigations have been incorporated in the developed heat transfer 


model. 


yeoe Previous Investigations 


Carroll et al. (1966) used an explicit finite difference method 
to predict the depth of frost penetration in a pavement structure during 
a single winter. Daily average air temperatures and a constant con- 
yvecrion coefficient of 2.7 Btu per sq ft per hr’ per Fo were aneorporated 
within the air-pavement surface boundary condition. The position of the 
constant temperature lower boundary condition was found to influence the 
predicted thermal regime. The investigators concluded that for accurate 
deep temperature simulation studies, temperatures well below the region 


of interest must be known or carefully estimated. 


Straub et al. (1968) modified the computer program used by 
Carroll et al. (1966) in order that extreme daily temperatures in 
asphaltic concrete pavement surfaces could be predicted. The air- 
pavement surface boundary condition was formulated using an energy 
balance approach similar to that presented in Section 4.5 of CHAPTER IV. 
The lower boundary condition of a constant temperature was positioned 
at a depth of 24 inches below the pavement surface. At this depth no 
diurnal temperature variations were measured. The intensity of solar 
radiation was found to have a great influence on extreme daily asphaltic 
concrete temperatures and the authors suggested that when initial 


temperature conditions are unknown approximate temperature gradients 
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can best be estimated at a time prior to sunrise when the influence of 
solar radiation on pavement temperature is a minimum. Initial pavement 
temperatures of + 10 F from measured temperatures were found to have 
little influence on predicted maximum daily surface temperatures. 
Ilowever, with increased depth such input errors resulted in signifi- 


cantly large deviations betwecn mcasured and predicted temperatures. 


Corlew and Dickson (1968) computed the thermal cnergy trans- 
ferred in asphaltic concrete paving mixtures at the time of paving 
operations. The air-pavement surface boundary condition incorporated 
the radiation, conduction and convection modes of heat transfer. The 
results of the study indicated that, during paving operations, temper- 
atures of asphaltic concrete paving mixtures may decrease more rapidly 
at the pavement-subsurface interface than at the air-pavement surface 


boundary. 


Dempsey and Thompson (1970) developed a heat transfer model for 
predicting thermal regimes in pavement structures subjected to frost 
action. The solar radiation heat flux was accounted for by means of a 
semi-empirical relationship between solar position and the time of day 
at a particular geographic location. The longwave radiation heat flux 
and the convection coefficient were approximated by Equations (IV-5) 
and (IV-8) presented in the previous chapter. The latent heat released 
upon freezing of subsurface layers was estimated by varying the heat 
capacity of the component layers in accordance with a frozen water 
content-tempcrature relationship presented by Yong (1965). All soil 


water freezing was assumed to occur between 32 F and 30 IF. 
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Ho et al. (1970) compared predicted and measured temperatures 
in various insulated pavement structures subjected to subfreezing 
climatic conditions. Mean daily air temperatures were employed as an 
upper surface boundary condition and the latent heat of fusion was 
considered by specifying a nonlinear frozen moisture content- 
temperature relationship for each component layer. Predicted sub- 
surface temperatures were found to be influenced by the position of 
the constant temperature lower boundary condition. The authors 
compared recorded pavement temperatures with those predicted from 
one- and two-dimensional analyses and concluded that if temperature 
variations in the vicinity of the center line of a pavement structure 
are of interest, a good approximation of the thermal regime can be 


obtained using the one-dimensional analysis. 


The results of these investigations indicate that the 
meteorological variables required in the formulation of an air- 
pavement surface boundary condition are dependent on the problem to 
be analyzed. When predicting asphaltic concrete temperatures major 
metcorological variables should be specified. lIlowever, when pre- 
dicting subgrade temperatures, such a boundary condition may be 
somewhat generalized without greatly influencing the accuracy of 
the predicted temperatures. In both short-term and long-term 
temperature prediction studies the constant temperature lower boundary 
condition must be positioned at a depth well below the region in which 


predicted temperatures are of interest. 
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5.4. Formulation of the Finite Difference Equations 


When using the finite difference technique to obtain an approxi- 
mate solution of the thermal regime existing in a pavement system, the 
system is divided into a series of nodal elements as shown in FIGURE 
V-l1. The cross-sectional area of each nodal element is assumed equal 
to unity and the nodal depth (Ax) is so chosen as to satisfy the mathe- 
matical stability of the resulting finite difference equations. Also, 
the nodal depth is selected so that the pavement layer interfaces are 
located at the centroids of the nodal elements and at any time (t) the 
temperature at the centroid of a nodal element is assumed to represent 
the average temperature of the element. The prediction of temperatures 
of nodal elements located within a pavement component layer, at an air- 
pavement surface boundary and at the interface of adjacent component 
layers, requires the formulation of three general finite difference 


equations. 


5.4.1. Finite Difference Equation for Temperatures 
within a Component Layer 


The one-dimensional equation of transient heat transfer by con- 
duction, can readily be expressed in explicit finite difference form 
with the aid of the space-time grid shown in FIGURE V-2. The nodal 
points in FIGURE V-2 represent the centroids of the nodal elements at 
which the temperatures are to be predicted. The second order partial 


derivative of temperature with respect to depth can be approximated by: 
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and the change in temperature with respect to time can be approximated 


by: 
Oly A eHAt em T(x yt) 
ot At 


Multiplying the former expression by the diffusivity of the component 
layer, a, which is equal to the thermal conductivity divided by the 
product of the heat capacity and total unit weight of the material, the 


resulting finite difference equation can be expressed as follows: 


Re [T(x-Ax,t) = Cet) ae a IT (xtax,t) = OCs 


YCAx 
At 


Te BAG me end (XC) | (V-1) 


The left hand side of Equation (V-1) represents the thermal 
energy transferred to and from a nodal element, during a time increment 
(At) whose centroid is located at a space and time coordinate of (x) 
and (t), respectively. The expression on the right hand side of 
Equation (V-1) equals the thermal energy stored in the element during 
the time increment (At). Assuming that an initial temperature gradient 
throughout a pavement structure is specified, the unknown temperature 
T(x,t+At) can be calculated from a knowledge of the preceding tempera- 


tures. Rearranging the terms of Equation (V-1) and solving for the 


unknown temperature yields: 


T(x,t+At) = [R] T(x-Ax,t) + [1-2R] T(x,t) 
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where: R = a At/(Ax)2. 


Carnahan et al. (1969) indicated that a sufficient condition 
for the convergence of Equation (V-2) is 0 < R < 1/2. The selection of 
a time and depth increment which satisfy this condition is dependent on 


the thermal properties of the pavement component layers. 


5.4.2. Finite Difference Equation for Temperatures 


at an Air-Pavement Surface Boundary 

An approximation of the thermal energy balance at an air- 
pavement surface boundary is given by Equation (IV-9) presented in 
Section 4,9 Of the previous Chapter.. Usine a central ditterence 
technique, the thermal energy transferred by conduction can be equated 
to the thermal energy transferred by radiation and convection by the 


following expression. 


ane CT x=Axst) --"TOktaxte) | = 0 


+ 
2Ax = Q 


rad conv 


The term Rai equals the algebraic sum of the shortwave and longwave 
yadiation heat fluxes of the energy baiance equation. Solving for the 
fictitious temperature T(x-Ax,t) in the above equation and substituting 


into Equation (V-2) yields: 


T(x, jerat) = [ZR ePGe xO) 81-2R] alx yt) 


2 
i a [Qrad * Qconv] (V-3) 


The longwave radiation heat flux, incorporated in the Oe term 


of Equation (V-3), is dependent on the fourth power of temperature. 
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To the author's knowledge, the stability of such a finite difference 
equation, incorporating a non-linear temperature function, can only be 
evaluated by adopting a trial procedure in which various time and 


depth increments are specified. 


5.4.3. Finite Difference Equation for Temperatures 


at the Interface of Component Layers 

The development of a finite difference equation to approximate 
temperatures at the interface of two layers is based upon the assump- 
tion of a continuous heat flux existing at the interface. With 
reference to the ith layer, FIGURE V-2, Taylor's expansion yields 


the approximate relationship 


% oT a Lock 
T(x-4x,, ©) = E(x, co - Ax ) + oh (et) 


al 


and the time derivative is approximated by 


(22) i Niece - T(x,t) 


2 
Solving the former expression for <— and substituting the result to- 
gether with the latter expression into the one-dimensional Fourier heat 


transfer equation yields: 


Ax a = Te Ee -_ T(x,t + TAX t= xXx) 


ox i i 


Similarly, for the eae layer: 
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In order for the heat flux to be continuous at the interface the 


following condition must be satisfied. 


oT oT 
ci (=) an (=) 
5b all 


Substituting the above relationship into the two previous equations and 


si 


solving for the unknown temperature T(x,t+At) yields: 


T(x,t+At) = be T (x+Ax,t) 


5 Soin 
aes 
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T(x-Ax,t) (V-4) 


irande tor? stabidity,. the 
CKie te Nae i 
I 


Ine Equation: (V-4),. 1 = (Ax)< ys Oa 


lA 


weemperature coefficient of T(x;,t) 1s positive when 7 
Although the selection of the time and depth increment is 
generally governed by the stability requirements of Equations (V-2), 
(V-3) and (V-4), the use of small mesh sizes increase the accumulative 
effecr of round-off errors. It 1s difficult to determine exactly the 
order of magnitude of the cumulative departure of the solution due to 


round-off errors, however, such errors are probably far outweighed by 


those errors introduced when defining thermal properties and climatic 
variables. 
Equations (V-2), (V-3) and (V-4) form the numerical basis of 


the analytical model developed for predicting pavement temperatures. 
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The equations are well suited for programming on a digital computer and 
thereby allow rapid solutions to large amounts of input data. The re- 
quired input data consists of meteorological variables together with 
structural, physical and thermal properties of the pavement component 


layers. 


payee Input Variables for Temperature Prediction 


A detailed description, flow diagrams and program listing of a 
computer program developed to predict thermal regimes in layered 
systems is presented in APPENDIX A. The basic computer program consists 
of a main program and six subroutines whose operations interact to 
calculate temperatures over a specified time period and at depths and 
time increments’ of) (Ax) ‘and’ (At), ‘respectively. \ ‘Three other sub- 
routines, whose operations are independent of the above six, permit 
computer plotting of the generated temperatures at specified depths as 


a continuous function of time. 


eae ee ae Meteorological Variables 


Meteorological variables such as air temperature, solar radia- 
tion, cloud cover and wind velocity are often unknown quantities with 
respect to a particular geographic location and time. Therefore, these 
variables must be estimated from data collected at nearby meteorological 
stations. The extrapolation of meteorological data from one location to 
another may result in errors between predicted and field temperatures. 
The only time dependent meteorological variables incorporated in the 


computer program were solar radiation values and air temperatures. 
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Major meteorological stations throughout Canada record the 
hourly direct and diffuse solar radiation received by the earth's 
surface and these values are published in the Monthly Radiation 
Summary by the Department of Transport - Meteorological Branch - 
Canada. In order to approximate the incident shortwave radiation on a 
pavement surface at time increments of (At) a linear relationship 
between hourly input solar radiation values was assumed. At each incre- 
ment of time, the shortwave radiation influencing pavement temperatures 
was computed by multiplying the incident shortwave radiation by the 
absorptivity of the pavement surface. The absorptivity of asphaltic 


concrete was assumed equal to 0.90. 


The net longwave radiation and convective heat transfer entering 
the air-pavement surface energy balance relationship were approximated 
by Equations (IV-4) and (IV-7) presented in CHAPTER IV. The emissivity 
of asphaltic concrete was assumed equal to 0.95 and the surface coeffi- 
Gient of convection was assumed equal to 2.7 Btu per sq ft pervhr per. 
Such approximations do not consider the effects of vapour pressure and 
cloud cover on the net longwave radiation or the influence of wind 
velocity variations on the heat flux resulting from conductive heat 
transfer. Changes to the computer program can readily be implemented 


when attempting to account for the influence of these climatic variables 
on pavement surface temperatures. 
In order to compute the net longwave radiation heat flux and the 


heat flux resulting from convective heat transfer, at time steps of 


(At), a linear air temperature-time relationship was assumed. 
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5.5.2. Structural and Physical Properties 


The number of layers in the pavement structure in which tempera- 
tures are to be predicted together with the thickness, moisture content 
and dry density of each layer must be specified. The moisture content 
of each component layer was assumed constant throughout the time period 
in which the temperatures were calculated. Dempsey and Thompson (1970) 
suggested that the assumption of a constant moisture content is reason- 
able for predicting temperatures in unfrozen soils, while greater 
errors may be expected in frozen soil temperature predictions because 
the diffusivity of frozen soils is more sensitive to moisture content 


changes than the diffusivity of unfrozen soils. 


Geo. - Lhermal Properties 


From the literature review presented in’ the previous Chapter, 
the thermal conductivity and heat capacity of asphaltic concrete paving 
mixtures have been found to vary within narrow limits and for practical 
purposes can be considered temperature independent. In the temperature 
prediction studies presented in Section 5.6, thermal conductivity and 
heat capacity values of asphaltic concrete paving mixtures were assumed 


to be 0.84 Btu per fir per ft per Foands0322 Btu per 1b pera, 
respectively: 

The thermal properties of subsurface pavement components are 
dependent on whether the material is in an unfrozen, frozen or freezing 
condition. In order to distinguish between these three conditions, a 


freezing temperature range was specified for each subsurface layer. 
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Unfrozen moisture content-temperature relationships, such as summarized 
in Section 4.6.4 of CHAPTER IV, aid in estimating the temperature range 
in which various soil-water systems may be considered to be in a 


freezing condition. 


The method used to incorporate a moving zone of freezing in the 
model is most readily described with reference to FIGURE V-2. When 
solving for the unknown temperature T(x,t+At) "while freezing" thermal 
properties were used if the preceding temperatures, T(x-Ax,t) and 
T(x+4x,t), were within the specified temperature limits of freezing of 
the particular component layer. If this condition was not satisfied, 
the unknown temperature was calculated using frozen or unfrozen thermal 


properties depending on the temperature T(x,t). 


The unfrozen and frozen thermal conductivities of subsurface 
materials were estimated from figures presented by Aldrich (1956) and 
Sanger (1963), illustrating the influence of moisture content and dry 
density on the thermal conductivity of various soils as determined by 
Kersten (1949). The thermal conductivity of a component layer in a 
state of being frozen was assumed equal to the average frozen and un- 


frozen conductivity of the layer. 


The heat capacities of subsurface materials in the previously 
mentioned states were approximated by means of Equations (IV-14) and 
(IV-15) and the apparent specific heat concept proposed by Williams 
(1968), discussed in Section 4.6.4 of the previous chapter. The 
apparent specific heat at the specified initial freezing temperature of 


a component layer was assumed equal to the total latent heat released 
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upon freezing divided by the total unit weight of the component mate- 
rial. At the lower temperature limit of the freezing zone the apparent 
specific heat was assumed equal to the frozen heat capacity of the 
material. A similar approach to the problem of accounting for latent 


heat released upon freezing was employed by Dempsey and Thompson (1970). 


As shown in FIGURE IV-3 of the previous chapter, apparent 
specific heat decreases rapidly with decreasing temperature. Therefore, 
during freezing of a subsurface material a linear relationship was 
assumed between the logarithm of apparent specific heat and temperature. 
This relationship was established using the previously calculated 
specific heat values determined at upper and lower temperature limits 


defining the freezing zone. 


5.6. Predicted Versus Recorded Temperatures 


In order to assess the predictive capability of the analytical 
model a comparison was made between predicted and recorded temperatures 
in four asphaltic concrete pavement structures. Three of the structures 
were those incorporated in the Manitoba test section, and are referred to 
as Structures B, G and D, shown in’ FIGURE Jfl-1 of CHAPTER- III. “Tempera- 
tures in each of these three structures were predicted over a time period 
extending from mid-November, 1967 to the end of March, 1968, and at 


depths corresponding to the depths at which temperatures were measured. 


Bihourly air temperatures measured at the test site and hourly 
solar radiation values recorded at the Winnipeg International Airport, 
located approximately 35 miles west of the test site, constituted the 


meteorological input variables. The depth and magnitude of the constant 
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temperature lower boundary condition was approximated by means of 
FIGURE V-3, which shows the maximum and minimum recorded temperatures 
at various depths in each of the three structures during the 1967-68 
winter. From this figure, a constant temperature of 50 F was imposed 
at a depth of 18 feet below the pavement surface of each structure. 
This depth is 6 feet below the maximum depth at which temperatures were 
recorded and approximately 9 feet below maximum recorded depth of frost 
penetration. Initial temperature gradients were determined from re- 
corded temperatures at the start of the prediction period when all 
component layers were in an unfrozen condition. Input parameters for 
the temperature prediction of Structures B, C and D are given in TABLES 


V-1, V-2 and V-3, respectively. 


The fourth pavement structure in which a comparison was made 
between predicted and recorded temperatures is located approximately 10 
miles east of Edmonton, Alberta. Input parameters for the temperature 
prediction of the structure are summarized in TABLE V-4. Within this 
structure, temperatures at the 1/2, 2-1/2, 5 and 7 inch depths in the 
asphaltic concrete pavement surface and at 1 foot intervals, to a depth 
of 9-1/2 feet, in the subgrade were recorded on approximately a weekly 
basis during the 1968-69 winter. The time period of temperature pre- 
diction extended from November, 1968 to April 1969, during which time 
an identical lower boundary condition to that imposed at the Manitoba 
site was assumed. The initial temperature gradient existing in the 
pavement structure on November lst was estimated and the hourly air 


temperatures and solar radiation values used were obtained from records 


wor 
16 easomt ora xbin ‘zou aah tibiies crabavod “taMp I OF 
sotutexanne. ‘lpeago mune ‘has cotixem odd ewote dotsdiw 4) 
Boel GGL any al id eamiriouede seus old to toss ML, 
bezoqut 2ew 7 2 te ee TEFAAGD 6 org i? eit aor 
serutovrae 189 to Scam thoueveg ot voled Jae? £20 “ 
oro ects teteqie? Hopi tx Tash mmixwm oft woled 396%, o al a 
faqx? to isqeb bolrosox minkxen woled jee? @ + eee bas 
~9T mex boninrstob stow daveatigny ieihiabiaain’ isidink ae 
. Lis aaiiw poled, eae fhetg dit Je Sxete old tu sorusageend : 
tot sxe Tanetae ah tap E4ibneo rosevine ‘ne ot otow a) 
ees eV ne . 


, oben, easy Soet sais B dakiw nt savtoutite ounbeke dtxot ee 
Ol fotemixongds Letase! et es ittpraqnss bebrooe: bre ‘Desout ; fide 
Srerei yc! bg a to) sratems sy tunel etusdiA . notaombit +0 anal 
Shits (senna ~ AV aap Kt Sesirammee O18 weuTouUtI2 onl? %0 in oth 
sid thi zig dit woae 32 D\E~ $ eS\i ond a eo1uIeTegAOT a nk 

ot 


dtqeb fe ‘od, ah peat soe i. 7s nts SostTHye Inomevaq ovoranog 2 
yiAaoW Bb aa pane Sioa } Lobrosse orow shia adie ait ni ook 
Poet oxuas 18 wie bot) omit sel abtniw e9-Bae! ons ae 


“8 “A o 4 : x t) ae 
ay fo rit aes arr =O ot J tre | on Babs taddev onl mos? babnasxe Ls 


ae fa 4 

5 . : 6) nobhtibaad yesbaued towol cians i, 
ener (aitiat ed? sbomueds aw stte 
iro cow Jal yedHoveK no sTutsuNs Gnemeved’ 


-seuley ndieeéhee tsfoe bre somite teqnes 
hal) eee Gere 


108 
taken at the Edmonton Industrial Airport and Stony Plain Meteorological 
Stations. These stations are located approximately 10 and 25 miles 


west of the pavement structure, respectively. 


Time and depth increments of 0.125 hours and 2 inches, respec- 
tively, were employed in each of the temperature prediction analyses. 
These values were found to be compatible with the stability require- 


ments of Equations (V-2), (V-3) and (V-4). 


Ses et Asphaltic Concrete Temperatures 


Comparisons between predicted and recorded temperatures in the 
asphaltic concrete surfaces of Structures B, D and C, during low 
temperature climatic conditions, are presented in FIGURES V-4 (A) and 
(B) and FIGURE V-5 (A), respectively. Since all meteorological input 
variables and asphaltic concrete thermal properties used for the 
temperature predictions of the Manitoba pavement surfaces were iden- 
tical, predicted temperatures at equivalent depths and times in each of 
the three structures were similar. However, recorded asphaltic concrete 
temperature distributions in each structure at any particular time were 
found to differ. As mentioned in Section 3.3 of CHAPTER III, Young et 
al. (1969) suggested that such differences may partially be accounted 
for by micro-climatic variables existing at the test site. In the 
prediction analyses, no attempt was made to account for the influence 
of micro-climatic conditions, such as differences in wind velocity at 


the air-pavement surface boundaries, on pavement temperatures. 


During the period of seasonal minimum ambient temperatures, 


predicted temperatures of Structure B were approximately 4 F higher than 
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recorded temperatures, while, during the same period predicted tempera- 
tures of Structure C were approximately 4 F lower than recorded 
temperatures. The difference between recorded and predicted tempera- 
tures of these two structures partially reflects the difference in 


measured temperatures in the structures at equivalent depths and times. 


Temperature predictions for Structure C were continued into 
the month of June, 1968. A typical comparison between predicted and 
recorded asphaltic concrete temperatures, during a time period in which 
ambient temperatures were approaching seasonal maximum values, is shown 
in FIGURE V-5 (B). Excellent agreement between predicted and measured 


temperatures during such a time period is evident. 


In order to obtain a measure of the accuracy of the predicted 
asphaltic concrete temperatures of Structures B, C and Dover the 
entire 1967-68 winter, linear regression analyses were employed to 
relate maximum daily recorded and predicted asphaltic concrete tempera- 
tures and minimum daily recorded and predicted temperatures in each of 
the three structures. The results of these analyses are presented in 


TABLES V-5 and V-6, respectively. 


A similar approach was used to compare the predicted and 
measured temperatures in the asphaltic concrete pavement surface 
located near Edmonton. However, in the latter analyses, predicted and 
recorded temperatures were compared at times corresponding to times 
of temperature measurement. The recorded-predicted temperature rela- 


tionships of the Edmonton structure are presented in TABLE V-7. 
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As observed from TABLES V-5, V-6 and V-7, maximum standard 
errors of estimate were associated with temperatures at and immediately 
below the air-pavement surface boundary. This result was expected 
Since many meteorological variables such as snow cover, wind velocity 
variations and cloud cover, which are known to influence the pavement 
temperatures, were not included in the formulation of the air-pavement 
surface boundary condition. Also, such errors may partially be due to 
the extrapolation of meteorological data from one locale to another 
and errors in assumed thermal properties of the pavements. Despite 
the many assumptions involved in the development of the analytical 
model, the comparisons between recorded and predicted temperatures 
suggest that the model enables temperatures in asphaltic concrete 
pavement surfaces to be predicted within limits of practical 


Slonitieance . 


5.6.2, Subgrade Temperatures 


Isotherms depicting the predicted and recorded thermal regimes 
which existed in Structures B, C and D and the Edmonton Structure are 
shown in FIGURES V-6, V-7, V-8 and V-9, respectively. The depths of 
the isotherms were calculated assuming a linear relationship to exist 
between temperatures at various thermocouple positions. Due to the 
relatively small temperature gradients in the subgrades, small differ- 
ences between predicted and recorded temperatures at a given depth and 
time may result in large variations between the positions of the re- 
corded and predicted isotherms. Therefore, it may be considered 


fortuitous that the distance between corresponding measured and 
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predicted isotherms at a given time was generally less than the 


distance between thermocouple positions in each of the subgrades. 


5.7. Summary 

In this chapter an analytical model for predicting thermal 
regimes in pavement structures has been developed and comparisons be- 
tween predicted and recorded temperatures of four pavement structures, 


located in Western Canada, have been presented. 


The close agreement between predicted and recorded tempera- 
tures in each of the four structures suggests that the model can be 
used to obtain reliable predictions of thermal regimes existing in 
various pavement structures in different geographical locations. Such 
temperature data aids in realistically assessing the influence of 
climatic environment on the response of pavement systems to various 
loading conditions. In relation to the present study of low temperature 
transverse cracking of asphaltic concretes known or predicted asphaltic 
concrete temperatures, constitute primary inputs for thermal stress 


computations and tensile strength determinations. 
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TABLE V-5 


RELATIONSHIPS BETWEEN MAXIMUM DAILY RECORDED AND 
PREDICTED ASPHALTIC CONCRETE TEMPERATURES OF 
STRUCTURES B, C AND D 


, : Coefficient Standard 
Structure PEE Vis BSS STON of Error of 
d Correlation 


Estimate (F) 


C 
D 
% Tops Top = Maximum daily recorded and predicted surface 
temperatures (F). 
Tops Top = Maximum daily recorded and predicted 2 inch 


temperatures (F). 
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TABLE V-6 


RELATIONSHIPS BETWEEN MINIMUM DAILY RECORDED AND 
PREDICTED ASPHALTIC CONCRETE TEMPERATURES OF 
STRUCTURES B, C AND D 


. : Coefficient Standard 
Linear Regression oe mel mc 
Correlation Estimate (F) 


Equation* 


Minimum daily recorded and predicted surface 
temperatures (F). 


Minimum daily recorded and predicted 2 inch 
temperatures (F). 
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TABLE V-7 


RELATIONSHIPS BETWEEN RECORDED AND 
PREDICTED ASPHALTIC CONCRETE TEMPERATURES OF 
THE EDMONTON STRUCTURE 


Standard 
Error of 
Estimate (F) 


Coefficient 
of 
Correlation 


Linear Regression 
Equation* 


Tp = Laake & 1.03 Tp ase )s) 
Tp = 1.36) + 1.01 Tp (hone 


Tp =f LORE O97 Tp 0.99 1.42 


* Predicted temperatures calculated by assuming a linear relationship 
between calculated temperatures at 2 inch depth increments. 
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FIGURE Y-2 DEPTH-TIME GRID OF A LAYERED SYSTEM. 
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ASPHALTIC CONCRETE TEMPERATURES OF STRUCTURE D. 
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ASPHALTIC CONCRETE TEMPERATURES OF STRUCTURE C. 


24 2 24 2 


\ Predicted ‘a—a Recorded «——e 


95 \ temperature *~4\ temperatures 
xs ! i BS 
* meen 
i U \ 
vo Surface ; 

85 +—-_—_ /} ~——- ! 
ris 
° 
- 
= 75 
oS 
o 
a 
E 
® 
te 

65h- 

rn 
7 eT 
SS ~ 
ao = - 
June 16, 1968 | June I7, 1968 | June 18,1968 | 


FIGURE V-5(B) COMPARISON EBETWEEN PREDICTED AND RECORDED 
TEMPERATURES OF STRUCTURE C DURING SEASONAL 
MAXIMUM AMBIENT TEMPERATURES. 


> cj be gr wer 
ae mate ea me “Ht PrEPIL ate | 


¥. ‘ . ; J , v4 m 
I | oe Pol in = AL yy 7 
a ‘e - + 
oe RN 
Ca Re ete ie a 


— . wf. J 

t ~s,/ y a . 7 S3D Ti? * 

: ¢ ’ ' = ‘ 

poe 5 4 

; : ‘.s 7 - 7 AF 

| ie ye 

. . a o——-» sufgisqma! babosoe 

i | we ; ane au ID Sgn! bs isloge 
t ; 

. 


Bae: yigunal. |} Gae! © yywine . BART Sveroetiol 


CAGHODSA GHA GSTDOARSS wadwTss NOG AMOS WEY 3 
3 SRUTIURTS IO SSHUTARSIMST STAHINOD KTIANTA 


Ps As! BS Si PS Ss) ; ps 
ie err ae : — atsed ra —— i} ; ~~ 
A : tre bby O94 ii +e befits) j 7 
emel | (! amuitpragmat je *\~———. sadtoveq OTS ere rae se 
F A 5 a 
| rer 


' 
1 adaeagl 


1A 


. x 
ant Wes y, ro | a 
| a f iy tx? 
a es | 


ie 4 ; , - 
 .  GaaRONFe ay eer VINE SD 0 ou 
 « AdAMORP eS Ost SUuHTe IO 2 Vm 


LARS HIT — 


124 


-g@ gunionyls dO 3qvYoaNsS SHL NI 3WI93Y 
TIVWWHYSHL G30N00aY GNW G3Ll9IdaYdd 3HL N3SM1458 NOSIYVdWOO po av Aer f= REO | 


“SDAA JO OWL 


G2 6c. Sc Oc’ SM OOF GS of Sse Oc S$ QO G o€ Ge O2 SI Ol G O€ Sd 
02 | 


g96| ‘Asonsqay 


0O0l 
SS ee Q3L90I03ud 
g3a0qu¥0034 


mil 


PRC TAI ADAP LO Sa RNR TR PH Se IRF SN EG oS PRE 


(sayou!) yideq 


i 


0 2 ra oe a 
SCT “ ¢ - ~) 


| Osa0n0aaa': 
OST SGass:) 


——— 


- i : er . ' = : = P aa 7 Sa a 
Ke ; 
5 


bee 
| Ss | | 
- A x 
f= j - ; — 2 os * = 
G3el, yiourdeay at a Pa S2eh,yiounot| Feel isdmssed 
a Ya 
eI if 
ae Se LES MND Ont ORONO Ss QeeeeeS (OAc cmae® ROcnemes Neeonees nn cen Os Oe See Oe earner 
Se 8S cs OS 2 Gg “8 Of cS. OS eC. oo c Ce “cS OS- cl Oo & . 06 


ay; ‘3 13595 ASAHT -“WSsWTSSE. UOetHASHOD g-Y SAMOA 
a G SHUTIURTS FO SCAMeEeUe SHAT Vi swMtean 


bZ5 


"2 SYNLONYLS 40 3QVYUOENS SHL NI SWIDSY 


IWAWYSHL G3dq4003SH ONV G3LdIGSYd SHL N33ML3E NOSINVdNOD 2Z-A 3yNDdIS 


JDaK jO owly 


G 62 Se O2 Gi Ol S Of G2 02° 7S1 = 0) S 


896] ‘Asoniqs4 §96] ‘Asonuor 


Og Sz 


O2 SI Ol 


aa 


eee OO OIC aed 
03040034 


CRA stair LAT BRON = z 7 Z ees 
“A “Gal? bare 


ia Mee a Sa we acs SP MOTE POM BA Cy ETB as 
ry ree = 02-20 .°-2 O--C= oS, Prxc3 EY? eO0s 
oS a ST RS 


ro ew a 


CEO OR DL 


; v 


Sete Nes a Eee oie eS eS 
De ORL RAP PB PI BET} 


902 aS 


Cig Ci eee eas Sraeere egos 6a. Ae Sy pres 
Eo wie eRe Tee ects fake ee ee 


VERE 


(seyoul) yydag 


: 
= ee | 
—" 
Ae 
a 


2 « = 


" 


be ia naa 


¥ 
@ 
¢ 
Z 


enews’ QSCAGDAA 
= — OST aAF | . 
ee dees - SS ee 


a Se — — 

t r a cs 

Lo =| 980, yiowdes 836i , vinunol TSE! wwdmsos0 1 

| % ae 

Se oe SE SR SY SES Goes Me Ma SS i Se es hoe eee! Ronee eee CHES 

@ es zs 0S 4 Ol -@ c€ && of Of) Ob 6 0F E8505 cf Oe Ss 

wey tovemiT 
JAMABHT GRdR0D7A GVA-C3TOUGIAS. SH WASWTae wOeIRAGMOD \-V SRMBte: 

= ) BQUTIUATe FO: SGAROsUa—SHT Mi SMIS3A = 
“i 7 


126 


‘q a3yunLonuis 30 3qVYOENS SHL NI AWIDSY 
TWWYSHL G30N003Y GNV GaLDIGSYd SHI NaaM138 NOSINYVdWOD 8-A 3YNSIS 


"409A ayy JO QWIL 
¢ 62 S$2@ O02 SI OI G o¢€ Se Oe sg OI S Og G2 O2 SGI Ol= -S. Of Se 


om 896| 'Asonsqe4 g96| ‘Asonuor | 196! 's4squedeg 


— ———— peyoipeid 
a 010904 


(seyoul) yydeq 


4 rans abies : i i 
WF isegcae Scetees Ost 


ee 


ages beizibs:4 


: as | aw a 2 a 
| 
[ 
' 
i 


Bee) yiowdst= soci, YiIounDL- : Tel sedmsosc “| 
es One Ln ee See Os Ce inne nee a ORS SS ee ee Ce 


: JAMABHT G209003R GMA GaTIOIRe SHT usBWIZ8NOsIRAIMOD B-V Balola 
3 qd 3auTsuate 30 30AnDBVe SHT Ul SMIOSA 


127 


SWISS 


S 


‘3UNLONYULS NOLNOWQ3 3HL JO 3q0VYOENS SHL NI 
3Y TIVWYSHL G304003Y ONV G3L0IG3ed 3H1 N3aM1358 NOSIYVdWOD 6 


"108K ayy JO OWIL 


G2 O02 O€ Se Oc SI Ol S Of G2 O2 si Ol S 


696] ‘ Asonuor 


————— pesoiperd 
aa =pTPsOIAY 


We are, 256— 


Boies. one phe 


A 34yndl4 


(saysu! ) ujdag 


ee 


ese} qntuids3 Csel -yreunol -- 
a Oo. @ as os. @- OQ 2 8 Sf Bs os “ah Ol 
.18y St Jo emit 


SMiJah UAMASHT OUSOROISR GUA C3STO!IOSRI SHT MSSWTad YOsinAIMOD e-V 280 
. AWHUTIUATe YOTKOMGS SHT 30 SGAROEUS SHT WI 


Xa 


CHAPTER VI 
MATERIALS CHARACTERIZATION 


6.1. Introduction 

During the past two decades various methods of determining the 
time and temperature dependence of the deformation characteristics of 
asphalts and asphaltic concrete paving mixtures have been developed. 
As mentioned in Section 2.4.3 of CHAPTER II, these methods involve 
direct laboratory testing techniques as well as indirect procedures 
which utilize nomographs from which the time and canperaeire dependent 
stiffnesses can be estimated from a knowledge of the physical properties 
of the material. Previous studies related to the craneverte cracking 
problem have made use of both direct and indirect procedures to 
characterize the behavior of asphalt materials subjected to low tempera- 
ture climatic environments. In the present study numerical methods used 
to predict thermal stresses in asphaltic concrete pavements require 
definition of the stiffness characteristics of the asphaltic concretes 
and in order to assess the low temperature fracture susceptibility of an 
asphaltic concrete paving mixture subjected to a particular temperature 
history, assuming a stress fracture criterion, the fracture strength- 


temperature relationship of the asphaltic concrete must be defined. 
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In this chapter the concepts which were used to characterize the 
stiffness and low temperature tensile strengths of asphaltic concrete 
paving mixtures incorporated in the Manitoba and Alberta test projects 
are reviewed and values of the necessary parameters to characterize the 
materials are presented. The stiffness values, together with laboratory 
determined tensile properties of asphaltic concrete cores obtained from 
the test sites, are briefly discussed in relation to observed differ- 
ences in low temperature behavior of the various pavement test sections 


included in the projects. 


6.2. Pavement Test Sections Included in Study 

A brief description of the Alberta test project has. been pre- 
Sented in Section, 3.2.1 of CHAPTER III. In the following chapter 
comparisons are made between observed and predicted times and tempera- 
tures of fracture of each of the three test sections having a pavement 
structure as shown in FIGURE III-1 (A). For these comparisons tensile 
strength-temperature relationships of asphaltic concrete cores, 
obtained from the test project at the time of construction and at a 
service age of 34 months, were determined using the tensile splitting 


test procedure reported by Christianson (1970). 


The design variables included in the Manitoba test project have 
been reviewed in Section 3.2.2 of CHAPTER III. In the present study 
observed and predicted low temperature fracture of ten of the twenty- 
nine pavement sections, during the 1967-68 winter, are compared. With 
the exception of asphalt and pavement structures, the design variables 


of cach of the ten sections are effectively the same. The asphaltic 
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concrete paving mixtures of the ten sections were those having optimum 
asphalt contents and a non-modified aggregate gradation. The sections 
under consideration include; three test sections of Structures B, C and 
D, shown in FIGURE III-1, which incorporated the low viscosity (LV) 
150-200 penetration grade asphalt, three test sections of Structures B, 
C and D which incorporated the high viscosity (HV) penetration grade 
asphalt, two test sections of Structures B and C which included the LV 
300-400 penetration grade asphalt and two test sections of Structures B 
and C which included the slow curing SC-5 liquid asphalt. Fracture 
strength-temperature relationships of asphaltic concrete cores obtained 
from the test project during the summer of 1968, and composed of the 
four asphalts, have been utilized in the fracture predictive analyses 


of the pavement test sections. 


Due to a limited number of asphaltic concrete cores obtained 
from the test projects, laboratory testing conditions were restricted 
and, therefore, indirect nomographic procedures were used to estimate 
the stiffness of the asphaltic concrete paving mixtures over time and 


temperature ranges necessary for the thermal stress predictions. 


6.3. Determination of Stiffness Modulus 

An extensive review and assessment of various procedures used to 
characterize the viscoelastic behavior of asphaltic concrete mixtures 
has been presented by Monismith et al. (1966). Included in this review 
were the concepts of stiffness proposed by Van der Poel (1954) and the 


concept of time-temperature superposition, both of which were used in 


the present investigation. 
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Van der Poel (1954) introduced the concept of a "stiffness 
modulus', which is time and temperature dependent, to characterize the 
response of asphalt materials to various loading conditions. This 


stiffness modulus, commonly shortened to stiffness, is defined as: 


S(t ‘ri eecan (VI-1) 
(Se 
where 
S(t,T) = stiffness at a particular time and temperature 
O, € = axial tensile stress and strain, respectively. 


On the basis of correlations obtained between results of creep 
and dynamic tests on various asphalts and results of routine tests on 
the asphalts, which included penetrations and ring and ball softening 
points, Van der Poel developed a nomograph from which the stiffness of 
an asphalt can be estimated as a function of time and temperature. This 
nomograph is shown in FIGURE VI-1. The required characterizing asphalt 
parameters for computing stiffness from the nomograph include asphalt 
penetration and penetration index. Assuming that the penetration at the 
temperature of the ring and ball softening point equals 800 dmn, as 
suggested by Van der Poel (1954), the penetration index, as defined by 
Pfeiffer and Van Doormaal (1936), can be computed by means of the 


following cquation. 


z - log P 
_ : - - s¢Q ion i = en (VI-2) 
R&B P 


where 


PI = penetration index 
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Pen = penetration of asphalt, dmm 
TREB = ring and ball softening point, C 
Tp = temperature at which penetration is determined, C. 


From comparisons of experimental values of asphalt stiffness 
with values read from the nomograph, Van der Poel found that nomograph 
stiffness values were generally within a factor of two of those deter- 


mined from direct tests. 


Heukelom (1966) presented a version of Van der Poel's nomograph 
with a slight correction for asphalts having very low penetration index 
values. Further modifications to the original method of determining 
stiffness have been suggested by Heukelom (1969) and McLeod (1969). 
These modifications, together with the use of the stiffness concept as 
applied to the characterization of asphalt paving materials, have been 
reviewed by Anderson and Haas (1970). Included in this review was the 
particular method used in the present study to characterize the time 
and temperature dependence of asphaltic concrete paving mixtures in- 


cluded in the test projects. 


Van der Poel (1954) recommended that FIGURE VI-1 not be used to 
estimate stiffness values of asphalts having relatively high wax con- 
tents since the ring and ball softening point of such an asphalt is 
influenced by the presence of paraffin wax. Recently, Heukelom (1969) 
suggested that the ring and ball softening point of waxy and/or blown 
asphalts be modified and developed a Bitumen Test Data Chart from 
which a modified softening point can be derived and used in association 


with Van der Poel's nomograph for estimating the stiffness of such an 


Get 
cd Terttigs Ro HoLresionsy) = ‘mite 
a TH EOG ginetion Llad base. grit , = ee 


barnierteteb 2k Sap aences gorde 35 onus andlor + gl aa 
‘ ; 


Peat t Tye eer to phetkey fescomrrodme bo actor ireqMos aac 
dqsinonan. teas Dayet sot coer ‘linkinaiee sit woud bhsst >On. 
-tsi9b szont to owt to Yoton 6 oldtiow viiatensg stew esuley 


‘ : i 
.efee? soeribh mort Be 


am 
1 


: 
; ‘ i 7 
iqgsryomon 2 feo Bins) To noxetsy if besser (0001) moledvet, — 


cabn} nroltertedneq Wol yrev gaivenl etenges tot #OLIsetTO9 ae . 7 


Si 
eg 


(edi) bosdol bas (eel) ‘moltowtialt yd beresggue aged evade ott. 


gisrraressb ta bonssa Terres “Oo 6Hhs oF ans IRIE Tne control 


28, FQSOOo avent4 ite oid 20 Saux sit (ein sentggo2 ene taas f ee 
ised svet ,elsisétam yaiveq Dleiqes to mnisexrivesosteds ods\os beth 
sit egwewokves 2tat nt bebulonl | (OVER) 2pah bas nosteten Ya errs 
: ’ 
sniy ott Ssbxedogseds oF Koya? yrecese ld at boew bodsem team 
-ri zetutxim yotver Stetonop citi signe Re sqnebneqeb - ; m 1b 
eiyatorg Jeeteds ait 


[aw 


os boas sd forms WM ‘aug oe babpsimain 9a pels) tsot rob - 
| hed as <== 7 
“100 *EW agen (isviaelay anivail ao tedaee Yo eoulsv eeanthive 


my : 
ef Jladqes Ae gira La aa! grinoszte hie ies ‘pnts ad iia 
(Q021} molsduel cet sinen9s KEW nitereg, a soTsZe1] ong a bi 


mwofd to\bis sexi to sels gilinottee ist: bre anit 9 a a gue 
owt teeoeed ica? loans iA stbeggt ova bre betta om od 23 isdqes 
noitsipgozes ni beeo bis haviveb, ad, RSS ‘oman belts iboa s dod 
we we 
7 8 ilove to 2eontxig ni a (Oy Neston: ae eb | 
”, bai 9 pc ung ee 


133 
asphalt. This modification, as applied to some Canadian asphalts, has 
been described in detail by Kopvillem and Heukelom (1969) and was used 


for estimating stiffness values of the paving mixtures. 


Properties of asphalts included in the Manitoba test project, 
sampled prior to plant mix operations and after being recovered from 
uncompacted field mixes, have been summarized by Young et al. (1969). 
Stiffness computations for the penetration grade asphalts at the 
Manitoba test project were based on index values, shown in TABLE VI-1, 
which were derived using as-supplied asphalt samples. The choice of 
using as-supplied asphalt properties rather than the properties of the 
asphalts extracted from uncompacted field mixes was made in view of 
Tesultsixeported by Shields et al. (1969) and Bright et al. (1969). 
These investigators suggested that asphalts recovered from uncompacted 
field mixes may exhibit premature hardening characteristics unless 
special sampling procedures are used and, therefore, the physical 
properties of such asphalts may not be indicative of the actual 


physical properties of asphalts within newly constructed pavements. 


Comparisons between estimated stiffness values of the asphalt 
cements included in the Manitoba test project, and determined using as- 
supplied and recovered asphalt properties, have been reported by Young 
et al. (1969). The low temperature stiffness values of a given asphalt 
cement, determined using the two sets of asphalt properties, were 
found to be approximately equivalent. This finding reflects the off- 


setting influence of the lower ring and ball softening points and 
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penetration index values of the as-supplied asphalts, as compared to 
the corresponding properties of the recovered asphalts, on estimated 


stiffness values of various asphalts. 


The properties of the asphalts incorporated in the Alberta 
test project and determined at various service ages have been reported 
by Anderson and Shields (1971). For stiffness computations associated 
with the Alberta test project both as-supplied properties and properties 
of asphalts extracted from pavement cores at a service age of 34 months 
were used with modifications to the index values associated with Van 
der Poel's nomograph in a manner suggested by Heukelom (1969). 
Properties of the asphalts used for stiffness computations are given in 


TABLE Vie2. 


Using Van der Poel's nomograph, together with penetration index 
and ring and ball softening point values given in TABLE VI-1 and VI-2, 
the stiffness of each asphalt was determined at loading times ranging 
from 10 to 100,000 seconds and over a temperature range of 40 to -40 F. 
The corresponding asphaltic concrete stiffness values were computed by 


means of the following equation, proposed by Heukelom and Klomp (1964). 


S C H 
mix 3 ake CuO. Vv (VI-3) 
S n L=Cy 
asphalt 
where 
os = mixture stiffness, kg per sq cm 
mix 


Sisphalt = stiffness of asphalt, kg per sq cm 
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C = volume concentration of aggregate 


defined as 


volume of aggregate 
volume of (aggregate + asphalt) 


4 x 10° 


=) 
i 


0.83 log 
asphalt 


Heukelom and Klomp indicated that [quation (VI-3) is valid for 
asphaltic concrete mixtures with C\, values ranging from 0.7 to 0.9 and 


air void contents in the order of three per cent. 


The volume concentration of aggregate values for the Manitoba 
asphaltic concrete paving mixtures were computed using specific gravity 
values and average asphalt contents of uncompacted field mixes, re- 
ported by Young et al. (1969), together with determined densities of 
Dophattic concrete cores obtained from the test project. The C, values 
of the Alberta pavements were computed using similar properties reported 
by Shields et al. (1969). For those paving mixtures having air void 
contents greater than three per cent the C, values were modified in 
accordance with the procedure proposed by Van Draat and Sommer and 
described by McLeod (1969). This modification involves correcting the 


determined C, value of such a paving mixture by means of the following 


equation 
C 
Ce (VI-4) 
MA 1+AV 
where 
C," = corrected volume concentration of aggregate of 


those mixtures having air void contents greater 


than three per cent 
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AV = difference, expressed as a decimal fraction, 
between the volume of air voids in the paving 
mixture and an air voids volume of three per 
cent for the same paving mixture. 

Physical properties of the asphaltic concrete paving mixtures 


associated with Manitoba and Alberta test projects and included in the 


present study are summarized in TABLES VI-3 and VI-4, respectively. 


Having estimated the stiffness modulus of each asphaltic 
concrete paving mixture over the previously mentioned time and tempera- 
ture ranges, the time-temperature equivalency hypothesis was used to 
establish a master stiffness modulus versus reduced time relationship 
for each of the mixtures. The use of this hypothesis, as applied to 
Materials characterization of asphaltic concrete, has been summarized in 
the review presented by Monismith et al. (1966). Briefly, according to 
the hypothesis, all basic time dependent response functions, such as 
stiffness modulus, are affected by a temperature change only within a 
corresponding uniform shift of the logarithmic time scale. To illu- 
strate this concept, stiffness curves as shown in FIGURE VI-2 may be 
shifted along the time axis to a selected reference temperature, Tp), 
and a reduced time, £, is defined such that stiffness curves for 
various temperatures superimposed when plotted..as, a function~of 1G. 
Materials exhibiting such behavior are termed "thermorheologically 
sample. 


The reduced time & is defined by the equation 


—E = t exp[f(T)] (VI-5) 
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where 
E = reduced time corresponding to a 
real time, t, at temperature T, 
f(T) = a temperature function giving the shift 


of the response curve. 


The temperature function f(T) is generally termed the shift 
factor, ans and equals the ratio of the time to observe a phenomenon at 
a temperature T to the time to observe the same phenomenon at the 


reference temperature To. 


6.3.1. Stiffness Modulus-Reduced Time Relationships for 
Manitoba Asphaltic Concretes 


The determined stiffness modulus versus reduced time relation- 
ships and the shift factor versus temperature relationships for the 
Manitoba asphaltic concrete mixtures are shown in FIGURES VI-3 and VI-4, 
respectively. Such relationships permit the stiffness of the asphaltic 
concretes to be readily evaluated at any time and temperature and con- 
stitute the primary inputs of the stress predictive methods used in the 


following chapter. 


As observed from FIGURE VI-3, at 0 F and a given reduced time, 
which is equivalent to real time at this reference temperature, the 
four Manitoba asphaltic concrete paving mixtures were found to exhibit 
large differences in estimated stiffness values. This would suggest 
correspondingly large differences between the magnitude of thermally 
induced stresses within the various pavement surfaces when subjected to 


similar low temperature climatic environments. 
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Although a realistic assessment of the low temperature fracture 
susceptibility of an asphaltic concrete paving mixture, based on a 
stress fracture criterion, requires a knowledge of both thermally 
induced stresses and tensile strengths of the mixture, a correlation 
was found between relative stiffness values of the Manitoba asphaltic 
concretes, shown in FIGURE VI-3, and transverse crack frequencies of 
the pavement sections in which the paving mixtures are incorporated. 
These frequencies, as reported by Young et al. (1969), are summarized 
in TABLE VI-5. Differences in crack frequencies between pavement 
sections which incorporated the same asphalt cement may partially be 
attributed to differences in asphaltic concrete thickness and subgrade 
soil type, as discussed in Sections 2.4.6 and 2.4.7 of CHAPTER II. 
However, those paving mixtures composed of the low viscosity asphalts 
exhibited the highest estimated stiffness values, at low temperatures 
and/or short loading times, and the highest crack frequencies, while, 
those mixtures composed of the high viscosity penetration grade asphalt 
and the SC-5 liquid asphalt exhibited the lowest stiffness values and 


no transverse cracks. 


6.3.2. Stiffness Modulus-Reduced Time Relationships for 
Alberta Asphaltic Concretes 


The stiffness modulus-reduced time relationships and the shift 
factor-temperature relationships for the "at construction" and 34 month 
asphaltic concrete mixtures of the Alberta test project are presented 
in FIGURES VI-5 to VI-9. As shown in these figures, each of the three 


asphaltic concrete mixtures exhibited an increase in estimated stiffness 
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values during the 34 month service period. Such increases can largely 
be attributed to the increase in density and resulting increase in 
volume concentration of aggregate values of the asphaltic concretes 
during this time period as shown in TABLE VI-4. Shields et al~ (1969) 
presented data showing that the most rapid rate of density increase of 


the asphaltic concretes occurred during the first year of service. 


The extent of transverse cracking at various service ages of 
the Alberta test project is summarized in TABLE VI-6. Anderson and 
Shields (1971) attributed the marked increase in cracking in all three 
sections during the third winter to the fact that this winter was con- 
sidered to be one of the most severe in Central and Northern Alberta in 


the past’ 75 years. 


6.4. Test Method Used for Determining Low Temperature 
Tensile Properties of Asphaltic Concrete Cores 


Laboratory procedures most extensively used to characterize the 
tensile properties of asphaltic concrete paving mixtures with respect 
to their low temperature fracture susceptibility have been the direct 
tension test, described by Haas (1968), and the tensile splitting test 
described by Anderson and Hahn (1968) and Christianson (1970). In the 
present study use has been made of the latter test method to determine 
the low temperature tensile properties of the asphaltic concrete cores 
obtained from the test projects. Briefly, the test method involves 
loading cylindrical asphaltic concrete specimens across a diameter in 
a compression frame, and within a chamber maintained at a constant low 


temperature. Output signals from a load cell and series connected 
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linear variable differential transformers attached to opposite faces of 
the specimen were continuously monitored on a two-channel recorder. 
Such recordings enabled the induced tensile stress and strain to be 


calculated at any time during the test. 


Assuming the specimen to exhibit elastic properties and also 
assuming a state of plane stress exists, the strain in the direction of 
a diameter at right angles to the applied load and at the center of the 
face of the specimen is equated to the normal stresses by means of the 


following equation. 


ey = = (0, - voy) . (VI-6) 
where 

€x = horizontal strain 

iS = modulus of elasticity 

Oo, = induced tensile stress 

Sate induced compressive stress 

) = Poisson's ratio. 


When employing loading strips of width less than one-tenth the 
diameter of the specimen the induced tensile and compressive stress at 


the center of the specimen is given by 


0 = 2P O = ePor 
x TDL 7, TDL 
where 
P = applied load 
D = diameter of specimen, and 


L = thickness of specimen. 
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Substituting the above relationships into Iquation (VI-6) 


yields 


fopbulei2e of wor 
sami a i pa) 


In order to calculate the strain due to the tensile stress, 2P/ DL, 
Poisson's ratio must be defined. With the exception of results re- 
ported by Monismith and Secor (1962) and Sayegh (1967), there is little 
published information defining the time and temperature dependence of 
Poisson's ratio, v, of asphaltic concrete paving mixtures. Monismith 
and Secor (1962) found that Poisson's ratio of an asphaltic concrete 
paving mixture varied from approximately 0.35 to 0.50 between tempera- 
tures of 40 and 140 F, respectively. Sayegh (1967) found Poisson's 


ratio to increase from 0.10 at low temperatures, high frequencies and 
small deformations to 0.50 at high temperatures, low frequencies and 
large deformations. In view of these results a value of 0.30 was assumed 
representative for the employed test conditions. Upon substituting 

this value into Equation (VI-7), the strain resulting from the induced 
tensile stress is approximately equal to one half the measured hori- 


zontal deformation over the one inch gauge length employed in the test. 


Due to the limited number of asphaltic concrete cores obtained 
from the test projects, testing conditions were restricted to a nominal 
‘cross-head loading rate of 0.06 in per min and temperatures ranging 
from 20 to -20 F. The test data was analyzed by means of computer 
programs developed by Christianson (1970), and in order to summarize 


much of the data, average tensile stress-strain relationships of cores 
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belonging to a given test population were determined by averaging 


induced tensile stress values at specific strain levels. 


6.4.1. Tensile Properties of Manitoba Asphaltic Concrete Cores 


The stress-strain relationships of asphaltic concrete cores 
obtained from the Manitoba test project are shown in FIGURES VI-10 to 
VI-13, and a summary of the tensile properties of the cores, at failure 
conditions, is given in TABLE VI-7. As observed from this table, each 
of the four asphaltic concretes exhibited a decrease in failure strain 
with a decrease in test temperature. The average failure strain- 
temperature relationships of the cores are shown in FIGURE VI-14. 
Within the range of test temperatures, the asphaltic concrete composed 
of the LV 150-200 penetration grade asphalt exhibited the lowest 
failure strains, while, cores composed of the HV 150-200 penetration 
grade asphalt and those of the LV 300-400 penetration grade asphalt 
yielded somewhat lower failure strains than cores composed of the SC-5 


liquid asphalt. 


A comparison between measured failure strains and crack 
frequencies of pavement sections given in TABLE VI-5 reveals that lowest 
failure strain values are associated with the asphaltic concrete which 
exhibited highest crack frequencies. As mentioned in Section 2.4.4 of 
CHAPTER II, Anderson and Shields (1971) reported that the expected 
failure strain of a paving mixture, determined by the tensile splitting 
test, may be related to the penetration and viscosity of the asphalt 
comprising the paving mixture, and suggested a failure strain value less 


than 10 x 107* in per in, at 0 F, to be indicative of a paving mixture 
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having a high potential for cracking under climatic conditions of 
Western Canada. The results of the present study are consistent with 


these findings and suggested strain value. 


The average tensile failure strength-temperature relationships 
of the Manitoba paving mixtures are shown in FIGURE VI-15. Considering 
the moderately rapid loading rate, together with the low test tempera- 
tures and observed brittle fracture mode of the specimens, rate effects 
on fracture strengths were believed minimized. Data presented by Tons 
and Krokosky (1963) and Haas (1968) indicate that variations in loading 
rate have little influence on the resulting low temperature tensile 
failure strengths of asphaltic concrete paving mixtures. Such results 
suggest that similar fracture strength-temperature relationships as 
shown in FIGURE VI-15 would have been obtained had test conditions not 


been limited to a single loading rate. 


As test temperatures were decreased from 20 F to 0 F the 
tensile failure strength of each of the asphaltic concrete paving 
mixtures increased. However, at a temperature of approximately -5 F, 
those cores composed of the LV and IIV 150-200 penetration grade asphalt 
cements exhibited a decrease in failure strength with a decrease in 
test temperature. Similar decreases in failure strengths with de- 
creasing temperatures have been reported by Tons and Krokosky (1963). 
These investigators attributed such reduced tensile strengths to the 
presence of an uneven distribution of stress in paving mixtures at low 
temperatures. Heukelom (1966) suggested that the tensile strength of 


an asphaltic concrete could be related to the stiffness of the asphalt 
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comprising the paving mixture. This suggested relationship is dis- 


cussed in the following section of this chapter. 


As shown in TABLE VI-7, at test temperatures of 20 and 0 F the 
asphaltic concrete composed of the LV 150-200 penetration grade asphalt 
exhibited the highest average failure stiffness. Failure stiffness 
values of the asphaltic concretes were also computed by means of Equation 
(VI-3), in which stiffness values of the asphalts were determined using 
Van der Poel's nomograph and a loading time corresponding to the time to 
failure of the asphaltic concrete cores. As shown in TABLE VI-7, 
stiffness values determined directly from tests were generally higher 
than stiffness values estimated by this indirect procedure. Differences 
between these two sets of stiffness values may partially be attributed 
to the assumptions involved in the analysis of the test results as well 
as those associated with the indirect stiffness computation method. 

Also, at low test temperatures measured deformations of individual cores 
belonging to a given test population were extremely small and often found 
to be a factor of two from the average of the population. Such measured 


deformations reflect large differences in computed stiffness values. 


6.4.2. Tensile Properties of Alberta Asphaltic Concrete Cores 


The results of the tensile splitting tests performed on 
asphaltic concrete cores obtained from the Alberta test project have 
been reported in detail by Christianson (1970). The average test 
stress-strain relationships of asphaltic concrete cores obtained from 


the test project and composed of asphalts obtained from Suppliers 1, 
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2 and 3 as denoted in TABLE VI-2, are shown in FIGURES VI-16, VI-17 
and VI-18, respectively. The tensile properties of the cores at 


failure are summarized in TABLES VI-8 and VI-9. 


From a comparison of crack frequency associated with the three 
pavement sections and the results of the tensile splitting test, 


Christianson (1970) concluded; 


a) An increase in crack frequency is accompanied by a decrease 
in failure strain and an increase in failure stress and 


failure stiffness. 


b) An increase in density with service life is accompanied 
by a decrease in failure strain, and an increase in 


failure stress, failure stiffness and crack frequency. 


c) In terms of relative comparisons, the pavement section with 
the highest crack frequency (that section which included 
Supplier ieee as shown in TABLE VI-6) also has the 
lowest failure strain and the highest failure stress and 


failure stiffness. 


d) In terms of service life, the low viscosity asphalt cement 
exhibits the greatest change in density, failure strain and 


failure stiffness, and highest crack frequency. 


The average tensile failure strength-temperature relationships 
of the asphaltic concrete cores are shown in FIGURE VI-19. These rela- 
tionships have been employed in the fracture prediction analyses 


presented in the following chapter. 
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O-5.6, [chisuiierStrengths of Mixes and Asphalt Stiffness 


Heukelom (1966) presented data which shows that the tensile 
strength of asphaltic concrete paving mixtures is a function of 
asphalt stiffness. The tensile strength-stiffness relationships for 
two asphaltic concrete paving mixtures, as presented by Heukelom (1966), 
are shown in FIGURE VI-20. Heukelom suggested that the curve marked 
Type I is an example of paving mixes with poor grading and/or compac- 
tion, whereas, the curve marked Type II represents paving mixes with 
better grading and/or compaction. Superimposed on FIGURE VI-20 are the 
tensile fracture strengths, obtained in this investigation, versus 
asphalt stiffness values of the asphalts comprising the mixes. The 
stiffness values were computed by means of Van der Poel's nomograph and 
using loading times corresponding to the fracture times of the asphal- 
tic concrete cores. The scatter in results may partially be attributed 
to the many variables, such as asphalt types and contents, air void 
contents and aggregate gradations, of the asphaltic concrete cores 
included in the study. At high stiffness values, corresponding to low 
test temperatures, the tensile strengths of the cores were generally 
less than those strengths which could be expected assuming the cores 
to be representative of Type I mixes. However, the measured tensile 
strengths are consistent with those given by Ilaas (1968). Using a 
direct tension test Haas reported tensile strengths of asphaltic 
concrete mixes composed of different Western Canadian asphalts of 
approximately 320 psi (22 kg per sq cm) at test temperatures of -10 


and -40 F. 
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6.6. Summary 


In this chapter concepts used to characterize the time and 
temperature dependence of asphaltic concrete paving mixtures incorporated 
within two field test projects in Western Canada have been presented and 
stiffness modulus-reduced time relationships for the paving mixtures have 
been developed. Results of tensile splitting tests on pavement cores 
obtained from the test projects have been summarized and briefly. dis- 
cussed in relation to observed differences in low temperature behavior of 


test control sections within the projects. 
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TABLE VI-5 


TRANSVERSE CRACK FREQUENCIES OF PAVEMENTS 
INCLUDED IN THE MANITOBA TEST PROJECT 
(Young et al., 1969) 


Asphalt Binder 


Structure 


(a) Cracks per mile. 
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TABLE VI-6 


TRANSVERSE CRACK FREQUENCIES 
OF PAVEMENTS INCLUDED IN THE 
ALBERTA TEST PROJECT 
(Anderson and Shields, 1971) 
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S(t), PSI. (log scale) 


Modulus, 


Time, t, (log scale ) 


FIGURE YI-2 EFFECT OF TIME AND TEMPERATURE ON THE 
TIFFNESS MODULUS FOR A THERMORHECLOGICALLY 
SIMPLE MATERIAL (after Monismith et al., 1966) 
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FIGURE VI-8 SHIFT FACTOR VERSUS TEMPERATURE RELATIONSHIPS 
FOR ALBERTA ASPHALTIC CONCRETES AT CONSTRUCTION. 
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FIGURE VI-9 SHIFT FACTOR VERSUS TEMPERATURE RELATIONSHIPS 
FOR ALBERTA ASPHALTIC CONCRETES AT 34 MONTHS. 
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FIGURE VI-I5 AVERAGE TENSILE FAILURE STRENGTH - 


TEMPERATURE RELATIONSHIPS OF MANITOBA 
ASPHALTIC CONCRETE CORES. 
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FIGURE VI-19 AVERAGE TENSILE FAILURE STRENGTH- 


TEMPERATURE RELATIONSHIPS OF 
ALBERTA ASPHALTIC CONCRETE CORES. 
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CHAPTER VII 


STRESS PREDICTION AND LOW TEMPERATURE 


FRACTURE SUSCEPTIBILITY OF ASPHALTIC CONCRETE 


or o RCrOdUuct1 ON 

In recent years various numerical methods have been used to 
predict thermal stresses in asphaltic concrete pavements subjected to. 
assumed temperature regimes. It is only recently that information has 
become available as to the real temperature regimes occurring in pave- 
Hencestructures prior to and at the time’ of crack initiation. In-enirs 
chapter such temperature data, obtained from the field test projects in 
Manitoba and Alberta, together with the stiffness and fracture strength- 
temperature relationships of asphaltic concrete mixes incorporated in 
the test projects and described in CHAPTER VI, is used to assess the 
results of various stress prediction methods. The stress prediction 
methods include pseudo-elastic and viscoelastic beam and slab analyses. 
Theoretical considerations for the formulation of the stress pre- 
diction methods are reviewed and numerical solutions of the stress 


equations are presented. 


Employing a postulated fracture criterion, computed stresses 


are compared with the fracture strength-temperature relationships of 
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the asphaltic concrete cores obtained from the individual field 

test projects to predict times of fracture. These times of fracture 
are then compared with times of observed low temperature cracking of 
the corresponding pavements in the field projects. Such comparisons 
provide a means of assessing the stress predictive methods in relation 
to observed differences in the low temperature fracture susceptibility 


of asphaltic concrete paving mixtures included in the test projects. 


7.2. Stress Predictive Methods 

A review of a number of different methods used to compute 
thermal stresses in asphaltic concrete pavements has been presented by 
Haas and Topper (1969). These analyses may be classified according to 
the boundary constraints on a pavement element. Considering an element 
of pavement as shown in FIGURE VII-1, the term ''beam analysis" implies 


the conditions 


ey To 0 Cy eye EVEI-1) 


& = ey = OU Cee ane a se (VII-2) 


For pavement structures the beam assumptions approximate conditions at 
the edge of a pavement while the slab assumptions are more severe than 


those which would be anticipated in the center of the pavement. 


Five analyses have been used in this study for stress compu- 


tations. These analyses are designated as: _ 


1) Pseudo-elastic beam. 
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2) Approximate pseudo-elastic slab. 
3) Viscoelastic slab. 
4) Viscoelastic beam. 


5) Approximate viscoelastic slab. 


A description of these analyses follows. 


7.2.1. Pseudo-Elastic Beam Analysis 
The stress equation for the pseudo-elastic beam analysis is 


t 


Halts = -f S(At,T)a, (T)AT(t) (VII-3) 


6 


where S(At,T) is the time and temperature dependent stiffness modulus. 
If the coefficient of the thermal expansion, a,, and S are independent 
of time and temperature this equation reduces to the simple elastic 
thermal stress equation. To evaluate stress numerically Equation 


(VII-3) is written as 


Gy (ti) = Gy(tq-4) * (49y)5 (VII-4) 
where t; 
CRO) -f S(At ,T)a,(T)dT(t) (Vil=5) 
aval 
ee 


Hills and Brien (1966) suggested that the stiffness modulus in 
Equation (VII-5) be determined for a loading time (At) equal to the time 
Incerva leer Perge Se over which the change of stress (Ao,); is computed, 
and at a mean temperature during this time interval. McLeod (1969) 


recommended that the stiffness be evaluated at a loading time (At) of 
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20,000 seconds. Fromm and Phang (1971) used stiffness values deter- 


mined at a loading time of 10,000 seconds. 


In this study the pseudo-elastic beam stress increments were 
computed using the trapezoidal rule with stiffnesses evaluated at the 
temperatures for times t; and t;_,, and a loading time (At) equal to 
the time step (t; -"ti_;). With the assumption that the coefficient of 
thermal expansion is temperature aereaen the approximation to 


Equation (VII-5) then becomes 
Mo 

(Ae Tee bea * Sic Ty) (VII-6) 

where the subscripts i and i-1 indicate quantities at times t, and ae 


Naz 


respectively. 


7.2.2. Pseudo-Elastic Slab Analysis 

A pseudo-elastic slab analysis may be evaluated to satisfy con- 
ditions (VII-2) by inserting the factor 1/[1-v(T,t)] within the integral 
sign of Equation (VII-3). As mentioned in Section 6.4 of CHAPTER VI, 
with the exception of results reported by Monismith et al. (1962) and 
Sayegh (1967), there is little published information defining the time 
and temperature dependence of Poisson's ratio, v, of asphaltic concrete 
paving mixtures. The approximate pseudo-elastic slab stress was there- 
fore evaluated by assuming Poisson's ratio to be a constant and equal 
to 0.30. The stress obtained from Equation (VII-4) was therefore 
simply multiplied by the factor 1/(1-v) to obtain an approximate 


pseudo-elastic slab stress. 
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7.2.3. Viscoelastic Slab Analysis 

Humphreys and Martin (1963) proposed a method of determining 
thermal stresses in a slab with temperature dependent viscoelastic 
properties. Using the time-temperature equivalence hypothesis of a 
thermorheologically simple material described in Section 6.3 of 
CHAPTER VI, these authors expressed deviatoric and hydrostatic stresses 
(Si j and o) in terms of the deviatoric and volumetric strains (i; and ¢) 
by the following relationships 


1E 
2 if ee — fe, jt! de (VII-7) 


=co 


a 
1J 


t 
i G,(E - &") = [eters ageit! Jdtt wqvires) 


= co 


a(t) 


where © is the reduced time introduced in the discussion of material 
characterization in Section 6.3 of CHAPTER VI, G, and G, are deviatoric 
and volumetric relaxation moduli and 6 is a pseudo-temperature defined 


as 
T(t) 
arena? (VII-9) 


Tet) 


Assuming the dilatational response to stress remains elastic, Humphreys 
and Martin (1963) presented a numerical technique for evaluating an 
effective modulus of a slab, R(&), from stiffness modulus-reduced time 


relationships such as those presented in CHAPTER VI. For a constant 
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coefficient of thermal expansion, a(T') = a Equations (VII-7), (VII-8) 


Oo? 
and (VII-9) yield 
t 


oF -50, f R[& - &']dT(t') (VII-10) 


fe) 
Monismith et al. (1965) used the trapezoidal rule to evaluate 


Equation (VII-10) numerically, according to the expression 


N 
o(t) = - *%o 5 [R(E-€4) + RCE-Eq_,)](Ty-Ti-,)— (VIT-11) 
2 


i=] 
where N is the number of time steps from 0 to t. Equation (VII-11) is 


the equation used in the present study to compute viscoelastic slab 


Stresses. 


7.2.4. Viscoelastic Beam Analysis 

The stresses in a viscoelastic beam are conceptually simpler 
than those in a slab and follow directly from the superposition 
principle once the time-temperature equivalency hypothesis has been 
introduced. Thus, for a constant a,, we may write 


t 
a(t) = a, f S(— - &*)dT(t') (VII-12) 
to 
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7.2.5. Approximate Viscoelastic Slab Analysis 


Assuming that Poisson's ratio remains constant an approximate 
estimate of stresses in a viscoelastic slab may be obtained by multi- 
plying the stresses predicted from Equation (VII-13) by the factor 
1/(1-v). This assumption for v, which is purely for analytical 
Simplicity, yields results considerably different than those obtained 


from Equation (VII-11). [Humphreys and Martin (1963). ] 


7.3. Computation of Thermal Stresses 

Thermal stress predictions based on the numerical solutions of 
the pseudo-elastic and viscoelastic stress equations were obtained by 
computer operations. Details of the computer programs developed for 
these analyses are contained in APPENDICES B and C, respectively. The 
primary input variables for the programs and the numerical integration 
procedures associated with the analyses are described in the following 


subsections. 


Vea Oe Input Variables for Thermal Stress Computations 


The primary inputs to the computer programs include 


1) a stiffness modulus-reduced time relationship for the 


asphaltic concrete, 


2) the corresponding shift factor-temperature relationship, 


3) a temperature field, 
4) a coefficient of thermal expansion, and 


5) Poisson's ratio, for slab analyses. 
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The pavement test sections included in the present study have 
been described in Section 6.2 of CHAPTER VI. Derived stiffness modulus - 
reduced time relationships and shift factor-temperature relationships of 
asphaltic concretes associated with these test sections are shown in 


FIGURES VI-3 to VI-9 of the same chapter. 


For computer operations, stiffness values are specified at 
small increments of reduced time and a linear relationship is assumed 
between the logarithm of stiffness and the logarithm of reduced time. 


Simvlarly, shirt factors, are specified at small temperature 


aps 
increments and a linear relationship assumed between the logarithm of 


ap and the input temperatures. 


For Structure A of the Alberta test project monitored tempera- 
tures during the First winter of Service were used for stress compu- 
tations while during the third winter, the winter of 1968-69, stresses 
were computed using a temperature history predicted by means of the 
analytical model described in CHAPTER V. For Structures B, C and D, 
the temperature histories Hatiabie for stress computations were those 


transcribed at 2 hour intervals from the continuously monitored field 


temperatures of the Manitoba test project. 


The results of a number of investigations related to defining 
the coefficient of thermal Brvenoien of asphaltic concrete have been 
reviewed in Section 3.3.3 of CHAPTER III. In the present study the 
coefficient of thermal expansion of all the asphaltic concrete paving 


mixtures was assumed temperature independent and equal to 1.5 x 107° per 


degree Fahrenheit. For slab analyses, Poisson's ratio of the asphaltic 
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concretes was assumed constant and equal to a value of 0.30. 


7.3.2. Numerical Integration for Pseudo-Elastic Analyses 

One of the difficulties in applying pseudo-elastic analyses 
based on the Hills and Brien approach of 1966 is that the predicted 
stress is dependent on the time interval used in the numerical evalua- 
tion. A comparison of stresses computed by Equations (VII-4) and 
(VII-5) for 15 minute and 2 hour intervals using the characterizing 
stiffness properties of the low viscosity 150-200 penetration eee, 
asphalt pavement of the Manitoba test project is shown in FIGURE VII-2. 
The maximum computed stress using the 15 minute time increment is 
approximately 50 per cent greater than that computed using the 2 hour 
time increment. This difference is primarily due to the change of 
Sean with loading time, since at a given temperature the stiffness 
of an asphaltic concrete increases with decreasing loading times. The 
previously mentioned methods of McLeod (1969) and Fromm and Phang (1971) 
eliminate the dependency on the time interval of numerical integration 
by eee eyinn loading times. , However, the results of the stress 
computations are directly influenced by the loading time specified. 


For these pseudo-elastic analyses stresses were computed using a 2 hour 


time increment. 


In addition to-the dependency of the pseudo-elastic stresses on 
assumed loading time, such methods do not make allowance for the relax- 
ation of stresses subsequent to the time interval in which stresses are 
Although thermal stresses within pavements are unknown, the 


computed. 


assumption of a constant stress over extended time periods is not 
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compatible with the known viscous bchavior of bituminous matcrials. A 
viscoelastic model of material characteristics would therefore be 


expected to produce more reliable results. 


7.3.3. Numerical Integration for Viscoelastic Analyses 

The modulus function, R(&), appearing in Equations (VII-10) and 
(VII-11) was evaluated numerically, for the asphalt concretes under 
consideration, by the method proposed by Humphreys and Martin (1963). 
Details of a computer program developed for this purpose are contained 
in APPENDIX D. A study of convergence of stress results with different 
time increments of numerical integration and the computation times on 
an IBM-360 computer are summarized in TABLE VII-1. The stresses shown 
are those predicted using the viscoelastic beam analysis and the 
material properties of the low viscosity 150-200 asphalt pavement in- 
corporated in the Manitoba test project and subjected to the temperature 
history shown in FIGURE VII-2. Since the viscoelastic stress equations 
require numerical computation over the entire previous temperature 
history to which the asphaltic concrete is subjected, long computer 
runs and large storage requirements are involved if stresses are to be 
predicted over extended periods of timc. Because the most recent 
temperature history is the most significant for stress evaluation, the 


following scheme of variable time increments was adopted: 
a) For times up to (t-24) hours a time increment of 2 hours 
was employed. 
b) For the following 22 hours a time increment of 15 minutes 


was employed. 
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c) For the 2 hours immediately prior to t, a time increment 


of 5 minutes was employed. 


This scheme, illustrated in FIGURE VII-3, made it reasonably economical 
to evaluate stresses over extended periods of time and, as shown in 
TABLE VII-1, the resulting accuracy is comparable to that for a 
constant 15 minute time increment. Details of the computer program 
which incorporates these variable time increments are presented in 
APPENDIX C, while a similar program, but one in which a constant time 
increment of numerical integration may be specified is contained in 


APPENDIX E. 


It should be noted that further economy in computing time could 
be obtained by deriving closed form expressions for the characterizing 


properties of the asphalt concretes. 


7.4. Computed Stresses and a Cracking Criterion 


Stresses at 1/2 inch depth intervals throughout each of the 
asphaltic concrete surfaces included in the present study were computed 
by each of Equations (VII-6), (VII-11) and (VII-13). Stresses for the 
approximate slab analyses were determined by multiplying the beam re- 
sults by 1/(l-v), where, as previously mentioned, Poisson's ratio was 


assumed equal to 0.30. 


For all pseudo-elastic analyses stresses were computed using a 
2 hour time increment, a time increment corresponding to the time inter- 
val of temperature input. Viscoelastic analyses were evaluated using 


the combination of 120, 15 and 5 minute increments previously described. 


fso imonoo> Vldenovsst? st Sbsm ani 
ni aviod2 28 | bag SMELT 10 -bot od : ) 29eZ29ETe - OT 

s 10% tadt Oo olden qtiye ak ‘arpa peor et 
MhLVSTY FSIUgOS. SAF Lo abana _rvore ont — odonin 2% 

i botne2srq srs 27 0SMesOnE omit Bidapeos ‘ees ae 
anis. tust2noo: & Ww fiw of sae ud (aietyone avita$e BS ‘ota’ 
ni honisinos ef pakeedoue ed Xai aobdengozet tao tromscr 2 


i i me while 


biues omiz gnisyemoo me monoa$ i tois beson od incl 
onizrivetosts do. ont ton OL 829 TQ? ee: gniviteb Va vn + 
mee “— 


: fa. ic 
“ST. msod orld gir .ig ba Lume ed ben imres: 


2g oljst e'noas fot Pini i 


188 


The analyses were initiated in late autumn or early winter when the 
temperature data indicated a uniform temperature of 30 F throughout the 


pavement surfaces, at which time the stress was assumed equal to zero. 


Thermally induced stresses within the low viscosity 150-200 
asphalt pavement of Structure D, for the period December 30, 1967 to 
January 4, 1968, computed by the viscoelastic beam analysis, are shown 
in FIGURE VII-4. It is apparent that the maximum tensile stress occurs 
at or near the upper surface and that a relatively high stress gradient 
exists during periods of maximum stress. Since this stress gradient is 
typical, the stress at the 1/2 inch depth was adopted as a measure of 
the severity of the stress condition to which the asphaltic concretes 


were subjected. 


The computed stresses at the 1/2 inch depth resulting from the 
pseudo-elastic beam, viscoelastic slab and viscoelastic beam analyses 
are shown in FIGURES VII-5 to VII-14, for the initial time periods of 
critical stress within asphalt pavement surfaces of the Manitoba test 
project. In addition, the fracture strength for the corresponding 
temperature, obtained from FIGURE VI-15 of CHAPTER VI,is also shown on 
these figures. Similar stress-time plots were determined for the 
asphaltic concretes incorporated in the Alberta test project. FIGURES 
VII-15, VII-16 and VII-17, show stresses computed at the 1/2 inch 
depth of the three Alberta test sections. The stresses shown are those 
computed using the pseudo-elastic beam analysis together with the pre- 


dicted third winter of service temperature history. 


eat ™: ») 


| i 
‘ie 


OUS-02t yr TPesB.iv wel, ena ae Roe eS 
os Yoel 08 aadmeosd bokteq ecton srutourye, ‘1 


ra 
gWorl2 STB. 2teyPEens acsdaizentsbgeky tit Kd banwqnie: nae Pa! 
; i Tool i 
21990 2eotte olLanes miitzem add anata ares ai ai ane 
WAN 4) 
tietbor ezasité gin Yleviteles 6 ia bine ‘eosin i al wig 
r ai bt 


2% ineitheiy <¢eoite erat oon 22 eon muta 20 seataei, 
to sivensm 6 es bosqobs: 28M rand soni s\ “ 76 eeonse aa 

7 ig ii 
RySe. 


esteraitos sitledqes arts Aye Gt not fino dvorsze eit to ¥ 


tee 
; 


edt mort anitivest Aigoh dont ae 09 agneorse bosucn ot al 
exisne mood > ite kedoe Ne ‘bets. geese einod ait alae | 


$297 BqODTEM, ott, *o aanttus 


naxbaogesrres oat 1 ghost 


Qo nwode O2fs ef IV satan Bb ataty: 
," re 
od? 16t honinrredely eer cine 


eMIUT .s99bong $e0t shad ofa oe 


dont: SX0 odd a pens aoaeoria! wort 


ri Shir hyp CGT 


For the purpose of comparing times of predicted and observed 
cracking the following criterion is postulated. Cracking is assumed to 
occur when the computed stress at the 1/2 inch depth equals or exceeds 
the fracture strength of the asphaltic concrete for the corresponding 
temperature. According to this criterion the predicted time of initial 
cracking of the asphaltic concretes may be obtained for each of the 
stress analyses by locating the first point of intersection of the 
fracture strengths with the computed stresses on the stress-time plots. 
It is also postulated that the degree of cracking is related to the 


amount by which the computed stresses exceed the fracture strength. 


Va fo" Comparison of Observed and Predicted Low 
Temperature Performance 


At the Manitoba test site crack detection circuits were used 
in an attempt to define initial times of fracture. Details of these 
circuits have been described by Deme and Fisher (1968). With the 
exception of SC-5 test sections and the asphaltic concrete pavement of 
Structure D which incorporated the high viscosity asphalt, crack 
detection circuits were included in all surfaces of the ten pavement 
sections under consideration. Unfortunately, the circuits did not 
perform entirely satisfactorily since, in some sections, visual crack 
surveys revealed cracking prior to the time of cracking indicated by 
the circuits. In the following comparisons initial cracking of the 
pavement surfaces was considered to have occurred when either the crack 


sensors indicated cracking or cracks were sighted visually, whichever 


came first. 
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During the first two winters that the Alberta test project was 
in service continuity circuits, consisting of a high conductivity silver 
base paint and placed in a longitudinal saw cut on the paved shoulder, 
enabled temperature conditions conducive to thermal fracture of the 
pavement sections to be defined. Details of these circuits have been 


described by Shields etal. (1969). 


7.5.1. Manitoba Test Sections 

A tabulation of the times in each day at which computed stresses 
intersect the fracture strengths of the asphaltic concretes included 
in the Manitoba test project is given in TABLES VII-2 and VII-3. 
A summary of the results obtained from the crack sensors and visual 
surveys, reported by Deme and Fisher (1968) and Young et al. (1969) is 
also included in TABLES VII-2 and VII-3. The following observations 
can be made on the comparison of predicted and observed behavior shown 


in these tables. 


7.5.1.1. Structure D, Low Viscosity 150-200 Asphalt Pavement 


Crack sensors indicated initial cracking at 2300 hours on the 
30th day of December, 1967 (23/30/12). The pseudo-elastic slab and 
approximate viscoelastic slab analyses indicated initial cracking at 
the same time as that recorded in the field. The viscoelastic slab 
analysis predicted initial cracking 6 hours earlier while the pseudo- 
elastic beam analysis predicted fracture 8 hours later. The visco- 
elastic beam analysis did not predict cracking during this temperature 


cycle. The plot in FIGURE VII-5 compares the computed stress histories. 
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All analyses predicted cracking from 4 to 8 hours earlier than that 
indicated by the crack sensors on the night of January 3-4. Cracking 
predicted during the intervening time period by all analyses, with the 
exception of the viscoelastic beam analysis, was not recorded by the 


crack sensors or observed visually. 


73.1.2. Structure B, Low Viscosity 300-400 Asphalt Pavement 


Crack sensors indicated initial cracking at 2:50/4/1. Again 
the pseudo-elastic slab and the approximate viscoelastic slab analyses 
indicated initial cracking within an hour of the recorded initial 
fracture time. The viscoelastic slab analysis predicted cracking 4 
hours earlier while the pseudo-elastic beam analysis effectively 
predicted cracking 5 hours later. The viscoelastic beam does not 
“predict cracking in this temperature cycle. The computed stress 
histories for this pavement are shown in FIGURE VII-6. All analyses 
with the exception of the viscoelastic beam analysis also indicated 
cracking during the night of January 5, 1968. Such predicted fracture 
is consistent with the increase in cracking of this test section de- 
termined by the January 6, 1968 visual survey, reported by Young et al. 


(1969). 


7.5.1.3. Structure B, Low Viscosity 150-200 Asphalt Pavement 


Visual surveys indicated initial cracking had occurred on 
December 31, 1967. All analyses predicted cracking in the late after- 
noon and early evening of December 30, 1967. Therefore, there is no 
incompatibility between the times of predicted and observed initial 


cracking. As shown in FIGURE VII-7, the computed stresses continued to 
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D2 
exceed the fracture strength on the succeeding days when the majority 


of cracking was observed. 


7.9.1.4. Structure C, Low Viscosity 150-200 Asphalt Pavement 


Visual surveys indicated initial cracking occurred on January 1, 
1968, while the majority of cracking occurred on the sixth and seventh 
days of the same month. All analyses, except the viscoelastic beam 
analysis, predicted initial cracking at sometime during the night of 
December 30, 1967, one day before cracking was observed. The stress 
predictions are shown in FIGURE VII-8 and, as observed from TABLE VII-2, 
all analyses predicted cracking over the time period during which the 


majority of cracking occurred. 


7.5.1.5. Structure B, High Viscosity 150-200 Asphalt Pavement 


No fracture was observed in this section of the test project. 
However, all stress analyses except the viscoelastic beam analysis 
predicted cracking at sometime during the period January 3 to January 6, 
1968. The predicted stress histories are shown in FIGURE VII-9. As 
shown in this figure and TABLE VII-3, stresses computed using the 
elastic beam analysis slightly exceeded the average tensile strength 
of the asphaltic concrete only once during the 1967-68 winter when a 


near minimum 1/2 inch depth temperature of -36 F occurred. 


7.5.1.6. Other Manitoba Test Sections 


Results of stress analyses for the five remaining Manitoba 
pavement test sections under consideration are shown in FICURES VII-10 


to VII-14. In each of these five sections no fracture was predicted. 
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This result conforms to field observations Summarized in TABLE VII-3. 


7.5.2. Alberta Test Sections 

Crack frequencies within the three pavement sections at the 
Alberta test site have been presented in TABLE VI-6 of CHAPTER VI. 
During the first two winters of service the majority of cracking 
occurred in the pavement section containing Asphalt Supply No. 1, and 
the temperature and.continuity circuit systems revealed that initial 
cracking occurred at temperatures summarized in TABLE VII-4. During the 
third winter, such systems were not operative and, therefore, fracture 


times and temperatures of each pavement section were not recorded. 


7.5.2.1. Results of Stress Analyses for First Winter of Service 

The results of the stress analyses for the first winter of 
service indicated no fracture for each of the three test sections. The 
maximum computed stress, determined by the viscoelastic slab analysis, 
and associated with the test section containing Asphalt Supply No. 1 was 
84 psi. This stress is considerably less than the corresponding 
measured tensile failure strength of the asphalt concrete shown in 
FIGURE VI-19 of CHAPTER VI. While this result is not consistent with 
the observed transverse cracking within this particular test section, 
the result conforms to field observations associated with the test 


sections incorporating asphalts obtained from Suppliers 2 and 3. 


7.5.2.2. Results of Stress Analyses for Third Winter of Service 
USS SE dial Mada een nde es cin Tesla SE SOs 
During January of the third winter, the winter of 1968-69, an 


extreme cold period of over 25 days was encountered. During this time 
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the predicted 1/2 inch depth temperatures at the test site never 
exceeded 0 F. Although the times of cracking were not recorded, the 
results of the stress analyses were compatible with observed low temper- 
ature performance, since all stress analyses predicted fracture of the 
test sections. Typical stress-time plots for the three test sections, 
obtained using the pseudo-elastic beam analysis, are shown in FIGURES 
VII-15, VII-16 and VII-17. The predicted temperatures of crack initi- 


ation of each section are summarized in TABLE VII-4. 


7.6. Evaluation of Results 

Considering the variability of the properties of asphaltic 
concrete mixtures and the approximations involved in the study of low 
temperature pavement behavior, a high degree of accuracy should not be 
expected of a predictive analyses. Before attempting to draw con- 
clusions from the comparisons of observed and predicted performance the 
nature of the approximations and assumptions involved are briefly 


reviewed. 


Ta On lis Approximations and Assumptions 


(a) Material Characterization: It was assumed that the 
stiffness of the asphaltic concrete mixtures could be determined by 
indirect methods. The coefficient of thermal expansion was assumed to 
be independent of temperature and equal to 1.5 x Waa per degree 
Fahrenheit. For the approximate slab analyses Poisson's ratio was 
assumed to be 0.30. The fracture strengths of the pavements were 


assumed to be represented by the strength-temperature relationships 
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determined from the average of laboratory tests for pavement cores 


obtained from the test project. 


(b) Nature of Loading: It was assumed that pavement cracking 
was a result of thermally induced stress and that all other loading 
mechanisms could be neglected. It was also assumed that interaction 
between the asphaltic concrete and the substructure did not influence 


thermally induced stresses. 


(c) Fracture Criterion: The stress fracture criterion 


previously described was assumed valid. 


In addition to the assumptions outlined above, each method of 
stress analyses. assumes the validity of the corresponding constitutive 
equation and that the error of the numerical integration procedure is 


within acceptable limits. 


7.6.2, Comparison of Stress Analyses 


The’ suitability of, any of the stress, analysis procedures” “tor 
use with the proposed fracture criterion, may be assessed (a) by com- 
paring times of predicted initial cracking with observed times of 
initial cracking, and (b) by correlating the time periods, during which 
subsequent cracking was observed with the time periods in which the 
predicted stress exceeded the fracture criterion employed. Strictly 
speaking, the stresses computed after initial cracking cannot be 


compared directly with the fracture criterion because the assumption 


that c, = 0 is no longer valid. 
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Using the time period correlation method, method (b) above, it 
appears that the viscoelastic beam analysis underestimates stresses for 
the fracture criterion employed. (See Structure D, LV. 150-200, TABLE 
VII-2 and FIGURE VII-5; Structure B, LV 300-400, TABLE VII-3 and 


FIGURE VII-6.) 


It ‘alsovappears that the viscoelastic slab ‘analysis overésti- 
mates stresses for the fracture criterion employed. (See Structure B, 
HV 150-200, TABLE VII-3 and FIGURE VII-9.) The viscoelastic beam and 
Slab analyses therefore appear to yield lower and upper bounds re- 
spectively, on computed stresses for the purpose of correlating with 
observed pavement behavior. The pseudo-elastic beam and approximate 
viscoelastic slab analyses appear to yield reasonable intermediate 
values of stresses for correlation purposes. The approximate pseudo- 
elastic slab ereivcie yields stresses similar to those of the approxi- 
mate viscoelastic slab analysis. Since none of these analyses can be 
rigorously justified on a theoretical basis the author has selected 


the simplest, namely, the pseudo-elastic beam analysis for a comparison 


with observed results. 

TABLE VII-S shows a correlation of times of initial predicted 
fracture, using the pseudo-elastic beam analysis adopted for this 
study, with the observed cracking of the sixteen different combinations 
of structure, materials and temperature history incorporated in the 
Study. In fourteen of these different combinations low temperature 


cracking was predicted at times approximating times of recorded and/or 


observed fracture. In the author's opinion this correlation is 
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sufficiently good to allow the low temperature fracture susceptibility 
of asphaltic concrete pavements to be assessed when subjected to 


thermal regimes similar to those experienced in Western Canada. 


Wands Summary 


Five different methods of stress computation were employed to 
evaluate thermal stresses in a number of different test sections 
incorporated in two test projects in Western Canada. Stresses within 
the particular asphaltic concretes were evaluated using stiffness 
properties presented in CHAPTER VI. Fracture strength-temperature 
relationships of asphaltic concrete cores obtained from the test 
projects were used for predicting times and temperatures of fracture. 
Observed cracking in the test projects was compared to that predicted 


by the various stress analyses using a postulated fracture criterion. 


To the author's knowledge, these test projects provide, for the 
first time, detailed records of field temperature data for which field 
observations of low temperature pavement cracking are also available. 
Such field data provides a basis on which to realistically assess 
various stress predictive methods. The correlation obtained between 
observed and predicted performance appears sufficiently good to 
indicate that reliable predictive analyses can be developed. The data 
currently available suggests that, for the present, a Suitable pseudo- 


elastic beam analysis can yield reasonable results. 
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TABLE VII-4 


TEMPERATURES DURING CRACK INITIATION 
AT ALBERTA TEST PROJECT 


PIyst . Third Winter 
Winter (Predicted) 
(Recorded) 
(F) 


Supplier 
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TABLE VII-5 


CORRELATION OF PREDICTED AND OBSERVED INITIAL TIMES OF 
FRACTURE USING THE PSEUDO-ELASTIC BEAM ANALYSIS 


Asphalt 


LV 150/200 Predicted fracture 8 hours 
after that recorded by 
— i 


crack sensors. 
vain 
| 
No. 1] 


Predicted fracture 5 hours 
after that recorded by 
crack sensors. 


Predicted fracture the 
evening of December 30, 
1967 - visual Surveys 
indicated cracking the 
following day. 


Predicted fracture the 
morning of December 31, 
1967 - visual surveys 
indicated cracking the 
following d 


A (at 34 months) Asphalt Supply : Fracture predicted and 
A (at 34 months) Asphalt Supply No. 2 | observed, however time of 
A (at 34 months) Asphalt Supply No. 3 | fracture was not recorded. 


HV 150/200 
LV 300/400 

ee No fracture predicted or 
HV 150/200 observed. 


construction) 
construction) 


GO PO 


Asphalt Supply No. 
Asphalt Supply No. 


Fracture predicted near 
seasonal minimum 
temperatures - no cracking 
observed. 


HV 150/200 


No fracture predicted - 


Asphalt Supply No. 1 
cracking observed. 


construction) 
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CHAPTER VIII 


GUIDELINES FOR REDUCING TRANSVERSE CRACKING 


OF ASPHALTIC CONCRETE PAVEMENTS 


8.1. Introduction 

During the past decade extensive data has been collected from 
field and laboratory studies related to the low temperature behavior 
of flexible pavements. From this data several approaches for mini- 
mizing transverse cracking within new pavements have been formulated. 
This chapter contains a brief review of these approaches and, 
employing the results of the stress prediction analyses contained in 
CHAPTER VII, a design guideline for reducing transverse cracking is 
developed. Comparisons are made between limiting low temperature 
design criteria specified by various agencies and individuals with 


Similar criteria derived from the proposed design guideline. 


8.2. Previous Design Approaches 


Many of the developed design approaches associated with the low 
temperature cracking problem have been summarized by Haas et al. (1970). 
In view of this summary, the content of this section is limited to a 
brief review of suggested approaches which incorporate the following 
concepts. 

1) setting limiting specifications on the asphalts, 
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2) setting limiting asphalt or asphaltic concrete mix 


stiffness values, 


3) predicting fracture temperatures of asphaltic concrete 


mixes, and 


4) predicting crack frequencies of asphaltic concrete 


pavements. 


8.2.1. Setting Asphalt Specification Limits 


From the results of field studies reviewed in CHAPTER II, the 
low temperature fracture susceptibility of an asphaltic concrete mix- 
ture appears primarily dependent on the asphalt used. Recognizing this 
dependency various highway departments modified their asphalt specifi- 
cations in an attempt to reduce the occurrence of transverse cracking. 
These modifications included setting higher penetration values and/or 
implementing minimum viscosity requirements. Typical of such modifi- 
cations were those introduced during the mid-1960's by the provinces of 


Saskatchewan and Alberta. 


Prior to 1963, the Saskatchewan Department of Highways asphalt 
specification called for a 150-200 penetration grade asphalt. Subse- 
quent modifications led to the present specification of minimum and 
maximum penetration and viscosity requirements at 77 F and 140 F, 
respectively. Until 1967, asphalt cements for use in highway con- 
struction in the Province of Alberta were graded according to 
penetzation=at 77 F. This grading resulted in asphalt cements of a 
given penetration grade exhibiting large differences in viscosity at 


140 F. In 1967, the Alberta Department of Highways introduced their 


current asphalt specification which incorporates a minimum viscosity 
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requirement at 140 F. The grade most commonly used calls for a 
minimum penetration of 250 at 77 F and a minimum vaScosity Of 275 
poises at 140 F. Anderson and Shields (1971) reported that this 
specification has resulted in a more uniform product from several 
distinct asphalt sources, the elimination of tender mixes and an 
apparent reduction in low temperature transverse cracking under normal 


winter conditions. 


While highway agencies have reported a reduction in transverse 
cracking through the use of softer asphalts, the influence of these 
asphalts on traffic load associated distress, such as rutting and 


fatigue, has not been fully evaluated. 


S.2625 Limiting Asphalt or Mix Stiffness Values 

From data collected from various field and laboratory investi- 
gations McLeod (1969) concluded that transverse cracking will occur if 
the stiffness of the asphaltic concrete paving mixture, at the minimum 
service temperature encountered, falls within the range of 1 x 10® to 
2 x 10© psi. McLeod suggested that these values were applicable for a 
dense, well graded mix and asphalt stiffness values derived from 
McLeod's 1969 modification of Van der Poel's stiffness nomograph at a 
loading time of 20,000 seconds. TABLE VII-1 shows McLeod's tentative 
design guide for maximum mix stiffness values associated with various 
temperatures. 

Using nomographic procedures for asphalt stiffness computations, 


Fromm and Phang (1971) determined the stiffness of asphalts comprising 


various mixtures known to exhibit different low temperature crack 
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frequencics. From these calculations the authors suggested a design 


approach that limits asphalt stiffness to 20,000 psi at a loading time 


of 10,000 seconds. 


In the spring of 1972 the British Columbia Department of 
Highways introduced an asphalt specification which incorporates a 
limiting asphalt stiffness. The limiting stiffness, 30,000 psi at a 
specified temperature and at a loading time of 7,200 seconds, is 
computed from Van der Poel's nomograph, FIGURE VI-1, using asphalt 


properties evaluated in a manner proposed by Heukelom (1969). 


Beceow Prediction of Asphaltic Concrete Fracture Temperature 


Basically, this approach is a variation of the limiting 
stiffness approach. Using stiffness values, thermally induced stresses 
within an asphaltic concrete paving mixture subjected to a given temper- 
ature history are predicted. The predicted stresses are compared with 
the fracture strength-temperature relationship of the mix. Such a 
comparison provides a means of predicting the fracture temperature of 
the mix. This approach, initially proposed by Hills and Brien (1966), 


was used in CHAPTER VII to assess the results of various stress 


predictive methods. 


8.2.4. Prediction of Crack Frequency 

Recently, Hajek (1971) developed a mathematical model for 
predicting crack frequency of asphaltic concrete pavements as a 
function of several of the variables described in CHAPTER II. The 


model is based on field data from the provinces of Ontario and Manitoba 
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and is of the following form. 


I = function (s,t,a,d,m) + E 

where 

I = cracking index (the sum of the full plus 
one-half of the half transverse cracks 
per 500 feet of two-lane roadway) , 

Ss = © stifiness modulus of the original” asphalt 
cement, kg per sq cm, according to 
McLeod (1969), 

t = thickness of asphaltic concrete pavement 
layer, inches, 

‘aa = age of pavement, years, 

d = subgrade soil type, classified as clay, 
loam or sand, . 

m = winter design temperature, C, and 

E.s.=»\ standand. error of. estimate. of 1. 


Ilajek concluded that the model enables crack frequency to be predicted 
with a reasonable degree of confidence and has the capacity to provide 
an immediate engineering guide for preventing or controlling low 


temperature cracking in Canada. 


8.3. Proposed Design Guideline 

In CHAPTER VII comparisons have been made between predicted and 
observed fracture times and temperatures of various asphaltic concrete 
For these comparisons results of different stress 


paving mixtures. 


predictive methods were used. The correlation between predicted and 


observed fracture conditions, shown in TABLE VII-5 of CHAPTER VII and 
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obtained using the pseudo-elastic beam analysis, suggests that such an 
analysis provides a reliable means of assessing the low temperature 
fracture susceptibility of an asphaltic concrete paving mixture sub- 

. jected to a particular temperature history. In view of this corre- 
lation, a design guideline for reducing the occurrence of transverse 
cracks, caused by thermally induced stresses exceeding the low tempera- 
ture tensile strengths of asphaltic concrete paving mixtures, is 
proposed. The development of the guideline involved the following 


computations. 


1) Using the SO. archers beam analysis and a loading time 
of 7,200 seconds, thermal stresses within each of the 
aavnalere concrete paving.mixtures, whose stiffness 
characteristics have been presented in CHAPTER VI, were 
oeiitad when subjected to five different temperature 

histories. These temperature histories included, the three 
recorded at the Manitoba test project during the 1967-68 
eeieews te inedaniced asphaltic concrete pavement tempera- 
tures at the Atenes test patie during the winter of 
1966-67 and the predicted third winter temperatures at 
hus ctesit sites thorrall stress predictions the coefficient 
of thermal expansion of Hak paving mirente was assumed 
constant and equal to a value of 1.5 x 107° per degree 


Fahrenheit. 


2) At temperatures of 0 F and less, and using 5 F temperature 


increments, maximum thermally induced stresses within each 
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of the asphaltic concrete pavements were transcribed from 
the results of the stress analyses, At corresponding 
temperatures and employing a loading time of 7,200 seconds 
stiffness values of the paving mixtures were computed. For 
these computations the method adopted in CHAPTER VII to de- 
fine the time and temperature dependence of the mixes was 
foeas The maximum predicted thermally induced stress 
values were then plotted as a function of mix stiffness 
values.. Individual stress and stiffness values of the low 
viscosity 150-200 asphalt pavement at the Manitoba test 
project are shown in FIGURE VIII-1, while, the bands in this 
figure represent the limits of the thermally induced stress- 
stiffness relationships of all asphaltic concrete mixtures 

Gales oes le 1A 
considered, 

Se The relationship between asphaltic concrete stiffness, 
asphalt stiffness and volume concentration of aggregate 
proposed by Heukelom and Klomp (1964), Equation (VI-3) of 
CHAPTER VI, is shown graphically in FIGURE VIII-2. Super- 
imposed on this ais are the limits of the maximum 
predicted Cpoeaite induced stress-mix stiffness relation- 


ships shown in FIGURE VIII-1. 


FIGURE VIII-2 provides a means of readily estimating thermally 
induced stresses within asphaltic concrete pavements subjected to low 
temperature climatic environments such as experienced in Western Canada. 


Such estimations assume that the unknown thermal stresses can be 
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approximated by use of the pseudo-elastic beam analysis and that the 
coefficient of thermal expansion of paving mixtures equals i75ex 10m. 
per degree Fahrenheit and is temperature independent. With reference to 
Equation (VII-6), computed thermal stresses are directly proportional to 
a temperature independent expansion coefficient and therefore the posi- 
tion of the thermally induced stress limits may be shifted along the 
upper horizontal axis of FIGURE VIII-2 when a coefficient other than 


that assumed in the present study is used. 


At a given temperature, the stress estimate requires a knowledge 
of asphalt stiffness at a loading time of 7,200 seconds, together with 
the appropriate volume concentration of aggregate value (C,) for the 
paving mixture, or a knowledge of the asphaltic concrete stiffness at 
the same loading time. In the former case, the importance of defining 
the Cy value of the paving mixture is shown in FIGURE VIII-3. From 
this figure, relatively small variations in Cy can result in large 
differences between predicted thermally induced stresses in a paving 


mixture whose asphalt stiffness has been defined. 


For design purposes, the correlations obtained between pre- 
dicted and observed low temperature performance of the various pavements 
studied in CHAPTER VII suggest that comparisons between predicted 
thermally induced stresses obtained from FIGURE VIII-2 and measured 
tensile failure strengths of paving mixtures enable the low temperature 


fracture susceptibility of asphaltic concrete mixes to be readily 


assessed. 
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8.4. Comparisons Between Proposed and Previous ly Suggested 


Guidelines 

In order to compare the proposed method of assessing the low 
temperature fracture susceptibility of paving mixtures with previously 
Suggested methods a fracture strength-temperature relationship of an 
asphaltic concrete mix must be assumed. From the results of the tensile 
splitting test presented in CHAPTER VI and the results of direct tension 
tests on asphaltic concrete specimens reported by Haas (1968), the low 
temperature tensile fracture strength of asphaltic concrete paving 
. mixtures composed of asphalts obtained from Western Canada crude sources 
have generally been found to be within the range of from 200 to 400 psi. 
For the following comparisons, a paving mixture having a fracture 


strength of 300 psi at temperatures less than 0 F has been assumed. 


From FIGURE VIII-2 and considering the fracture strength of 300 
psi, low temperature cracking could be expected if the stiffness of the 
sphalerc concrete paving mixture exceeded 1.7 x 10® psi. This value 
corresponds favourably: with the stiffness limits of 1 x 10° to 2 x 10° 
psi Pereertcd'h McLeod (1969). If the volume concentration of the 
mix is 0.88 or 0.86 the maximum stiffness of the asphalt comprising 
the mix ne be limited to 12,000 and 22,000 psi, respectively, 

These asphalt stiffness values approximate the limiting value of 20,000 
psi at a loading time of 10,000 seconds suggested by Fromm and Phang 
(1971) and the limiting stiffness of 30,000 psi at a loading time of 


7,200 seconds adopted by the British Columbia Department of Highways. 
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While previously suggested limiting asphalt and asphaltic 
concrete mix stiffness values for minimizing low temperature transverse 
cracking are in close agreement with values derived from FIGURE VIII-2, 
this figure enables the influences of variations in asphalt and mix 


stiffness and volume concentration of aggregate on computed thermal 


stresses which have been correlated with observed field behavior to be 


evaluated, 


8.5. Summary 

In this chapter suggested design approaches for reducing 
transverse cracking of asphaltic concrete pavements have been briefly 
reviewed and a guideline based on the correlation obtained between 
predicted and observed times and temperatures of fracture of various 
pavements studied in CHAPTER III has been developed. Using this 
developed guideline the influence of asphalt mix stiffness, asphalt 
stiffness and volume concentration of aggregate on maximum predicted 


thermally induced stress can be readily deduced. 
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TABLE VIII-1 


MAXIMUM MIX STIFFNESS FOR SELECTING 


ASPHALT GRADE (McLeod, 1969) 


Minimum Stiffness Modulus i 
Temperature Cracking Cracking 
) Expected Eliminated 
1,000,000 500,000 


700 ,000 300 ,000 
400 ,000 200 ,000 


100 ,000 50,000 
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FIGURE WIIL-3 VOLUME CONCENTRATION OF AGGREGATE VERSUS 


MAXIMUM PREDICTED THERMALLY INDUCED STRESS 
FOR VARIOUS ASPHALT STIFFNESS MODULI. 
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CHAPTER IX 


SUMMARY, CONCLUSIONS AND RECOMMENDATIONS 


AS Ie ee Summary and Conclusions 


This study has considered one manifestation of non-traffic load 
associated distress of flexible pavements, namely, transverse cracking 
caused by thermally induced stresses exceeding the tensile strength of 
paving mixtures when subjected to low temperature climatic 


environments. 


Various factors found from previous studies to be related to the 
transverse cracking problem were reviewed. Principal findings of these 
Studies are: 

1) Transverse cracking of asphaltic concrete pavements can lead 
to significant performance losses and increased maintenance 
costs. In addition to being a source of maintenance ex- 
penditures, the reduction in pavement service life caused by 


transverse cracking may be of much more importance in an 
economic sense. 

2) Documented results of detailed field test projects, located 
in Western Canada and designed to study the low temperature 
response of asphaltic concrete pavements, indicate that the 


majority of transverse cracks initiate at a time when 
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a) 
pavement surface temperatures arc approaching or at 


seasonal minimum values. 


3) The low temperature fracture susceptibility of an asphaltic 
concrete mix is primarily dependent on the asphalt used. 
This dependency is reflected by differences in crack 
frequencies of pavements composed of asphalts obtained from 
different crude sources, incorporating asphalts of various 
grades as specified in terms of viscosity and/or pene- 
tration and including asphalts which exhibit different low 


temperature stiffness properties. 


4) The extent of cracking is not only dependent on the asphalt 
used but also on such variables as asphaltic concrete 
thickness, pavement age and subsurface materials. 
Generally, crack frequency is found to increase with de- 
creasing asphaltic concrete thickness, increasing pavement 
age and is often found to be highest within pavement 


structures having sand subgrades. 


The importance of temperature as related to the behavior of 
asphaltic concrete paving mixtures has long been recognized, however, 
it has only been in recent years that continuous temperature regimes 
within pavement structures subjected to low temperature climatic con- 
ditions and during periods of recorded transverse cracking have been 


collected. Analyses of such temperature data, obtained from field 


studies in Western Canada, has yielded specific information on time 


variations of temperature occurring in several Dav ete ee aes 
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subjected to climatic temperatures approaching -40 F. For the data 
analyzed, major findings are: 

1) At equal depths within asphaltic concrete surfaces of 
different thicknesses, temperatures and temperature 
variations are nearly identical when the pavements are 


subjected to similar climatic environments. 


2) Temperature gradients in asphaltic concrete surfacings 
increase with a decrease in daily minimum temperature, 
reflecting the existence of large stiffness and thermally 


induced stress gradients. 


3) Rates of surface temperature change are generally less 
than 3 degrees Fahrenheit per hour during winter months 
and at depths of 10 inches rates of asphaltic concrete 


temperature change are reduced to approximately 0.5 degrees 


Fahrenheit per hour. 


4) The thermal regime existing in a subgrade, subjected to 
prolonged subfreezing temperatures, is primarily dependent 


on the physical and thermal properties of the subgrade. 


In the absence of recorded pavement temperatures a realistic 


assessment of the response of asphaltic concrete pavements to low 


temperature climatic environments is dependent on a method whereby 


pavement temperatures can be reliably predicted. The application of 


principles of heat transfer in which the influence of climatic factors 


and pavement matcrial properties on pavement temperatures cnabled such 


a predictive method to be developed. Prior to its development, heat 
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transfer concepts applied to pavement structures and thermal properties 
of pavement material properties were reviewed. Principal findings of 
this review may be stated as follows: 

1) A rational approach to defining the influence of climatic 
variables on pavement surface temperature involves an 
energy balance procedure in which surface temperatures are 
related to various meteorological parameters. The accuracy 
of such a relationship is dependent on the method used to 
describe the influence of changing climatic conditions on 
surface temperatures, which will generally be governed by 


available meteorological data at a given locale. 


2) Documented results of studies related to thermal properties 
of pavement component materials enable approximations of 
the thermal properties of specific component materials to be 


made. 


Employing heat transfer concepts presented in the previously 
mentioned review, together with a finite difference approximation of the 
one-dimensional heat transfer equation, an analytical model for pre- 
dicting temperatures in pavement structures was developed. Formulation 
of the finite difference equations was presented and input parameters 
for the model were described. The parameters include structural, 
physical and thermal properties! of each component layer, meteorological 
data and a nonlinear relationship between unfrozen moisture content and 


temperature for each subsurface layer. Comparisons between predicted 


and recorded temperatures in various pavement structures were given. 
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The close agreement which was found to exist between predicted and 
recorded temperatures suggests: 
1) The model can be used to obtain reliable predictions of 
thermal regimes in asphaltic concrete surfacings under 


actual climatic conditions. 


2) The model enables thermal regimes existing in subgrades 
during periods of low temperature climatic environment to 


be approximated. 


Using recorded and predicted asphaltic concrete temperatures, 
thermally induced stresses within pavement surfaces of field test 
projects located in Western Canada were computed by various stress 
prediction methods. Such computations require definition of the time 
and temperature dependent stiffness characteristics of the asphaltic 
concretes and in order to assess the low temperature fracture suscepti- 
bility of the paving mixtures, assuming a stress fracture criterion, 
the fracture strength-temperature relationships of the mixtures must be 
defined. Concepts used to characterize the time and temperature 
dependence of the asphaltic concrete mixtures were described and re- 
sulting stiffness modulus-reduced time relationships were presented. 
The low temperature tensile properties of pavement cores obtained from 
the test projects were determined by means of the tensile splitting 
test. All tests were conducted using a loading rate of 0.056 inches per 
minute at test temperatures ranging from +20 to -20 F. The results of 
these tests were summarized and briefly discussed in relation to 


observed differences in low temperature behavior of test control 


Yo enol to fbeny atdsiier cletdo of been od oe mi 
vob egaisdtie a oltishges me coo ale 
‘atte t-nttedbte >isemet use 


a 


| a 
| paar yaa te, cinta iv 
2sbatgdv2 ai gnivaixs som Niger Lemyant oidnine ee ~ . 

; » we, 
oo Insmnorivne sissmil> swasteqmet wal 7 ao iraq gaieub be, 


amit 


( semmsdnidoqnet sio13009 =i tare 26 soso ang be vobavast poe x 7 


Gin we fi 
i ‘ a 
teot bie to: ane Tranteveag viata ‘eecconne beoubme 


re} 
key 


a 


zeavte 2uoitsy yd tine 9yaW sbenmnd avareel A 
omtt 33 YO noffiniteb erivper a foge niet Be 
ott fadqes ait to egiseiretosysi> seontbite mebnsgeb muds: 
~faquoeue Susoet? etiteredmes wol ot vendes os alibi 
git? @tutont2 eeott2 5 grnimrads .keuitxin ghtiveq: edt 20 ¥ 
en ont to eqidencitsior outs teynes~titgnstse wane <t 
, pruiterer deal omit of3 ov ire7onmBdo OF bode esqsonOd ~ 

| aotuTKin osoTaneD dM ages edd Yo ban 
pisefon omit becuase Ae 


259 


sections within the projects. Principal findings associated with this 


summary and discussion are: 


1) 


2) 


3) 


In terms of a relative comparison between different paving 
mixtures subjected to near identical low temperature 
histories, the paving mixture having the highest stiffness 
and lowest failure strain values exhibits the highest crack 


frequency. 


The low temperature tensile failure strengths of asphaltic 
concrete mixtures composed of asphalts obtained from 
Western Canadian crude sources range from approximately 
200 to 400 pounds per square inch and for a given mix an 
increase in density contributes to an increase in tensile 


failure strength. 


Differences in asphalt grade, as expressed in terms of 
penetration and/or viscosity, reflect differences in low 
temperature behavior. For a particular penetration at 77 F, 
an increase in viscosity at 140 F, can be expected to 
produce a decrease in crack frequency and for a particular 
viscosity at 140 F, an increase in penetration at 77 F, can 


be expected to produce a decrease in crack frequency. 


The stress prediction methods included in the study were 


pseudo-clastic and viscoelastic beam and slab analyses. Theoretical 


considerations for the formulation of the prediction methods were re- 


viewed and numerica 


Computed stres 


1 solutions of the stress equations were presented. 


ses were compared with fracture strength-temperature 
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relationships of the asphaltic concretes to predict times and tempera- 


tures of fracture. These predicted fracture conditions were then 


compared with times and temperatures of observed cracking of the 


corresponding pavements in the test projects. These comparisons per- 


mitted a realistic assessment of the various stress predictive methods 


from which the following conclusions are made: 


rE) 


2) 


3) 


a) 


Stress predictive methods coupled with the use of recorded 
or predicted temperature data and appropriate materials 
characterization provide a means of reliably assessing the 
low temperature fracture susceptibility of asphaltic 


concrete paving mixtures. 


Maximum and minimum thermal stresses are computed using 
the viscoelastic slab and beam analyses, respectively, while 
the approximate viscoelastic slab analysis yields inter- 


mediate stress values. 


Thermal stresses computed using the pseudo-elastic analyses 
are dependent upon the specified loading time used in the 
numerical integration of the stress equations. This 
dependency is primarily due to changes in asphaltic concrete 


stiffness with loading time. 


Using a constant time increment for numerical integration of 
the viscoelastic stress equations necessitates long computer 
runs and large storage requirements if stresses are to be 
predicted over extended periods of times. These require- 


ments can be reduced by employing a scheme of variable time 


increments. 
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5) The correlation between observed and predicted low tempera- 
ture performance of the asphaltic concretes suggests that, 
for the data currently available, the low temperature 
fracture susceptibility of an asphaltic concrete paving 
mixture can be reliably assessed using stresses computed 
using the pseudo-elastic beam analysis together with the 


assumed fracture criterion. 


In view of the correlation obtained between predicted and 
observed times and temperatures of fracture of the asphaltic concrete 
pavements included in the study, a simplified guideline for reducing 
the occurrence of transverse cracking caused by thermally induced 
stresses exceeding the low temperature tensile strengths of asphaltic 
concrete mixtures is proposed. Using this guideline the influence of 
variations in asphalt mix stiffness, asphalt stiffness and volume 
concentration of aggregate on thermal stresses computed by the pseudo- 


elastic beam analysis can be readily deduced. 


9.2. Recommendations 


Recommendations arising from this investigation and directed 
towards developing a solution to the low temperature cracking problem 
are listed as follows: 

1) Before any of the stress predictive methods and fracture 

criterion presented in this study may be applied with full 
confidence, further comparisons between observed and pre- 


dicted fracture times and temperatures of different pavement 


surfaces must be made. In the meantime the predictive 
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techniques should be implemented within the framework of 
present design, construction and material selection 
strategies so as to minimize the occurrence of low tempera- 


ture cracking. 


Future studies related to the low temperature cracking 
problem should be directed towards defining subgrade 


influences on crack initiation as well as crack frequency. 


Efforts should be directed towards developing effective 
pavement maintenance and rehabilitation strategies for the 
many miles of existing cracked pavements in areas subjected 


to low temperature climatic conditions. 
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APPENDIX A 


COMPUTER PROGRAM FOR PAVEMENT TEMPERATURE PREDICTIONS 


A.1l. Introduction 

This appendix contains details of a computer program developed 
for predicting temperatures in a one, two or three layered pavement 
Structure. Temperature calculations are based on finite difference 
equations which have been formulated using an explicit procedure and 
assuming a one-dimensional heat transfer problem. These equations are 


as follows: 


1) Temperatures within a component layer 


T(x,t+At) = a, [T(R-Ax sty FIT bet] 
2aAt 
‘ [ : att lot) (A-1) 


2) Temperatures at the interface of component layers 


T(x, ttat) = E2Ki47 Sty P(CxrAanye) 


i I 


Tet) 


2K; At 
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T(x-Ax,t) (A-2) 


3) Temperatures at an air-pavement surface boundary 


2aAt b> 2aAt T t 
T(x,ttAt) = er TGeax,t) (1 at (x,t) 
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(Quad . oA (A-3) 
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where 

T = temperature, 

au = depth, 

t = time, 

Ax and At = depth and time increments, respectively, 

O = diffusivity, 

K = thermal conductivity, 

Q ad = algebraic sum of short and long wave 


radiation heat fluxes, 


Q = heat flux resulting from convective 
heat transfer, 


i = refers to the layer 
I. = (Ax)? (Yi Ci + Yio Ci+i), where 
y and C are the dry density and 


mass heat capacity of the layers, 
respectively. 


The program consists of a main program and nine subroutines. A 
listing of the program, compatible with the operations of an IBM-360 
computer, is contained in this appendix. Within this listing many of 
the parameters used are defined and major program operations outlined. 
In the following section computations in the main program and each of 


the nine subroutines are described. 


ies Program Operations 


A generalized flow diagram for the entire program is shown in 
FIGURE A-1. From this figure, the operations of seven of the nine sub- 
routines intcract directly with those of the main program. Three of 


the subroutines, those titled Graphs, Graphl and Graph2, have becn 
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developed for Calcomp plotting of the calculated temperatures. When 
time-temperature plots are not required these three subroutines, to- 
gether with those statements in the main programereferring. to? Plot(s)e 


Xlimit, and Graphs, should be deleted. 


meal. Main Program 

The operations of the main program include specifying all input 
variables and the calling of subroutines for temperature calculations 
and plot displays. A flow diagram for the main program is shown in 
FIGURE A-2. With reference to this figure, operations within program 


involve the following sequence of events. 


1) Input variables are specified. 

These variables are defined in the program listing and 
are those used when calculating temperatures in two or three 
layered structures. For a one layered structure the number of 
layers (NL) must be specified as 2 and the fictitious uppermost 
layer, of an assumed thickness (TOL), is assigned identical 


thermal and physical properties as the layer under consideration. 


2) Subroutines Boundary, Nodal and Depth are called. 


Operations of these subroutines are given in subsections 


A'?.2 "At2.5 and AlZ.4, respectively. 


3) Initial temperature gradient is specified. 
A temperature TM(M,J) is assigned to each node point 
when time (J) is equal to l. The number of nodes (KOA) equals 


(total thickness of the structure/depth increment, Ax) + 1. 


sito bani anit ynistona dh 
bart?) Sabie hoviuped don ome ‘so 
(214084 9a) anc, igor nim odd ai ears 


.botolab od bluqde, 
3 ge" wah 


+. 


ell 2, a4 nes % or 


suqod Mie goivtioage: ohitfont mexyoty isa ot Yes <thoteacebap 
enoigalio tes rusted 03 sani pueda to ariliss. oat vee 


ry 


mt rewcaghe ai fetgorg niam oft tot mest st Ey wolt A sete 
TS TQOTH edt anoivsteqe .sroglt 2 idd oo ee Pe 


birsifliae 


+ | venave Re enreapie-yotwot fod 

| Aun 

host ropge ‘ome! getdakray aga (t 

Ape gaitelhl | me xgorq ‘ont pi Ronited oth asideitey Seon vl, 


AL ; aD 
at is owt TE gotaitersqmaet snisateniae wesw boat se0Hs 918 


rofl 9A stqrsbrte oe re el 
co bien ait bin 2 arariagi mo pada 

{notinobt Romper’ ef . (OT) evoiietay bemmeas & eet 

-aabie) gous) en wwsishie lead 
oe oy 4): ee ee i: 
169. 65 oy a son — IsboK .Yvtsbaod cosbrssalt : te ; ee 


. dy ae it 206 ow bee ee 
Peni tuordy2 oeadis to : 


pi «Yiov isosqeot ~b. S.A sme cah wet 


pone eel 
Vertes 
“7 eal 


A4 


The temperature of the deepest node (M = KOA) is assumed 


constant throughout the time period of temperature calculations. 


4) Counters LLA, LLB and LLP are defined. 

These counters, which are used throughout the entire 
program, direct program operations within subroutine Output 
where various commands for displaying the calculated tempera- 


tures are specified. 


5) Climatic inputs are specified. 

The date, air temperature and solar radiation values 
are read. The time of initial specified temperature and 
radiation values must correspond to the time at which the 
initial temperature gradient is defined. Solar radiation 
values are expressed in units of langleys per hour. However, 
radiation values in units such as Btu per hour may be specified 
provided that modifications are made within subroutine tempera- 


ture as outlined in subscction A.2.7. : 


6) Subroutine Airsun, Temperature and Graphs are called. 


Operations within these subroutines are described in 


Subsections A.2.6, A.2./7 and 2220, respectively. 


7) Temperature gradient is defined. 

The last computed temperature gradient within a 
structure and for a given day's climatic input become equal to 
the first temperature gradient for the following day. 
Operations then return to sct 5 where the date, air temperature 


and radiation values for the following day are read. 
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A.2.2. Subroutine Boundary 
In subroutine Boundary multiplication constants of finite 
difference Equations (A-1), (A-2) and (A-3) are calculated. An outline 


of operations within this subroutine is shown in FIGURE A-3. 


1) Diffusivity of each layer is calculated. 

Using both unfrozen and frozen thermal properties the 
diffusivity of each layer is computed. These values are then 
multiplied by the mesh ratio (At/Ax2) to obtain constants Re 


and Re for each component layer. 


2) Latent heat released upon freezing of the lower layers 
US calculated: 

Latent heat released upon freezing of the uppermost 
layer is assumed equal to zero. For each of the lower layers, 
latent heat released upon freezing (Btu per cubic foot) is 
cvaluated by determining the product of, (1.44) (moisture 
content) (dry density) and the per cent of water frozen (PERWF) 
at temperature (F) which defines the limit of the "while" 


freezing condition. 


3) Constants of the finite difference equations are computed. 
Thermal conductivity, heat capacity, dry density values 
and time and depth increments specified in the main program are 


used to compute multiplication constants of the finite 


difference equations. 
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4) Temperature-apparent specific heat relationships of the 
lower layers are calculated. 
Apparent specific heat values at the temperature 
limits defining the "while" freezing condition, temperatures B 
and F, are denoted as CWF and CF, respectively. During 
freezing, a linear relationship is assumed between the logarithm 


of these apparent specific heat values and temperatures B and F. 


A.2.3. Subroutine Nodal 

In this subroutine nodes at which temperatures are to be dis- 
played in the output are determined. Temperatures at all nodes in the 
uppermost layer are displayed. For example, if the depth increment 
(DELX) equals 2 inches and the thickness of uppermost layer equals 6 
inches then calculated temperatures at depths of 0,2,4 and 6 inches 
are shown in the output. The depth increment at which temperatures 
in the lower layers are displayed is controlled by the factor KMZ 


defined in the program listing of this subroutine. 


A.2.4. Subroutine Depth 


Subroutine depth calculates the depths of the nodes determined 
in subroutine nodal. These depths are shown in the output in the form 
of headings under which the calculated temperatures are printed. 

Using subroutines Graphs, Graphl and Graph2 plots of these temperatures 


as a function of time can be obtaincd. 
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soe Subroutine Output 


The operations in this subroutine, in which title headings, 
depths and temperatures are written, are described in the program 
lasting. = In addition to this output, commands are given for the last 
computed temperature gradient to be punched on cards. This temperature 
gradient, together with appropriate climatic variables, can be used as 
input to the main program when temperatures during the succeeding time 


period are to be calculated. 


A.2.6. Subroutine Airsun 

In this subroutine air temperature (TA) and solar radiation 
values (HS) are calculated at time increments of At (DELT). For these 
computations both air temperature and radiation values are assumed to 
vary linearly between the times in which they are specified in the 


main program. 


A.2./. Subroutine Temperature 
This subroutine calculates temperatures in a pavement structure 
at time and depth increments of DELT and DELX, respectively. Program 


operations of this subroutine arc outlined in the flow diagram shown in 
FIGURE A-4. 
1) Surface temperature. 
At time J+1 and depth M, M=1, the surface temperature 
TM is calculated by means of Equation (A-1). Hleat transfer by 
convection (CONV) is approximated by the product of convection 


coefficient, defined in the main program, and the difference 
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A8 
between air and pavement surface temperature at time (J). Heat 
transfer due to longwave radiation (RAD) is approximated by the 
product of emissivity (EM), Stefan-Boltzmann constant (SD) and 
difference between the fourth powcr of air and surface tempcra- 
ture. Shortwave radiation heat transfer is computed by 
multiplying the absorptivity of the surface (A) by the intensity 
of solar radiation (Il) at time J. Since solar radiation values 
are specified in the main program in units of langleys per hour, 
a constant of 3.69 is used to obtain equivalent values in units 


Of Be, perahour. 


2) Temperature at depths. 
At time J+l the temperature at each node in the 
structure is determined. This computation is dependent on one 


of the three following conditions. 
a) whether the node is located at the interface of 
adjacent layers which are in an unfrozen or frozen 
States; 
b) whether the node is located within a layer which 


is in an unfrozen or frozen state, or 


c) whether the node is located within a layer(s) 


which is in a "whilc' freezing condition. 
kor each of these conditions the location of the node is dc- 


fined by the magnitude of M with respect to the number of 


nodes in each layer. Unfrozcn, frozen and freezing conditions 
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A9 
are defined by the magnitude of temperatures TM(M-1,J), TM(M,J) 
and TM(M+1,J) with respect to the specified temperatures B and 


F defining the freezing zone. 


3) Time interval for temperature display is specified. 

In the program listing a time interval of 2 hours has 
been selected for temperature display. This time interval is 
defined by the ‘summing factor SDEL which increases in magnitude 


by a factor of DELT with every increase in counter J. 


4) Call Output. 

If the factor SDEL, described in step 3, becomes equal 
to or greater than the seclected time interval of temperature 
output subroutine Output is called. If this condition is not 
satisfied operations return to step 1, where temperature 


calculations for the following time step commence. 


Subroutine Graphs 


In this subroutine temperatures to be plotted are converted 


from array to vector notation. These temperatures are equivalent to 


those displayed in the output. ‘This subroutine also calls for plot 


subroutines Graphl and Graph2. 


are plotted. 


Subroutine Graphl 


Subroutine Graphl establishes a grid on which the temperatures 
In the subroutine the scales selected for this grid are: 
3 inches horz. = 24 hours 


2? inches vert. = 10 degrees F. 


by, 4 


(Leer itl, ee ieee een as 


bre @ aaa 


Ta | 
at 


an ne. > ’ } te ee “eyes \¥ : * 


een Sid caeai Sry seroqnes tot ck iy | 
aid erpod S: te oso - ® gnrsert maTgZorg oie fog 
73 ii OM ° iv q’ 


ait Levey anid) aia Uatek i stusterxvoqmes 108 bas ah (ees Maes. 


obs Lianne eshorant | (foil d402. 10t903 _ gevhaanite ott eh tn 
iG Lads anh jeeRonont cove age iw Taad eo pees) 


sie 8 a 


fsups ele Bs (oz i petitiozab ol AGE voanet eat a o4 ia 
eae | 


srutstsqi’® th diesel aukt hesooloe ody wats batgehie 10 7 
ee 
Ton ek i doisibron er at bolls ef tuq3yO aaa suqgue eo! 
oT iserbquey eAevite 2 qate ot Panter ae: 
| t q amet ssitodo 7; iat aor 4 ive Y 


_bontgnntod gate oars Sele pants ano tts te tgs” , 
if J 
;% a th pera | (ait . a 


nn ais bpszolg od of cru seroqHed On ae ets 


; Ss20adt +7 
os RS mh eomuteT MOT O20 -noi 280m, rosoov" ° 7 i 
toly aise Pan nied Pale: sit diveredue ein? FO INO: ads nk ~ 


-Silqjnx9 ‘ ges 


i Be 


i wiedinta 
as , 
nia 


a) on wo wabatos eslns2 ody on ? I xd 


ae 


A10 
A description of the structure is written in the upper lefthand corner 


of the grid and the ordinate is labcled. 


A.2.10. Subroutine Graph2 


In subroutine Graph2 calculated temperatures displayed in the 
output are plotted and for a given depth are shown as a continuous 
function of time. Jlaving plotted one day's temperatures operations 


return to main program and step 5 described in subsection A.2.1. 
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THIS PROGRAM UTILIZES A FINITE OILFFERENCE TECHNI QUE TO_ 
CALCULATE THE THERMAL REGIME IN A ONE»TWO OR THREE LAYERED 
PAVEMENT STRUCTURE SUBJECTED TO VARIGUS CLIMATIC ENVIRCNMENTS. 
A ONE OIMENSIGNAL HEAT CONOUCTION PROBLEM IS ASSUMED. 


THE PROGRAM CONSISTS OF A MAIN PROGRAM ANO NINE SUBROUTINES e 
SIX OF THE NINE SUBROUTINES INTERACT TO CALCULATE AND WRITE OUT 
THE GENERATED TEMPERATURES.THE FUNCTION OF THE THREE REMAINING 
SUBROUTINES »sTHOSE ENTITLED GRAPHS>»GRAPHL AND GRAPH2sI1S TO PLOT 
THE GENERATED TEMPERATURESeWHEN PLOTS ARE NOT REQUIRED THESE THREE 
SUBROUTINES COUPLED WITH THOSE STATEMENTS IN THE MAIN PROGRAM 
REFERRING TO PLOT(S)»sXLIMIT AND GRAPHS SHOULD BE REMOVED. 


MAIN PROGRAM 
We Me te Gm he ok ek 
COMMON KU (3) sKF (3) «CU(3)sCF(3).DDEN( 3) o> wc( 3) » TOL( 3) e+ RUC 3)sRF( 3) 
1CWF (4) eWOEN( 3) »eNOP(3) »HS(29) eH(200)2TA(30)2TM(1152200) »sPERWF( 2), 
2TAT(200) »TITLEI (20) »TITLE2(20).IDEEP( 30) »,SLOPE(3) ».81(3) .G6(30). 
3BUF (2048) » DELXeKOoKMZs NL »NOBC eNPBeLNPBeNPSoLNPSeKOAs TOTAL» 
4LAST ol_KL »>KOLKL oKLK sKOLKL1I 2® KO] o CAs CAF o CFF » CGF eCB»CBF se CC+CCFs CDsCDF >» 
5CE »CEF »«CHF eCIFeBeF eDELT sJJJolLLitsRReIILstLPsArsEMsUCeNOTI eo Xe 
6TYME( 30) sMeNDeLLBsHRollLAeJeoINC - 
REAL KUeKF 
INTEGER HReX 
CALL PLOTS (BUF 38192) 
CALL XLIMIT( 250-0) 


THE FOLLOWING PARAMETERS ARE READ INS 


ND —-— THE NUMBER OF DAYS FOR WHICH THE THERMAL REGIME IS 
TO BE CALCULATED. 

NU — THE NUMBER OF LAYERS IN THE PAVEMENT STRUCTURE.s 

THE NODAL POINT NUMBER THAT DEFINES THE UPPERMOST DEPTH 

IN THE LOWER LAYERS AT WHICH TEMPERATURES ARE TO 

COMMENCE DISPLAY IN THE OUTPUT. 

THE NUMBER CF TEMPERATURES TO BE PLOTTED PER DAY 

AT A GIVEN DEPTH IN THE STRUCTURE (IF X=13 TEMPERATURES 

ARE PLOTTED ON A BIHOURLY BASIS»MIONIGHT TO MIDNIGHT. 

IF X=25 TEMPERATURES ARE PLOTTED ON AN HOURLY BASIS.) 

THE NUMBER OF AIR TEMPERATURES SPECIFIED PER DAY. 

(IF JJJ=25 AIR TEMPERATURES ARE READ IN HOURLYs 

MIONIGHT TO MIONIGHT o) 

THE NUMBER OF SOLAR RADIATION READINGS SPECIFIED PER 

DAYe(SIMILAR TO X AND JJJ ABOVE.) 

THE TIME INCREMENT IN HOURS AT WHICH TFE TEMPERATURES 


ARE TO BE CALCULATED. 
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THE DEPTH INCREMENT IN INCHES AT WHICH THE TEMPERATURES 
ARE TO BE CALCULATED. 

THE TEMPERATURE AT WHICH FREEZING COMMENCES IN THE 
SUBGRADE. 


THE TEMPERATURE AT WHICH A GIVEN PERCENTAGE OF WATER 
IN THE SUBGRADE IS FROZEN. 


REAO(Se21)INOoNL eo INCs Xo JIJsLLLL oe DEL Te DEL Xo Bo F 
1 FORMAT(615+4F1003) 


THE SURFACE 


PS 
EM = 
uc - 


READ(5+2)A0EMsUC 
2 FORMAT(3F10-3) 


COEFFICIENTS ARE READ IN3 


ABSORPTIVITY 
EMISSIVITY 
CONVECTIGN COEFFICIENT (BTU/HR/DEGREE F/SQO FT) 


THE FOLLOWING PARAMETERS ARE READ IN FOR EACH LAYER GF THE STRUCTURE; 


KU — 


KF - 


CUD a= 


cF - 


VAS 9 


ToL - 


DO 4 JT=1 NL 


THE UNFROZEN THERMAL CONOUCTIVITY OF THE LAYER. 
(BITU/HR/FT/OEGREE F) 

THE FROZEN THERMAL CONDUCTIVITY OF THE LAYERe 
(BTU/HR/FT/SOEGREE F) 

THE HEAT CAPACITY OF THE LAYER IN AN UNFROZEN 
STATE se {(BTU/SLB/DEGREE F) 

THE HEAT CAPACITY OF THE LAYER IN A FROZEN 
STATE. (8TU/LBG/OEGREE F) 

THE ORY DENSITY OF THE LAYER-e(LB/CUBIC FT) 
THE MOISTURE CONTENT OF THE LAYERe 

(PERCENT OF ORY UINT WEIGHT) 

THE THICKNESS OF THE LAYERe{ INCHES) 


READ(5e3)KU(C JT) e0KE CIT) »CUC JT) »CF( UT) eODEN( JT) »WC( JT) »TOL( JT) 


3 FORMAT(7F10e3) 
4 CONTINUE 


THE PERCENT OF WATER FROZEN AT TEMPERATURE F IS READ IN FOR 
EACH LAYER WITH THE EXCEPTION OF LAYER ONE» THE PAVEMENT.CIF PERWF=0.0 
TEMPERATURES ARE CALCULATED ON THE BASIS OF THE UNFROZEN 
PROPERTIES OF THE LAYER.IF PERWF=100-0 TEMPERATURES ARE 
CALCULATED ON THE BASIS OF THE FROZEN PROPERTIES OF THE LAYERSIF 
PERWF>0e0 OR < 100.0 LATENT HEAT IS CONSIDERED IN A MOVING ZONE OF 
FREEZING BETWEEN THE SPECIFIED TEMPERATURES 8 AND Fe) 


DO 6 NULF=2eNL 


READ( 5S) PERWF(NLF ) 
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FORMAT (1F10e3) 
CONTINUE 


TITLE! APPEARS ON EACH PAGE OF OUTPUT ANO AS A TITLE ON THE PLOTS. 


READ(Se7)(TITLEICI )»~1=1,20) 
FORMAT (20A4) 

CALL BRODY 

CALL NOOAL 

CALL DEPTH 


THE INITIAL TEMPERATURE GRADIENT IS SPECIFIED.THE GRADIENT 
CONSISTS OF A TEMPERATURE FOR EACH NODAL POINT. THE TEMPERATURE OF THE 
OEEPEST NODE IS ASSUMED CONSTANT THROUGHTOUT THE TIME PERIOD IN 
WHICH THE TEMPERATURES ARE BEING CALCULATED. 


J=1 

READ(5.8)(TM(Ms J) »M=1 oKOA) 
FORMAT (SFLOe1) 

LLA=0 

LLB=0 

LLP=1 

OO 150 III=ieNO 

HR=0 


DAILY DATA IS READ INe 


TITLE2 — THE DATE 
TA —- AIR TEMPERATURES(DEGREES FAHRENHEIT) 
HS - SOLAR RADIATION VALUES(LANGLEYS PER HOUR) 


READ(5e9)(TITLE2( 1) ol=1 020) 
FORMAT (20 A4) 

CALL OUTPUT 
READ(Se2 10)( TAC J) » J=he JIS) 
FORMAT (1L3F601) 
READ(5e11)CHS( J) o J=Ro LLLL) 
FORMAT (13F601) 

CALL AIRSUN 

CALL TEMP 

CALL GRAPHS 

J=1 

DO 12 M=1+eKOA 
TM(Me J) =TM(MeNOT I) 

CONTI NUE 

CONTINUE 

LLB=1 

CALL OUTPUT 

CALL PLOT (0.02000:999) 
STOP 

END 
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SUBROUTINE BRDY 


SUBROUTINE BOUNDARY 
EERE RE KEK EERE KEKE KE 


SUBROUTINE BOUNDARY CALCULATES THE HEAT CONDUCTION FARAMETERS 
OF THE ENERGY BALANCE EQUATIONS AT THE ALR-SUREACE AND ADJACENT 
LAYER INTERFACES. 


COMMON KUC3) sKF (3) »CU(3) »CF (3) sDDEN(3) »WC(3) sTOL(3) »RU(3)sRF(3)o0 
I1CWF (4) »WOEN(3) »NOP(3) 5HS( 29) »H( 200) 6 TA( 30) 3 TM(115%200) »PERWF( 2)» 
2TAT(200) »>TITLEL (20) eTITLE2 (20) » IDEEP( 30) »SLUPE(3).B81(3)sG6(30)> 
IJBUF (2048) »sDELX »s KOe KMZ o NL »NOBC »NPBeLNPBsNPS pL NPS sKOAs TOTALS 
4LAST »LKLeKOLKL »KLK sKOLKL1+KO1 » CAs CAF s CFF » CGF o CB eCBF eCC CCF sCOsCOF, 
SCE oCEF o CHF »CI FeBeF eDELTsJJJoLLLL ea RReIII ol lL PeAsEMsUCesNOTIoXo 
6B TYME (30) eMsNDeLULBeHRetlAv J» INC 

REAL KUsKF»NOP 

TOTAL=0-0 


THE OIFFUSIVITY OF EACH LAYER IS CALCULATED. 


DO 2 JT=1>,NL 

RUCIT) =CKUC IT) 7 (CU JT) X®DOENC JT) )) FOELT/S( (CDELX/ 126 )**2) 
RFC IT )=CKFC IT) /CCFE (IT) ®DOEN(C IT) D) ¥OELTZ( (DELX/126) *¥2) 
TOTAL=TOL( JT ) +TOTAL 

CONTI NUE 


THE LATENT HEAT IS CALCULATED IN ACCORDANCE WITH THE PERCENTAGE 
OF WATER FROZEN AT TEMPERATURE F SPECIFIED IN THE MAIN PROGRAMe 


DO 4 NLF=2».NL 

WOEN(NLF ) =ODEN( NLF )¥ (10 tC WC(NLF)/100¢) ) 

CWF (NLF )=( 1244*%WC( NLF) *(PERWF(NLF) 7/1002) *DDEN(NLF ))/WOEN(NLF) 
CONTINUE 

NOP(1)=TOL( 1) /DELX#1 


KO=NOP(1) : 
RR= (26 *(KUC1)7(0CU(1 ) *¥ODEN(1))) *DELTI/S(KUCL)*(DELX/126.)) 


TNP=TOTAL/SDELX+t+1le 
KOA=TNP 
THE INTERFACE CONSTANTS OF THE FINITE OILFFERENCE EQUATIONS 
ARE DETERMINED. 


IF (NL-2)1008s6 
CA=((DELX/ 126 )**2) #(DDEN (2) *CU( 2) #DDEN( 3) *CU03)) 


CAF=( (DELX/12 6) #¥*2) #( DDEN( 2) *CF (2) +D00EN( 3) #CF( 3D) 
CFF=( (DEL X/12 2) ##2) #( ODEN (2) * CF (2) FODEN( 3) *CU03)) 
CGF=( (DELX/12 6 )**2)*(ODEN( 2) *CU(2) #ODEN ( 3) *CF( 3)) 


CB=2-*KU(3)*D0ELT 
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CBF=2- *KF(3)*D0ELT 
8 CC=2-*KU(1) *OELT 
CCF=26*KF (1) *DELT 
CO=2 6 *¥KU( 2) ¥DELT 
COF=2 6-* KF (2)*DELT 
CE=( (DELX/1L20 )**2)*(DDEN(1 )* CUCL) #DDEN(2) *CU2) ) 
CEF=( (DOELX/12 6 )*¥*2) *( ODEN (1) #CF (1) tOD0EN( 2) *CF( 2) ) 
CHF=( (DEL X/12 « )**2) *(ODEN( 1) *€CF (1) #DD0EN( 2) *CU(2)) 
CIF=((DELX/124)**2)*€(DDEN(1)*CU( 1) (1DDEN(2)*CF(2)) 
10 CONTINUE 


THE RELATIONSHIP BETWEEN APPARENT SPECIFIC HEAT AND TEMPERATURE 
IS DETERMINEO FOR EACH OF THE LOWER LAYERSe 


OO 12 I0=2.NL 
SLOPE (10) =(0-434*ALOG( CWF (IO) )—00434*ALOG( CF(ID)))7006.434*ALOG(B) 
1-0-434*ALOG(F) ) 
BI(10,)=0 2434¥*ALOG(CWF( ID) )-SLOPE(ID) *06434*ALOG(8) 
12 CONTINUE 
RETURN 
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SUBROUTINE NODAL CALCULATES THE NODAL POINTS FOR WHICH THE 
GENERATED TEMPERATURES ARE TO BE DISPLAYED IN THE OUTPUT. 


COMMON KU(3) eKF (3) eo CUC3)sCF(3)sDDEN(3) »WC(3) eTOL{ 3)» RUC 3) sRFC 3)» 

1CWF (4) »WDEN(3) eNOP(3) sHS(29).H( 200)» TAL 30) 5 TM( 1150200) sPERWF(2) » 
2TAT(200) sTITLEL (20) oe TITLE2(20).IDEEP( 30) »SLOPE(3).81(3)+G6(30), 

IBUF (2048) »DELX » KO» KMZ 4 NL »>NOBC »NPSeLNPBeNPS oLNPSdXOAc TOTAL 2 
4LAST >sLKL»KOLKL eKLK »KOLKLI1 »KO1 » CAs CAF » CFF s CGF »CB oe CGF » CC »CCF 2 CD eCOF 2 

SCE » CEF » CHF oCIF eBeF sDELT sJJJoilll»RReITITLeLL PsAsEMsUCeoNOTI9 Xo 

6TYME ( 30) oMesND oL LB eHReolLLAeoJoI NC 

REAL NOP 

TEMPERATURES AT ALL NOOAL POINTS IN LAYER ONE ARE OISPLAYED 
IN THE OUTPUT. 


KO — THE NUMBER OF NODES IN LAYER ONEe 
KMZ — THE NOOAL INCREMENT IN THE LOWER LAYERS AT WHICH 
TEMPERATURES ARE TO BE WRITTEN CIF KMZ=12 AND OELX=2-0 
INCHES sTEMPERATURES ARE WRITTEN AT 24 INCH INTERVALS.-) 
NOB ~— THE NUMBER OF NODES IN LAYER TWOe6 
NOBC — THE FIRST NODE IN LAYER TWO AT WHICH TEMPERATURES ARE 


WRITTENe 
LNPB - THE LAST NODE IN LAYER TWO AT WHICH TEMPERATURES ARE 
WRITTEN? 
NPS — THE FIRST NOOE IN LAYER THREE AT WHICH TEMPERATURES ARE 
WRITTENe 


SUM=—-DELX 

DO 1 NP=1.K0 
SUM=DELX+SUM 

CONTINUE 

KMZ=12 

NOP (2) =TOL(2)/D0ELX 
NOB=NOP (2) 

NOBC=KO+FINC 
NPB=NOBtKO 

ABZ=KMZ 

AB=ABZ *¥DELX 

IF (NPB e«LT.- NOBC) GO TO S 
STOT=TOL(1)+#TOL(2) 

DO 3 NB=NOBC +NPBeKMZ 
AC=STOT—-SUM 

IFCAC eLTe AB) GO TO 7 
SUM=SUM+AB 
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LNPB=NB 

NRE=NP6B-NB 

IF CNRE el Te KMZ) GO TO 7 
CONTINUE 

NOBC=KO+1 

LNPB=KO+l 

SUM=SUM4+DELX 

CONTINUE 

IFCNL eLTe 3) GO TO 11 
NPS=LNPBt6 

DO 9 NS=NPS eKOAcKMZ 
AS=TOTAL-SUM 

IFCAS eLT. AB) GO TO 11 
SUM=SUM+AB 

LNPS=NS 

NER=KOA-—NS 

IF (NER el Te KMZ) GO TO i2 
CONTINUE 

RETURN 

ENO 
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SUBROUTINE DEPTH 


SUBROUTINE DEPTH 
MEK ER EEE 


SUBROUTINE OEPTH CALCULATES THE DEPTHS WITHIN THE PAVEMENT 
STRUCTURE OF THOSE NODAL POINTS DETERMINED IN SUBROUTINE NODALe 


COMMON KU(C 3) 2KF (3) -CU(3) »-CF(3) 2.0DEN(3) »wC( 3) »TOL( 3) ~-RUC3) RFC 3). 
LCWF (4) »eWOEN(3) sNOPC( 3) otH1S(29)+H(200)5TA( 30) 5TM( 1152200) sPERWFL 2) » 
2TAT(200) sTITLEL (20) eTITLE2 (20) 2 IDEEP( 30) eSLOPE(3).B1(3)»,G(30)>5 
3BUF (2048) e-DELX » KO eKMZ».NL »NOBC »NPBoLNPBsNPS eLNPSsKOAsTOTALs 
4L AST »sLKL eKOLKL oKLK »KOLKLI »KO1l » CAs CAF sCFFeCGF»C8» CBF» CC»CCFeCDsCOF » 
SCE » CEF » CHF oCIF eBeF eDELT »pJJJolLLLisRReIII sil PsAsEMsUCs NOTIsX» 
6TYME( 30) eMeNDolLLBoHRetlAeJso INC 

I=1 

JUD=DELX 

DO 2 NP=1+sKO 

IDEEP( I )=JUD*NP-—JUD 

I=I+1 

CONTINUE 

LAST=KO 

LKL=0 

DO 4 NB=NOBCeLNPB» KMZ 

IDEEP( I )=JUO*NB—-JUD 

=[+1l 

LKL=LKL +41 

CONTINUE 

KOI=KO+rl 

KOLKL=KO+4#LKL 

LAST=KOLKL 

IF(NL el Te 3) GO TO 8 

KLK=0 

O00 6 NS=NPSeLNPSeKMZ 

IDEEP (I )=JUO*NS—JUD 

I=I+l 

KULK=KLK¢41 

CONTINUE 

KOLKLI=KOLKLt1 

LAST=KO4#LKL+KLK 

RETURN 

END 
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SUBROUTINE OUTPUT 


SUBROUTINE OQUTPUT 
BERK EK EEK KEKE KEE RK 


SUBROUTINE OUTPUT WRITES OUT THE DEPTHS CALCULATED IN 
SUBROUTINE DEPTH ANO THE TEMPERATURES AT THESE DEPTHS e THIS 
SUBROUTINE ALSO CALLS FOR THE PUNCHING OF THE LAST TEMPERATURE 
GRADIENT WITHIN THE STRUCTUREe THIS TEMPERATURE GRADIENT BECOMES THE 
INITIAL TEMPERATURE GRADIENT TO BE SPECIFIED IN THE MAIN PROGRAM 
FOR THE SUCCEEDING TIME PERIOD IN WHICH THE THERMAL REGIME IS TO 
BE CALCULATED. 


COMMON KU(3) KF (3) » CU(3)eCF(3)».DDEN(3) oWC( 3) eTOL(3)sRU(3)»RE( 3) 
1CWF (4) »sWOEN(3) »NUP( 3) »HS(29) sH( 200) »TA( 30) 5 TM( 115-200) »sPERWFC( 2) » 
2TAT(200) sTITLEL (20) eTITLE2( 20) » IDEEP( 30) »SLOPE(3).B1I(3)2G(30)>, 
3BUF (2048) »sDELX » KO sKMZ NL »NOBC »NPBsLNPB eo NPS oL NPS oKOAs TOTAL» 
4LAST eL KL »« KOLKL »KLK oKOLKL1+KO1 »CAsCAFs CFF sCGF »CB»CBF » CC 2CCF » CD sCOF » 
SCE «CEF e CHF eCIF 8 sF .OELT »p JJJeLittl» RReIfI sil P»eAsEMeUCsNOT [eX 
6TYME( 30) +MeND oLLB oHR pLLAsJsINC 

INTEGER HR 

IF(LLA eEQe 1) GO TO 27 

IF(LLB eEQe 1) GO TO 35 

IFCLTI-LLP) 23.123 

LLP=LLPt+3 


HEADINGS FOR THE OUTPUT ARE WRITTEN. 


WRITE (C6e3)(TITLEL(C I) sf=1.20) 

FORMAT CIH1L 0 /25X p20A4S/ ) 

WRITE(695) 

FORMAT (50Xe15HOFPTH IN INCHES//) 

IFCNL-2) 11leolle7 

WRITE(6 09) 

FORMAT (2X5 3HHRe eo 1 SXeTHLAYER 1223Xe7HLAYER 20625Xe7HLAVER 3//) 
GO TO 15 


WRITE (6213) 
FORMAT (2X9 3HHRe cL BXeTHLAYER 1235Xs7HLAYER 2/7) 


THE DEPTHS WITHIN THE STRUCTURE ARE WRITTEN’ 


15 WRITE(6017)( IDEEP(I )oI=1eKkO) 


2) 


FORMAT (SXe71S) 

WRITE (6019) (IOEEP( 1)» 1=KO1 »KOLKL) 
FORMAT ( 1LH+e40Xs9I5) 

IF(NL eLT. 3) GO TO 23 

WRITE (6021) (IDEEP( I)» T=KOLKL1 »LAST) 


FORMAT (1H*»60Xe1015S) 
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THE DATE IS WRITTENe 


WRITE (6025) (TITLE2(1) »1=1 520) 
FORMAT (/420X»20A4/) 
J=0 


THE GENERATED TEMPERATURES ARE WRITTEN AT A TIME INTERVAL 
SELECTED IN SUBROUTINE TEMPERATURE « 


WRITEC 6229) HRoe( TM( Meo JS+1) »M=159KC) 
FORMAT (3X9 T2s7F501l) 

WRITE C6031) ¢ TMC Ms J41) o> M=NOBCeoLNPBs KMZ) 
FORMAT (1H4+040Xo9F Sel } 

IF CNL eL Te 3) GO TO 43 

WRITE C6233) (TM( Ms J+] ) » M=NPSeLNPS sKMZ) 
FORMAT (1H¥260X,10F5e1) 

GO TO 43 


° THE LAST TEMPERATURE GRAOIENT IN THE STRUCTURE IS WRITTEN AND 
PUNCHED ON CARDS. 


WRITE (6537) 

FORMAT(1H1 o5X%es25HLAST TEMPERATURE GRADIENT ) 
00 39 M=1eK0As8 

L=M+7 

IF(tL e6GTe KOA)L=KOA 

WRITE( 6041 )0 TM( Nol) sN=Mol) 
WRITE(7 041 )(TM( Nol) oN=MoL) 
CONTI NUE 
FORMAT (80 1X sF90e2)) 

RETURN 

END 
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SUBROUTINE TEMP 


SUBROUTINE TEMPERATURE 
RK TORR Rk tok kk a mR tk dk & 


SUBROUTINE TEMPERATURE CALCULATES THE TEMPERATURES WITHIN THE 
PAVEMENT STRUCTURE AT EACH NODAL POINT AND AT EACH TIME INCREMENT 
AS OEFINED BY DELX AND DELT RESPECTIVELY. 


COMMON KU(3).KF (3) eCU(3) »CF( 3) s0ODEN(3) eWC(3)sTOL(3) »RU(3) »RF (3). 
LCWF (4) » WOEN( 3) eNOP( 3) 9HS(29)0t(200)s TA( 30) o TM( 1159200) »sPERWF(2) » 
2TAT( 200) »TITLE1 (20) »sTITLE2(20) >» IDEEP( 30) »SLOPE ( 3) G1 (3) eG(30)> 
3BUF (2048) eDELX » KO eKMZ oe NL »NOBC oe NPB oe LNPB eNPS oL NPS oKUAs TOTAL» 
4ALASTslL KL» KOLKL oKLK eKOL KL! »KOL » CAs CAF» CFF» CGF »CB»CBFsCCeCCF»sCOsCOF » 
SCE 9 CEF » CHF eo CIF 0B oF eDELToJJJolLiLL es RReoIII vLLP»sAsEMs UCe NOTI Xo 
6TYME(30) pMeoNDeLLBoHReLLAs Js INC 

REAL KU.»KF 

INTEGER HR 

OL=24 -/DELT 

NOT=0L 

NOTI=NOT+H¥1 

JFC=KOA-1 

SDEL=0.-0 

SO=.1714E-908 


THE TIME OURING THE DAY AT WHICH THE TEMPERATURES ARE BEING 
CALCULATED IS DEFINED BY Je 


DO 67 J=1-.NOT 

M=1 

CONV=UC#*¥(TAT(J)—TM(MsJ)) 

RAD=EM*SD# (COC TAT( J) 4460 6) ©#4)— CC TM(Ms J) +460 6 )**4) )D 

BC= (3 -69*A#H( J) 2) +CONVtRAD 

TM(Ms J#1) =(1 0-2 eo FRUCL) PFTMO Me JD #02 0*RUC 1 ID¥TM(M4+10659)4+(RR*¥BC) 


TM(KOAs J+1)=TM(KOA2 J) 


THE TEMPERATURE AT EACH NOOAL POINT IN THE STRUCTURE IS 
DETERMINED AT TIME Je 


DO 61 M=2,JFC 


IF (M—KO) 290193 
IF CTM(M-1eJ) eb Te B eANDe TM(MtleJ) eL Te B) GO TO 13 


LECTM(M—-1¢J3) ol Te B eANDe TM(M41eJ) ¢GTe B) GO TO 15 
TFC TM(M—-1eJ) eGTe B eANDe TM(M#+1eJ) oh Te 8) GO TO 17 
GO To 19 

IF(NL elE&e 2) GO TO 9 


M-NPB)7e1195 
fe oa ik F eANDe TM(M—12J) eL Te BeANDe TM(M4+16J) eGTeF 


1eANDe TM(Nt1L 6 J) elite) GORtG 51 
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IF (CM eGTe NPB eANDe TM(M5J) «GTe 8) GO TO 31 
¢ IF(M eGTe NPB eANDe TM( oJ) eL Eo BIGO TO 37 
7 TFCTM(M—1LeJ) eGTe F eANDe TM(M—-10J5) eh Te B eANDo TM(M41,J) oGTe F 
1eANDOs TM(M#1eJ) eLTe B) GO TO 43 
IFC TM(MsJ)—-B) 39239,33 
9 TF(M e¢GTe KO eANDs TM(M—15J5) «GTe F eANDe TM(M—19J) eh Te B eANDe 
LTTM(M4#1.J3) eGTe F eANDe TM(MtieoJ) eb Te 8B) GO TO 43 
TF(M «GTe KO eANDe TM(M,J) 2GTo F) GO TO 33 
IFCM eGTe KO eANDe TM(MeJ) eLEo F) GO TO 39 
11 TFCTMOM=-15J5) eh Te B eANDe TM( Mtl, J) el Te 8B) 60 TO 21 
IFCTM(M—19J) el Te B eANDe TM(M4+15J) eGTe 8B) GO TO 23 
IFC TM(M—10J) eGTe B eANDe TM(Mt15/) eL Te 8) GO TO 25 
GO TO 27 
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TEMPERATURES AY NODAL POINTS AT INTERFACE BOUNDARIES ARE 


CALCULATED. 


13 TM(Ms-Jt1)=TM(Me J) ® (1 e—CCF/SCEF-—COF/CEF ) +TM(M—- 155) *( CCF /SCEF) +TM(M4F1 5 
1J)*(COF/SCEF) ; 
GO TO 61 

15 TMOMs U4+1)=TMOM0 J) & (1 eo —CCE/CHF-CO/ CHF ) + TM(M—1 0 J) #( CCF/SCHFE D+ TM(MF1 oJ 
1)* (CD/CHF) 
GO TO 61 

17 TM(M,J41)=TM( Me J) ¥( 1 eo -CC/CIF—COF/C IF) tTM(M—1 9 JS) *(CC/CIF D+ TM(M41 53) 
1*(COF/CIF) 
GO TO 61 

19 TMOMe S+L)=TM(Me J) (Le -CC/CE-CD/CE) t+TM(M—1 465) #(CC/CE) tTM(M4+1.I)%(CO 
1/CE) 
GO To 61 

21 TM(MeJ#1)=TM(MeJ) # (1 --COF/CAF-CBF/CAF D+ TM(M—19I)*( COF/CAF) + TM (M41 » 
1J)*(CBF/CAF) 
GO TO 61 

23 TM(Ms Jt1)=TM( Me J) * (Le -COF/SCFF-CB/ CFF) +TM (M—1 6 J) #( COF/CFF ) +tTM( M41 9 J5 
1) *(CB/CFF) 
GO TO 61 

25 TM(MsJ+1)=TM( Me J) # (1 e—CD/CGF—-CBF/CGF) +TM(M—1 6 JS) *(CO/CGF) TM (MEL J) 
1*(CBF/CGF) 
GO TO 61 

27 TMUMe J+1) = TMC Moe J) ® CL e—CO/CA—CB/ CA) TM (M1 9 JS) RC CO/CA)D TM (M4159 *( CB 
1/CA)D 
GO TO 61 

29 NZ=1 
GO TO 35 

31 NZ=3 
GO TO 35 

33 NZ=2 


TEMPERATURES AT NODAL POINTS IN A FROZEN OR UNFROZEN ZONE 


WITHIN A LAYER ARE CALCULATED. 
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Nanna 


35 TM(Me J#1)=(RUCNZ) #TM( M1 9 J) 4 (1 eo —-20*RU(NZ) )*TM( Ma J) RUC NZ) OTM (MHL 
1J) 
GO TO 61 

37 NZ=3 
GO TO 41 

39 NZ=2 

41 TMC Me J#1)=CRE (NZ) ¥TM(M—1 05) 40 10 -20%RFE (NZ) )*#TM( Me J) FRE (NZ) ¥TM(MEL » 
lJ) 
GO TO 6l 

43 NZ=2 


TEMPERATURES AT NOOAL POINTS WITHIN A ZONE OF FREEZING ARE 
CALCULATED. 


IF(PERWF(NZ) eEQe 1006) GO TO 45 
IF(PERWF(NZ) e€Qe 0.0) GO TO 47 
GO TO 49 
45 RWF=RF(NZ) 
GO TO 59 
47 RWF=RU(NZ) 
GO TO 59 
49 BF= (SLOPE (NZ) #.434*ALOG(TM(MsJ))481(NZ))/( 00434) 
CBWF=EXP( BF) 
RWE=((CKUCNZ) +KF (NZ) )/20) 70 CBWE*DDEN(NZ) )) XDELT/( (DELX/ 126) #*2) 
GO TO 59 
51 NZ=4 
IF (PERWF(NZ-1) e£Q- 1006) GO TO 53 
IF(PERWF(NZ—-1) eEQe 060) GO TO 5S 
GO TO S7 
53 RWF=RF(NZ—-1) 
GO TO 59 
55 RWFE=RUC(NZ-1) 
GO TO 59 
S57 BF=(SLUPE(NZ—1) #6434¥*ALOG( TM( Me J) )#BI (NZ—-1)970 06434) 


CBWF=EXP( BF) 
RWE=( CC KUCNZ—-1) #KF (NZ—1))720)97(0CBWF*ODEN(NZ—1))) *®DELT/( (DEL X/126) 


1**#2) 
s9 TMCOMs J#1) =RWE*XTM(M—1L oJ) tC Lo -20*RWE DE TM(Me JD tRWE*TMOMEL JS) 
61 CONTINUE 

SDEL=SDEL+#*tDELT 


THE TIME INTERVAL AT WHICH THE TEMPERATURES ARE TO 6E 
WRITTEN IS SELECTED.(THE TIME INTERVAL SELECTED IS 2 HOURS.) 


IF(SDEL eGEe 2.) GO TO 63 
GO TO 67 
63 HR=HR+2 
LLA=1 
CALL OUTPUT 
65 SOEL=0.20 
LLA=0 
67 CONTINUE 
RETURN 
END 
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SUBROUTINE GRAPHS 


SUBROUTINE GRAPHS 
CREEK EEE KEK KE 


THIS SUBROUTINE CONVERTS THE TEMPERATURES TO BE PLOTTED FROM 
ARRAY TO VECTOR NOTATION AND CALLS FOR THE PLOTTING OF THESE 
TEMPERATURES 


COMMON KU(3)2KF (3) 2CU(3).CF(3)2.00EN(3) »wC( 3) -TOL( 3) sRUC(3)sRFC 3) 
1CWF (4) oe WOEN( 3) »sNOP( 3) »sHS(29) 3H( 200) + TA( 30) sTM( 1153200) sPERWF(2) » 
2TAT(200) eTITLEL(20),TITLE2(20) >» IDEEP( 30)+eSLOPE( 3) »B81(3) »G( 30)» 
3BUF (2048) sDELX eKOxsKMZeNL »NOBCeoNPBslLNPHsNPSoLNPS oKOAs TOTAL» 
4L AST oL KL oe KOLKL oKLK sKOLKL1 eKOl » CAs CAF eo CFF » CGF CBs CBF CC oe CCF s CD »CDF » 
SCE » CEF » CHF oCIF »BoFsDELT eJJJoLitl eRRoIITI aL lL P»rsAxsEMsUCsNOTI 9X 
6TYME( 30) eMeaNOeoLLBeHRelLLAs Jo INC 

IFC ITI eGTe 19 GO TO 2 

CALL GRAPHI 

CONTINUE 

JUDE=DELX 

00 10 M=1 ,LAST 

IF(M eF£Qe 1) GO TO 4 

1=( IOEEP(M) +JUDE )/JUDE 

GO TO 6 

f=1 

KK=1 

TU=1./0ELT 

JU=TU 

OO 8 J=1eNOTI».JU 

GCKK)=TM(T 2 J) 

KK=KK+41 

CONTINUE 

KK=KK-1 

CALL GRAPH2 

CONTINUE 

RETURN 

END 
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SUBROUTINE GRAPHI 
SUBROUTINE GRAPHI 
EEK REECE EK EAR EK KK K 


SUBROUTINE GRAPH1 ESTABLISHES A GRID ON WHICH THE TEMPERATURES 
ARE PLOTTED. 


COMMON KU(3) sKF (3) eCUC3) »CF(3) sSODEN(3) »WC( 3) »eTOL( 3) »RU(3) »RFL3)5 
ICWF(4) sWOEN(3) »NO PC 3) »HS( 29) 0H( 200). TAC 30) 5 TM( 1155200) sPERWF( 2) 6 
2TAT(200) oeTITLEL (20) »TITLE2(20). IDEEP( 30) »SLOPE(3).81(3)+G(30)> 
3BUF (2048) »DELX » KO »KMZ » NL »NOGC pNPBeL NPB NPS oLNPS oKUOA sv TOTAL ® 
4LAST »LKL sKOLKL oKLK sKOLKL1 oKO1 » CAs CAF o CFF s CGF e CBe CBF s CC eo CCF sy CO sCDF » 
SCE » CEF » CHF oCIF oBoF sDELT pep JJJolLLitsRRelII oll PrsAsEMsUCeNOTI 9X» 
6TYME(30) eMeNOoLLBsHRelLLAs Js INC 


ALEN — THE LENGTH OF THE PLOT. (SCALE33 INCHES HORZe=1 DAY) 
THE MAXIMUM NUMBER OF DAYS THAT CAN BE PLOTTED IS 
. CONTROLLEO 8Y XLIMIT LN THE MAIN PROGRAM. 


ALEN=ND¥*3 
CALL PLOT(4.00.0-¢-21 .-3) 
CALL GRID(0c0»0e¢0sALENs 20001513) 


TITLEL PLUS STRUCTURE DESCRIPTION IS WRITTEN. 


CALL SYMBOL (120 2025e0s0eS5eTITLE1»0203860) 

CALL SYMBOL ( 3280224030 0035% 14H4 INCH SURFACE 200214) 
CALL SYMBOL (3280s 230320¢e35e18H6 INCH BASE COURSE 50-0:18) 
CALL SYMBOL (3¢8052203200e35e13HSAND SUBGRADE20-0013) 


THE COORINATES OF THE GRIO ARE LAGELED FROM —-SO DEGe F TO +80 
DEGe Fe (SCALE22 INCHES VERT.=10 OEGREES F) 


TEMP=-50.0 

Y¥=-0Oe2 

OO 14 [=1-5 

CALL NUMBER (—1 0402 V¥ 000405 TEMP90e0s~1) 
Y=Y¥+2.-0 

TEMP=TEMP4106¢ 

CONTINUE 

CALL NUMBER (—02602V20040+ TEMP 20 00 e-1) 
DO 16 [=18 

Y=Y+2.0 

TEMP=TEMP4106¢ 

CALL NUMBER (—160¢00 0 V2 00400 TEMP e 000e-1) 


CONTINUE 
CALL SYMBOL (- 2 000 70702006002 1HTEMPERATURE (DEGe Fe) s90e0s21) 


CALL PLOT (020210+00-3) 
RETURN 
END 
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SUBROUTINE GRAPH2 


SUBROUTINE GRAPH2 
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SUBROUTINE GRAPH2 PLOTS THE CALCULATED TEMPERATURES AT 
PREVIOUSLY DEFINED DEPTHS AS A CONTINUGUS FUNCTION OF TIMEe 


COMMON KU(3) sKF (3) 2CU(3)+CF(3)eODEN(3) »WC( 3) TCL (3)2RU(3) + RFCS)» 
1CWF (4) »pWDEN(3) »NOP (3) »HS( 29) eH( 200) 2 TAC 30) 0 TM( 1155200) sPERWE(2) >» 
2TAT(200) ,TITLEL (20) eTITLE2(20) -IDEEP( 30) »SLOPE(3).81(3)-+G(30), 
3BUF (2048) »sDELX » KO» KMZ » NL »NOBC »NPBeLNPBeNPS »oLNPSoKUAs TOTAL » 
4L AST oLKL »KOLKL oKLK eKOLKL 1 oKO1 » CAs CAF o CFF 9 CGF os CB ep CUF » CC 0 CCF o CD eCOF 9 
SCE e CEF o CHF eCIF sB oF oDELT eJJJollLits RRo Il LolLLPsAr»sEMsUCsNOTIsXos 
6TYME( 30) eMeND oL LB oHR epLLAsJeINC 

INTEGER X 

IF(M «EQe 1)GO TO 18 

GO TO 20 


THE DATE AND TIME OF DAY ARE WRITTENe 


CALL SYMBOL (0 2650-77-90 e0elLSeTITLE2 2020080) 
CALL SYMBOL (-0.102—-8220200e16s2H24 5000 22) 
CALL SYMBOL (0 e952—-S8e202 0 lL6e1HBs0e0el) 
CALL SYMBOL (1690 1—-8e209001l692H1690 2002) 
CALL SYMBOL (2090 9—-8020000e1622H24202022) 
CALL SYMBOL (1 632-8035» O00ol2s4HHRS02 00004) 
CONTINUE 

TYME (1)=060 

DO 26 [=2X 

IF (X e€Qe 13) GO TO 22 

TYMECT)=I-1 

GO TO 24 

TYME(L)=( I-11) *2 

CONT (NUE 

CONTINUE 

TYME(X#+1)=020 

TYME(X+2)=8-0 

G(X+1)=0-0 

G(X#+2)=5-0 


TEMPERATURES AT A GIVEN DEPTH ARE PLOTTEDe 


CALL L INE( TYME 0GoX 010000) 
IF(M «FQ. LAST) GO TO 28 
GO TO 30 

CALL PLOT(320:0-0+-3) 
RETURN 

END 
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SUBROUTINE TEMPERATURE 
SUBROUTINE GRAPH 2 
SUBROUTINE AIRSUN 


10 


SUBROUTINE GRAPHS 
: TT 


2) 


MAIN 


SUBROUTINE GRAPH | 


SUBROUTINE OUTPUT 


PROGRAM 


SUBROUTINE BOUNDARY 


SUBROUTINE DEPTH 
B=|SUBROUTINE NODAL 


FIGURE A-| GENERALIZED FLOW DIAGRAM OF TEMPERATURE 
PREDICTION PROGRAM. 
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FIGURE. 


A-2 


Input Variables 
are specified 


Subroutines Boundary, Nodal 
and Depth are called 


Initial Temperature 
gradient is specified 


: Counters LLA, LLB and 
LLP are defined 


5 
DO 150 III = 1,ND 
Climatic inputs are 
specified 


Subroutines Airsun, 
Temperature and Graphs 
are called 


7 
DO 12 M = 1,KOA 
Temperature gradient 
is defined 


END 


FLOW DIAGRAM FOR MAIN PROGRAM. 
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FIGURE. 


Main Program 


Diffusivity of each 
layer is calculated 


Latent heat 
released upon freezing of 
the lower layers is 
calculated 


Constants of the 
finite difference 
equations are computed 


Temperature-apparent 
specific heat relationships 
of the lower layers are 
calculated 


Return to Main Program | 


A-3 FLOW DIAGRAM FOR SUBROUTINE BOUNDARY. 
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Temperatures 
at interfaces 


FIGURE . 
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DO 67 J = 1,NOT 


Surface temperature 


; , 


Temperature at depths 


Temperatures in 
a freezing layer 


Temperatures 
within a layer 


3 Time interval for temperature 
display is specified 


Call Output : 
Return to Main Program 


A-4 FLOW DIAGRAM FOR SUBROUTINE TEMPERATURE. 
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APPENDIX B 


COMPUTER PROGRAM FOR THERMAL STRESS 


PREDICTIONS USING PSEUDO-ELASTIC ANALYSES 


Bl 
APPENDIX B 


COMPUTER PROGRAM FOR THERMAL STRESS PREDICTIONS 


USING PSEUDO-ELASTIC ANALYSES 


Pols - Introduction 
This appendix contains details of a computer program for calcu- 
lating thermal stresses in an asphaltic concrete subjected to a known 
temperature change. For these stress computations pseudo-elastic beam 
and slab analyses are used. 
The stress equation for the pseudo-elastic beam analysis is 
t 
o, (t) = -f SCAG Waa tigdi ay) (B-1) 
t 
O 
where S(At,T) is the time and temperature dependent stiffness modulus 


and a, (T) the coefficient of thermal expansion. To evaluate stress 


numerically Equation (B-1) is written as 
= Ao). B-2 
0, (t;) Galt, qdesi Cos): (B-2) 


where 
ti 


(Ads )iieteay - if S(At ,T)a,(T)dT(t) (B-3) 


Lie 


Using the trapezoidal rule, with stiffnesses evaluated at times t; and 
tj_, and a loading time equal to the time step t; - tj_,, and assuming 


a) temperature independent, the stress increments in Equation (B-3) are 


computed by the expression 
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B2 
a 
(Ao); = > (81+8}-.) (Tq-Ti-1) (B-4) 


Approximate pseudo-elastic slab stresses are obtained by 
multiplying Equation (B-2) by the factor of 1/(1-v), where Poisson's 


ratio, v, is assumed time and temperature independent. 


B.2. Program Operations 


The computer program consists of a main program and three sub- 
routines, titled Readin, Linear and Stress. Operations within each of 


these four components are given in the following subsections. 


Be2.d. Main Program 


A flow diagram for the main program is shown in FIGURE B-1. 
Within this program all subroutines are called and all output commands 


are introduced. 


B.2.2. Subroutine Readin 


In this subroutine all parameters necessary for the numerical 
solution of Equation (B-2) are specified. These parameters are 


defined in the program listing. 


B.2.3. Subroutine Linear 


Linear relationships are assumed between the logarithms of 
stiffness and the logarithms of reduced time values, and between the 


logarithms of shift factor values and temperatures, specified in sub- 


routine Readin. In subroutine Linear the coefficients of these rela- 


tionships are computed. Such relationships enable the stiffness of the 


asphaltic concrete to be readily evaluated. 
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Temperatures (TT) at time increments of tj - tj; (DELT Pare 
computed. For these computations, a linear temperature-time relation- 


ship is assumed over the time period INCTIM defined in the main program. 


B.2.4. Subroutine Stress 

In this subroutine thermal stresses (ST) are calculated at time 
increments of DELT. Using a loading time equal to DELT and knowing the 
temperature TT at each time step stiffness values are computed. Upon 
substituting these values into Equation (B-4), or the corresponding 
pseudo-elastic slab equation, the stress increment (Acx); (STR2) is 
calculated. In the program listing the initial stress o,(t,) (ST) is 


equated to zero. 


Program operations return to the main program where tempera- 
tures and thermal stresses are written in time increments of NT defined 


in subroutine Linear as (INCTIM) (60) /DELT. 
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THIS PROGRAM USES A MODIFIED ELASTIC APPROACH TO 
CALCULATE THE THERMAL STRESSES IN AN ASPHALTIC CONCRETE 
SUBJECTED TG A KNOWN TEMPERATURE FIELD. 


STRESSES ARE CALCULATED AT EQUAL TIME INTERVALS 
OVER THE TOTAL TIME PERIOD IN WHICH THE TEMPERATURE FIELO 
IS DEFINEDe 


THE PROGRAM CAN BE USED TO DETERMINE THE THERMAL STRESSES 
ASSUMING ONE OF THE FOLLOWING CONOITIONSS 
1) RESTRAINED BEAMsE=FUNCTION( REDUCED TIME) 
2) RESTRAINED SLAB8,E=FUNCTION(REDUCED TIME) 


MAIN PROGRAM 
€HRKEREREEE REE 
COMMON TITLEL( 20) oT(15) sATC15) 0E (15) XT 015) sAL (15) sA2(15) eB1l{15)>5 
1B2(15) sTEMP( 1300)» TT(3000) 2€1( 3000) »ST( 3000) sNUMDAY»NUMEXI »NUMATT » 
2DELT » ALPHA eUs CODE s NEXE1 eNATTI oNNoNT oMT eSTR1 ®NUMINCo INCTIM 
CALL READIN 
CALL LINEAR 
DO 1 MT=1.NN 
CALL STRESS 
CONTINUE 
WRETE( 60 390 TITLELC TE) o1=1220) 
FORMAT (6X%» 20A4) 


TEMPERATURES AND STRESSES ARE WRITTEN AT TIME INCREMENTS 
CORRESPONDING TO THOSE AT WHICH THE TEMPERATURES WERE SPECIFIED. 


DO 7 J=1 »sNNeNT 

WRITE (605) TT(JS) o STC JS) 
FORMAT(12X%e2F1001) 
CONTINUE 

STOP 
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SUBROUTINE READIN 


SUBROUTINE READIN 
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THIS SUBROUTINE READS IN ALL DATA NECESSARY FOR THE 
GENERATION OF THERMAL STRESSES. THIS DATA CAN BE OIVIDED INTO 
THE FOLLOWING THREE CATEGORIESS 
THE ASSUMPTION OF A RESTRAINED GEAM OR SLAB,» 
CONDITIONS 1 OR 2 MENTIONED PREVIOUSLY. 

MATERIAL PROPERTIES « 
TEMPERATURE AT A GIVEN DEPTH IN THE ASPHALTIC 
CONCRETE AS A FUNCTION OF TIMEe 


1) 


2) 
3) 


COMMON TITLE! (20) 6T(15) sATCISI SECIS) sXIC15)eA1(15) »A2(15),681015)> 
182015) 6 TEMP(1300)sTT(3000) -E1( 3000 )eST( 3000) sNUMDA YsNUMEXI»NUMATT y 
2DELT»s ALPHA» Us CODE sNEXI1 »sNATTIL eo NNoNT eM eo STRIeNUMINCs INCTIM 


CODE 


U 


NUMDAY 


NUMEXI 


NUMATT 


INCT IM 


OELT 


ALPHA 


E 


XI 


AT 
T 
TEMP 


THE INPUT PARAMETERS ARE DEFINED AS FOLLOWS: 


CODE=1-e0 DENOTES CONDITION te 

CODE=2-e0 DENOTES CONDITION 2e 

POISSONS RATIGes GNLY ITF CONDITIGN 2 IS SPECIFIED 
A VALUE OF POISSONS RATIO MUST BE READ INe 

THE NUMBER OF DAYS IN WHICH THE TEMPERATURE 
REGIME IS READ INe 


THE NUMBER OF VALUES SPECIFIED FROM THE STIFFNESS MODULUS-— 


REDUCED TIME RELATIONSHIP. 
THE NUMBER GF VALUES SPECIFIED FROM THE SHIFT FACTOR— 
TEMPERATURE RELATIONSHIP. 

THE TIME INCREMENT IN WHICH THE TEMPERATURE 

REGIME IS SPECIFIEOD.( HRS) 

THE TIME INCREMENT IN WHICH THE THERMAL STRESSES 

ARE TO BE CALCULATED. (MINS) 

THE COEFFICIENT OF THERMAL EXPANSION AND CONTRACTION 
OF THE ASPHALTIC CUONCRETEs(IN/IN/DEG F) 

THE STIFFNESS MODULUS OF THE ASPHALTIC CONCRETE«(PSI) 
THESE VALUES ARE TO BE READ IN IN DESCENDING ORDER 

OF MAGNITUDE. 

REDUCED TIME.(SECS) THESE VALUES ARE SPECIFIED IN 
ASCENDING ORDER OF MAGNITUDE’ AND MUST BE COMPATIBLE 
WITH THE E VALUES. 


SHIFT FACTORS» 
TEMPERATURES COMPATIBLE WITH THE SHIFT FACTORS.(DEG F) 


THE TEMPERATURES AT A GIVEN OEPTH WITHIN THE 
ASPHALTIC CONCRETEe(DEG F) 


READ(S,e1) (TITLEI( I) f= 220) 


1 FORMAT (20A4) 
READ(5»2)CODE 
2 FORMAT(1FS-0) 
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IF (CODE-240)80408 
READ(5e6)U 

FORMAT(1F562) 

GO TO 10 

U=0.0 
READ(S»12)NUMDAY »NUMEXI »NUMATT» INCTIMsDELT sALPHA 
FORMAT (4110 .1F 10025 1E Gel) 
READ(5e14) (ECL) »f=1 eNUMEXI) 
FORMAT(11E7 62) 
READ(5s16)¢(XE(1)+T=lsNUMEXI) 
FORMAT(1LLE7*1) 
READ(5518) CATCL) s T=leNUMATT) 
FORMAT(9E7e2) 
READ(S5s 20) (TCL) of =t eNUMATT) 
FORMAT (9F8.0) 
NUMINC=24/INCTIM 
NUMTEM=NUMINC#NUMDAY #1 

DO 24 L=1».NUMTEM 
READ(5 22) TEMP(L) 
FORMAT (1F 1020) 

CONTINUE 

RETURN 

END 
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SUBROUTINE LINEAR 


SUBROUTINE LINEAR 
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A LINEAR RELATIGNSHIP £S ASSUMED BETWEEN LOG(E) AND 
LOG(XT)*¢AND LOG(AT) AND T VALUES. THIS SUBROUTINE CALCULATES 
THE A ANO 8 COFFICIENTS OF THE RESULTING LINEAR EQUATIONS. 
ALSOsTEMPERATURES OF THE ASPHALTIC CONCRETE AT TIME 
INTERVALS OF DELT ARE CALCULATED. 


COMMON TITLEL( 20) sT(IS)sAT(LSISGECLSVeXT(15 VeAL (15) sA2Z{15) sB1C15)6 
182015). TEMP(1300)5TT( 3000) ,E1( 3000 )2ST( 3000) seNUMDAYsNUMEXI »>NUMATT » 
2DELT »s ALPHAs Us CODE seNEXIL »NATTLeoNNo NT oMTe STRIeoNUMINCs INCTIM 

NEX L1=NUMEXI-1 

NATTL=NUMATT-1 


THE A ANDO 8 COEFFICIENTS FOR THE LOG(E)—-LOG(XI) 
EQUATIONS ARE EVALUATED. 


00 1 J=1.NEXI1 
A1(J)=-1LeO*ALOGIOCECJISE( I+1)9 

B1¢€ J) =ALOGIO(E(J))—AIL CI) FALOGLIO(X ICI) D 
CONT I NUE 


THE A ANO 8B COEFFICIENTS FOR THE LOG(AT)—-T EQUATIONS 
ARE EVALUTEDe. 


00 3 J=1sNATTI 
A2(J)=CALOGIOCAT( JS) JD—ALOGIOCATCJ#+1)9)70T( 5 )-TO 41) 9 
B2(J) =ALOGIOCAT( JI) )-A2( ID FT CS) 

CONT INUE 

STR1I=0.0 

TIME=INCTIM 

NT=(TIME*60 -)/OELT 

NTL=NT +2 

N=NUMDAY* NUMI NC 

NN=N#€NT#1 


TEMPERATURES (TT) AT INTERVALS OF DELT ARE CALCULATED. 


TTC L)=TEMPCI) 

L=2 

OO 7 M=1eN 

00 S LM=2,NT1 

TT(L)=TEMP (OM) +0 CTEMPCM#L)-—TEMPOM) DZNT)*CLM—21) 
L=L+l1 

CONTINUE 

CONTINUE 

RETURN 

END 
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SUBROUTINE STRESS 


SUBROUTINE STRESS 
MM Fe de te a he ek ee 


SUBROUTINE STRESS CALCULATES THERMAL STRESSES AT TIME 
INTERVALS OF DELT. 


COMMON TITLEL(20)6T(15) sATCIS) .EC15) »X1415) eAL (15) 5A2(15) »sBLC15)> 
182015) »-TEMP(1300)sTT( 3000) .E1{ 3000) 6 ST( 3000) sNUMDAYs NUMEXI » NUMATT » 
2DELT »sALPHA.U,» CODE sNEXE1 »sNATT1 »NNoNToMTeSTR1»NUMINC» INCT IM 

DO 1 L=1.6NATTI 

LECTTOCMT) eh Te TIL) sANDe TTCMT) «GE. T(L+1)) GO TO 3 

CONTINUE 

A=AT(1) 

GO TO 5 

A=EXP((A2(L) *TTOMT)4#820L))700434) 


THE REDUCED TIME (XI1) AS A FUNCTION OF DELT AND 
TEMPERATURE IS CALCULATED. 


XI 1=(DELT*60.)7A 

{FCXT1 ek Te XI(1)) GO TO 9 

IF(XI1 eGT. XTC NUMEXIT)) GO TO 11 

DO 7 L=1eNEXT1 

IFOXEL eGEe XI(L) eANDe. XI1 LE XI(L#1)) GO TO 13 
CONTINUE 

ELCMT)=EC1) 

GO TO 15 

E1(MT)=ECNUMEXT D 

GO TO 15 
EL(MT)=EXPC(B1 (LL) FAL(L) €ALOGIO(C XI 1) 970.434) 
IF (MT eEQ-e 1) GO TO 17 


THE THERMAL STRESS INCREMENT (STR2)+2TOGETHER WITH THE 
TOTAL STRESS (ST) ARE CALCULATED. 


STR2=(100/(100-U) ) #( ALPHA/2 00) HCE L(MT) FE LC MT— 1) RC TTC MT~1L)—TTCMT)D 


STRI=STR2+tSTRI 
ST(MT)=STRI 
GO TO 19 
ST(MT)=020 
RETURN 
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FIGURE. B-| 


Call Subroutines 
Readin and Linear 


Call Subroutine 
Stress 


Write Title 


DO! 7 Jo = TV SNNSNT 


Write Temperatures 
and Stresses 


FLOW DIAGRAM FOR MAIN PROGRAM. 
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APPENDIX C 


COMPUTER PROGRAM FOR THERMAL STRESS 


PREDICTIONS USING VISCOELASTIC ANALYSES 


WITH VARIABLE TIME INCREMENTS 
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APPENDIX C 


COMPUTER PROGRAM FOR THERMAL STRESS PREDICTIONS USING 


VISCOELASTIC ANALYSES WITH VARIABLE TIME INCREMENTS 


Cod. Introduction 

This appendix contains details of a computer program developed 
for calculating thermal stresses in an asphaltic concrete with tempera- 
ture dependent viscoelastic properties. The program can be used for 
stress computations involving three different analyses. These analyses 


are designated as: 


1) Viscoelastic slab, 
2) Viscoelastic beam, and 


3) Approximate viscoelastic slab. 


Assuming the coefficient of thermal expansion of an asphaltic concrete 


to be a constant value, stress equations for the viscoelastic slab and 


beam analysis are 
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C2 


where 
Qo = coefficient of thermal expansion 
R(€) = a modulus function 
S(&) = stiffness modulus 
T = temperature 
< = time 


Evaluating Equations (C-1) and (C-2) by the trapezoidal rule, results in 


N 

Ox(t) = - 90D) [RCE+E,) + R(EEy_ 1) CUy-Ty_2) (C-3) 

i=] 
and 
N 
O 

ox(t) =- + [S(é-£4) + $(E-€4-.)] (T4-TH-1) (C-4) 

i= 4 


where N is the number of time steps from 0 to t. 


Approximate viscoelastic slab stresses are obtained by multi- 
plying Equation (C-4) by the factor 1/(1-v), where Poisson's ratio, v, 
is assumed constant. 


For stress computations at a time t the viscoelastic stress 
equations require numerical integration over the entire previous 
temperature history to which the asphaltic concrete was subjected. 
Therefore, when stresses are to be predicted over extended periods of 
time long computer runs and large storage requirements are involved. 
Since the most recent temperature history is the most significant for 


stress evaluation, the following scheme of variable time increments is 
c > “ 


included in the program. 
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1) For times up to (t-24) hours a time increment of 2 hours is 


employed. 


2) For the following 22 hours a time increment of 15 minutes 


is employed. 


3) For the 2 hours immediately prior to t, a time increment 


of 5 minutes is employed. 


C.2. Program Operations 


The program consists of a main program and seven subroutines. 
Computations within each of these components are outlined in the 


following subsections. 


C.2.1. -Main Program 


A flow diagram for the main program is shown in FIGURE C-1. 
Within this program all subroutines are called and all output state- 
ments are given. In the listing, output display consists of 


temperatures (TT) and computed stresses (ST) at 2 hour time intervals. 


C.2.2. Subroutine Readin 


In this subroutine all parameters necessary for the numerical 
solution of the viscoelastic stress cquations are specified. These 


parameters are defined in the program listing. 


C.2.3. Subroutine Linear 


For the stress computations, stiffness values at times J and 
temperatures TT are evaluated by assuming linear relationships between 


logarithm of stiffness and logarithm of reduced time values, and 
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C4 
between logarithm of shift factor values and temperatures specified in 
subroutine Readin. In subroutine Linear the coefficients of these 


relationships, designated as Al, Bl, A2 and B2,,are calculated, 


C.2.4. Subroutine XI15M 

In this subroutine the reduced time vector ~ (XI15M) ee ee 
with the input temperatures (TEMPS) is calculated at time Steps Of ts 
minutes. For this computation a linear temperature-time relationship 
is assumed over the time period INCTIM. Temperatures at 15 minute 
intervals are computed and used to evaluate the corresponding shift 
factors AT15. Then, using the trapezoidal rule, the reduced time 


vector XI15M is calculated. 


Co2m5 nr Subroutine 'Str2 


In this subroutine thermal stresses are computed for times up 
to (t-24) hours using time increments of 2 hours. The limit of this 
time period is defined by the parameter LM2, while, the limits of the 
time periods over which 15 and 5 minute time increments are used are 
defined by the parameters LM15 and LMXI, respectively. Stiffness 
values El(J) are determined at reduced time values Of ((e-6")) 
[XI15(LMXI)] - [XI15(J)], where J increases in increments of 2 hours, 


and the thermal stress ST1 at the time limit LM2 is computed. 


€.2.6, ‘Subroutine Stris 


In subroutine Str15 thermal stresses are calculated over the 


time period LM2 to LM15. During this period, a time increment of 15 


minutes is employed for stress evaluation and operations proceed in a 


manner similar to that described for subroutine Str2. 
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C.2.7. Subroutine XI5M 

The reduced time vector —€ (XI5M), for temperatures over the 2 
hour time period immediately prior to the time of stress prediction, is 
computed. For this computation a linear temperature-time relationship 
is assumed and temperatures at 5 minute intervals together with the 


corresponding shift factors (ATS) are determined. The reduced time 


vector XI5M.is then computed. 


C.2.:8.. / Subroutine stro 

Using the reduced time vector XI5M and temperatures computed in 
subroutine XI5M, thermal stresses during the 2 hour time period are 
calculated. The predicted thermal stress ST at time LMXI represents the 
sum of stress increments computed in each of the three stress sub- 
routines. Operations return to the main program where calculations 


commence for stresses at subsequent times and temperatures. 
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IN THIS PROGRAM VISCOELASTIC ANALYSES ARE USED TO 
CALCULATE THERMAL STRESS IN AN ASPHALTIC CONCRETE SUBJECTED 
TO A KNOWN TEMPERATURE FIELD. 


STRESSES ARE CALCULATED USING A NUMERICAL 
INTEGRATION PROCEDURE WHICH INCORPORATES TIME STEPS 
OF 2 HOURS»®15 MINUTES AND 5 MINUTES. 


THE PROGRAM CAN BE USED TO DETERMINE THE THERMAL STRESSES 
ASSUMING GONE OF THE FOLLOWING CONDITIONS; 
1) RESTRAINED BEAM,E=FUNCTION{REOUCED TIME) 
2) RESTRAINED SLAB »E=FUNCTION(REDUCED TIME) 
3) RESTRAINED SLAB»sR=FUNCTION(REOUCED TIME) 


M MAIN PROGRAM 


Ae He he ok Kk tk ek 


COMMON TITLEL (20) 2T(15) sATI(ISI6E(15)>XL015) oA1(1S5) sA2015)581015)>0 
182015) .TEMP(1300).TT(3000) »ATIS( 3000) »XT15(3000) »E1(3000)»ST{(3000) 
22TA(50) sAT5(50) 2X15(50) »NUMOAYsNUMEXT »>6NUMATT »NUMINCsINCTIMsAL PHA,» 
BNATTI oNEXT 1 eCODE oUseSTl ol MoMT sLM2eLMISoLMX I 

DOUBLE PRECISION DXIsXI15eXI5 


CALL 
CALL 
CALL 
oo 1 
CALL 
CALL 
CALL 
CALL 


READIN 
LINEAR 
XIT15M 
MT=14 LM 
STR2 
STRIS 
XISM 
STRS 


CONTINUE 
WRITE(S es 3V(TITLEL( I) of =1 520) 


FORMAT( 6X 20A4) 
L=105 


TEMPERATURES AND STRESSES ARE WRITTEN 


ODO 7 MZ=14.eLM 
WRITE (655) TT(L) oST( NZ) 
FORMAT (6X 92F1061) 


L=L+8 


CONTINUE 


STOP 
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SUBROUTINE READIN 


SUBROUTINE READIN 
Me HR Re a i tee 


THIS SUBROUTINE READS IN ALL DATA NECESSARY FOR THE 
THERMAL STRESS PREDICTIONSe THIS DATA CAN BE DIVIDED INTO 
THE FOLLOWING THREE CATEGORIES; 

THE ASSUMPTION OF A RESTRAINED BEAM OR SLAB, 
CCNOITIONS 1+2 OR 3 MENTIONED PREVIOUSLY. 
MATERIAL PROPERTIES. 

TEMPERATURE AT A GIVEN DEPTH IN THE ASPHALTIC 
CONCRETE AS A FUNCTION OF TIME. 


1) 


2) 
3) 


COMMON TITLEL (20) oTC1S) eAT( LS) sEC 15) o XT (15) sAL (15) eA2(15) eBI(15)> 
1B2(15) 6 TEMP( 1306) + TT(3000) sATI5{ 3000) eX 1153000) sE1(3000)+ST{ 3000) 
29TA(50) 6 ATS(50) »X15( 50) »NUMDAY > NUME XI eNUMATT sNUMINC op INCTIMs ALPHAs 
BNATT1 sNEXI1 sCODEsUsST! sL Mo MT sLM2sL MLS eLMXI 

DOUBLE PRECISION DXIsXI15sXI5 

THE 


CODE 


U 


NUMDAY 


NUMEXI 


NUMATT 


INCTIM 


ALPHA 


3 


ne | 


AT 
¥ 
TEMP 


INPUT PARAMETERS ARE DEFINED AS FOLLOWS: 

CODE=1-e0 DENOTES CONDITION le 

COOE=2e0 DENOTES CONDITION 2e 

CQDE=3-0 DENOTES CONDITION 3e 

POISSONS RATIOe ONLY [TF CONDITION 2 IS SPECIFIEO 

A VALUE OF POISSONS RATIO MUST BE READ INe 

THE NUMBER CF DAYS IN WHICH THE TEMPERATURE 

REGIME IS SPECIFIED. 

THE NUMBER GF VALUES SPECIFIED FROM THE STIFFNESS MOOULUS-— 
REOQUCED TIME RELATIONSHIP. 

THE NUMBER OF VALUES SPECIFIED FROM THE SHIFT FACTOR— 
TEMPERATURE RELATIONSHIP. 

THE TIME ENCREMENT IN WHICH THE TEMPERATURE 

REGIME 1S SPECIFIEDe{HRS) 

THE COEFFICIENT OF THERMAL EXPANSION AND CONTRACTION 
OF THE ASPHALTIC CONCRETEe(IN/IN/DEG F) 

THE STIFFNESS MODULUS OF THE ASPHALTIC CONCRETE (PSI) 
THESE VALUES ARE SPECIFIED IN DESCENOING ORDER 


OF MAGNITUDE « 
RECUCEO TIMEe{SECS) THESE VALUES ARE SPECIFIED IN 


ASCENDING ORDER OF MAGNITUDE AND MUST BE COMPATIBLE 
WITH THE E VALUES. 


SHIFT FACTORS. 
TEMPERATURES COMPATIBLE WITH THE SHIFT FACTORS.(OEG F) 


THE TEMPERATURES AT A GIVEN DEPTH WITHIN THE 
ASPHALTIC CONCRETE«(DEG F) 


READ(Se LIC TITLEL( LT) ol=1220) 


1 FORMAT (20A4) 
READ (5e2)CODE 
2 FORMAT(1F5-0) 
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SUBROUTINE LINEAR 


SUBROUTINE LINEAR 
ROR HO A ek te tee took te ke 


A LINEAR RELATIONSHIP IS ASSUMED BETWEEN LOG(E) ANDO 
LOG(XI)sAND LOG(AT) AND T VALUES. THIS SUBROUTINE CALCULATES 
THE A AND 6 COEFFICIENTS OF THE RESULTING LINEAR EGUATIONSe 


COMMON TITLEL(20)6T(15) sATCIS) cEC LISI eXLCLSI2AL(1IS) sA2(15) BIC1IS)» 
162015) » TEMP(1300).TT{ 3000) sATI5( 3000) e-X115(3000) -£1(3000),ST( 3000) 


2eTA(SO) eATS(S0) 2X 15050) »sNUMDAY » NUMEX LeNUMATT »eNUMINCs INCTIMs AL PHA, 
3NATT1I »NEXTL sCODE eUsST 1 »sL Mo MT pL M2eLMISoLMXI 


DOUBLE PRECISION OXI eXL1i5.XI5 
NEX I L=NUMEX I-1 
NATTI=NUMATT—1 


THE A ANC 8 COFFFICIENTS FOR THE LOG(E)—LOG(XI ) 
‘EQUATIONS ARE EVALUATED. 


OO i J=1sNEXElL 
A1L(J)=—-1.0#ALOGIO(CE( JI SEC J+1)D 
B1(J)=ALOGIO(CE(I))—-AL( J) XALOGIO(XI(I)) 


CONTINUE 


THE A ANO 8&8 COEFFICIENTS FOR THE LOGCAT)I—¥F EQUATIONS 
ARE EVALUTEDe 


DO 3 J=1eNATTI 
A2( JS) =C(CALOGIOCAT(C I) D—ALOGLOCAT( JFL) DISC TOC II-TC S41 DD 


B2(J)=ALOGIOCATI I) —A2C ID*TC ID 
CONTINUE 

RETURN 
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SUBROUTINE XI15SM™M 


SUBROUTINE Xf15M 
RRRRER DEREK EH KE 


THIS SUBROUTINE CALCULATES THE REDUCED TIME VECTOR(XT15M) 
ASSOCIATED WITH THE INPUT TEMPERATURE FIELOe THIS VECTOR IS 
DETERMINED AT TIME INTERVALS OF 15 MINUTES OVER THE 
COMPLETE TIME PERIOD IN WHICH THE TEMPERATURE FIELD IS 
DEF INEDe 


COMMON TITLE? (20) 0T(15) 6ATL15)sEC15) oXI( 15) sA10 15) sA2015) oBILI5) © 
1B2015) oo TEMP( 1300) +TT(3000) 6AT15(3000)+X115(3000).E1(3000)sST( 3000) 
2% TA(SO) sATS(50) »X£5(S0) eNUMDAY»sNUMEXI sNUMATT sNUMINCe INCTIMs ALPHA 
SNATTLeNEXI1LsCODEsUoST1 sLMoMT sLM2sLMISoLMXI 

OOUBLE PRECISION DXIsXIt15exXI{5 

ST1=020 

LM=NUMDAY*12¢1 

TIME=INCTIM 

NT=TIME*60 e/156¢ 

NTI=NT+1 

N=NUMDAY *NUMINC 

TV¥(1)=TEMP(1) 


TEMPERATURES (TT) AT 15 MINe INTERVALS ARE CALCULATED. 


L=2 

00 3 M=1eN 

DO 1 MM=2,NT1 

TT(L) =TEMPIM)+¢ (TEMPO M41 )—-TEMPO(M) D/NT )*(MM—-1) 
L=L+1 

CONT [NUE 

CONTINUE 

NN=N#*¥NT#1 


SHIFT FACTORS (ATIS) AT 15 MINe INTERVALS ARE CALCULATED. 


DO 9 L=1.NN 

DO 5S J=1sNATTI 

IFC TT(L) oh Te TKI) 
CONTINUE 
AT15(L)=ATO1) 

GO TO 9 
ATIS(LI=EXPC(CA2Z(I) HTT (LL) +8205) 97064934) 


CONTINUE 


eANDe TT(L) eGEe T(It+l1)) GO TO 7 


THE REDUCED TIME VECTOR (x£15) IS CALCULATED. 


XT15¢(1)=0-90 
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DO 11 t=2.NN 

XTIS(LI=(9O0 6/2 cO) (1 COSATIS(LIF1-O/ATIS(L—-1))+XT1S5S(L—-1) 
CONTINUE 

RETURN 

ENO 


C10 


ihe a Vn! L : 
re : iv 
ti Pa ie: 
i 7 ¢ Ae 
pI Ohe 
' . eR ; - 
Ch ape aie Ae 
td vi ij 
APO j a 
{ : ’ 
Muah oY 
Van _ 
en 
Linden 
Vive | 4 ial 
A : 4 ; 
i} A? ald a) = 
‘ ' ite, “G 
‘ 


7 rm v a 
| bk Wes" 4 21 00 ; 
Ti < DION TME LCL Hd et TAN 2194 hee TAA. te Coe GN. G8ET= CITEDTK 
1 : " ~~ wn BUAITHOD 4 ae | 

) Ye, MAUTSR > 
a ana ee ( 
7 } i 7 


} 


ri 
{ 
: ne 
> 4 cs j 
- 
q 
it ‘ 
a _ 
1 
’ apy 
\ 4 : 4am 
| ¢ “«( 
7 
-— ag f ¥ 
7 , 
| ( We pe 
. es Od 


| c . pe , ieee, ‘ f , 
pr : tar : l ; it : Py jon ot * ' 


OoPAD AH) 
7 
7 


ANNA ANAANARAA 


ananAA 


10 


Cit 


SUBROUTINE STR2 


SUBROUTINE STR2 
Se Ok aR dk gk A tk 


SUBROUTINE STR2 CALCULATES THERMAL STRESSES DURING THE 
TIME PERIOD PRIOR TO (T-24HRS-) «THESE STRESSES ARE CALCULATED 
USING A 2 HOUR INTERVAL. 


COMMON TITLEL (20) oT (15) sATKC 15) 0E015) XT 015) SAL (15) sA2(15)2B1C15)¢ 
182(15) 6 TEMP( 1300) oe TT( 3000) eAT15( 3000) sX115(3000) .E1(3000)+ST( 3000) 
2+ TACSO) © ATS(50) »X15( 50) »NUMDAY »NUMEXI sNUMATT »NUMINC » INCTIMs ALPHA 
BNATT1 oNEXI1 CODE sUoSTI ol Me MT ol M2oL M15 sLMXI 

DOUBLE PRECISION DXI»X115eXI5 

LM2=(MT-13)*841, 

LM1S5=(MT—2) *841 

LMXI=(MT-1) #841 


THE CHANGE IN RKEOUCED TIME (OXI) AS A FUNCTION OF TIME 
ANDO TEMPERATURE IS EVALUATED. 


DO 10 J=1.lM258 

OXIT=XTiS(LMXII—XTIS(SI) 

IF(DOXEI eLEe XI(1))GO TO 4 

ODO 2 L=l»NEXT1 

ILF(DOXT eG&e XI(L) eANDe OXI eLEe XI(L+1)) GO TO 6 
CONTINUE 

E1€CJ)=ECNUMEXT D 

GO TO 8 

E1(CJ)=E(1) 


GO TO 8 
EL( J)=DEXP((B1(L) FAL(L) #OLOGIO( OXI ))702434) 


IF(J eEQe 1) GO TO 10 D 


THE THERMAL STRESS LNCREMENT (STR2) » TOGETHER WITH THE 
TOTAL STRESS (ST1)e0IS CALCULATED. 


ST2=(CODE/ (1 o—U )) FALPHA* (EL ( J-8 DFE L( I) DECC TTC I—8)—-TT(0I))7260) 


ST1=ST2+STt 
CONTINUE 
RETURN 
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SUBROUTINE STRI5 


SUBROUTINE STRI5S 
PRGERKAEEKKE EERE 


SUBROUTINE STR2S5 CALCULATES THERMAL STRESSES OURING THE 


TIME PERIOD (T-24HRS.) TO (T—2HRS.)eTHESE STRESSES ARE 
CALCULATED USING A 15 MINUTE TIME INTERVAL e 


COMMON TITLE1L( 20) oT(15) oAT( 1S) eE(15) oX1(15) sAL(15) sA2(15)eB1(15) 9 
1B2C0 15) » TEMP(1300)+sTT( 3000) .AT15(3000) 0X 115(3000) sEL(3000).ST{( 3000) 
.2sTA(S0O) sATS(50) oX1S5(50) »NUMDAY» NUMEXI eNUMATT sNUMINCs ENCTIMs ALPHA 
BNATTI»NEXIL»sCODEs UeSTI ol MoMT aL M2eL MIS oLMXI 


OOUBLE PRECISION OXIeXI15eXIS 


THE CHANGE IN REDUCED TIME (OXI) AS A FUNCTION ‘OF TIME 
AND TEMPERATURE IS EVALUATED. 


DO 11 J=LM2,LM15 

DXT=X1T15(LMXI)—XI15 (9) 

IFCOXI wLEe X1(1)) GO TO 5S 

DO 3 L=1sNEXI1 

LF(DXI eGEe XI(L) eANO. OXE eLEe XI(L4+1)) GO TO 7 
CONT ENUE 

&1(J5)=E(NUMEXI) 

Go TO 9 

E1(J)=E(1) 

GO TO 9 
E1(J)=DEXP((B1(L) +#A1(L) #OLOG10( OXI) 270.434) 
IF( J eEQ@e LM2) GO TO 11 


THE THERMAL STRESS INCREMENT (STR2) » TOGETHER WITH THE 
TOTAL STRESS (ST1)sIS CALCULATED. 


ST2=(CODE/ (1 --U) ) ¥ALPHA¥ (ELC J—1LIFELC III EC OTTO I-1)-TT0 I) 97260) 
ST1=ST2+STl 

CONTINUE 

RETURN 

END 
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SUBROUTINE XI5M 


SUBROUTINE XISM 
REE ERK AE KEK KK 


THIS SUBROUTINE CALCULATES THE REDUCED TIME VECTOR(XI5SM) 
ASSOCLATED WITH TEMPERATURES 2 HOURS IMMEDIATELY PRIOR TO 
THE TIME AT WHICH STRESSES ARE TG BE DETERMINED. 


COMMON TITLEL(20)oT( 15) sATC 15) sE(15) »XT(15) eAL (15) sA2Z(15) sBIC15) 6 
182015) sTEMP(1300).TT{ 3000) sATI5( 3000) »X115( 3000) .£1(3000).ST( 3000) 
22TA(SO) eATS(50) XI 5(50) »NUMCAY » NUMEX I sNUMATT eNUMINCs INCTIMs ALPHA®s 
SNATTLIeNEXI1e¢CODEsUoSTI oLMoMT ol M2,LMIS oL MXI 

OOUBLE PRECISION DXfsXI15eXI5 

TAC 1)=TT(LM15S) 


TEMPERATURES (TA) AT INTERVALS OF S MINe ARE CALCULATED. 


DO 2 L=2,25 
TACL)=TT(LMIS) + (CTT (LM1548)-TT(LM15) 9/240) *(L-1) 


2 CONTINUE 


SHIFT FACTORS (ATS) AT INTERVALS GF S MINe ARE CALCULATED. 


ATS5(1)=ATIS(LMIS) 

OO 10 L=2:25 

O00 6 J=lsNATTI 

IE (TACL) eh Te TKI) eANDo TA{L) -GEe T(J+1)93 GO TO 8 
6 CONTINUE 

ATS(L)=AT(1) 

GO TO 10 
B ATS(L) =EXP( (A205) *TACL) 82 (5) 9702434) 


10 CONTINUE 


THE REDUCED TIME VECTOR (XIS) IS CALCULATED. 


XES(LI=XTIS(LMIS) 

DO 12 L=2%25 

X15 (L)= (300 0/200) FC Le OSATS(L +1 oOZATS(L~1) D#XIS(L—1) 
12 CONTINUE 

RETURN 

END 
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SUBROUTINE STRS 


SUBROUTINE STRS 
BHR RM HK KKK KEK & 


SUBROUTINE STRS CALCULATES THERMAL STRESSES FOR THE TIME 
PERIOO(T-2HRS-) TO THE TIME AT WHICH STRESSES ARE TO BE 
DETERMINED. 


COMMON TITLEI(20) oT(1S) sATCIS) sEC15) oX1T{15) pAL (15) eAZ(1S) sBI(15).6 
182015) 6 TEMP(1300).TT( 3000) 0AT15( 3000) © X115( 3000) 0£1( 3000) »eST{ 3000) 
2eTA(S50) sATS(50) oXI15S(50) »NUMDAY » NUMEX I »NUMATT »NUMINCs INCTIMsALPHAd 
SNATT Lt eNEXI 1 eCODE »sUsST1 ol MoMT sLM2oLMISoLMXI 

DOUBLE PRECISION DOX1IeX115.XI5 

00 10 J=1525 


THE CHANGE IN REOUCED TIME (DXI) AS A FUNCTION OF TIME 
- ANDO TEMPERATURE [IS EVALUATED. 


DXI=XT1IS(LMX1I—-XIS( J) 

IF(DOXI eLEe XI(1)) GO TA 4 

DO 2 L=1eoNEXI1 

TFCDXI eGEe XI{L) eANDe DOXE eLEe XI(L+1)) GO TO 6 
2 CONTINUE 

E1(J)=E (NUMEXT) 

GO TO 8 

4 ELC JI=E(1) 

GO TO 8 
6 E10 J) =DEXP((81(L)+AL(L) F#OLOGIO (OXI ))70 2434) 


8 IF( J eEQe 1) GO TO 10 


THE THERMAL STRESS INCREMENT (STR2) »eTOGETHER WITH THE 
TOTAL STRESS (ST1)01S CALCULATED 


ST2=(CODE/ (1 o-U )) HALPHA® (E1 (J-1 FEL (I) # CC TAC I-“1I-TACI) 7200) 


ST1i=ST2+STi 
10 CONT INUE 

ST(MT)=ST1 

ST1=0-0 

RE TURN 
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APPENDIX D 


COMPUTER PROGRAM FOR CALCULATING 


THE MODULUS FUNCTION R(é) 
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APPENDIX D 


COMPUTER PROGRAM FOR CALCULATING 


THE MODULUS FUNCTION R(é) 


D.1. Introduction 

Humphreys and Martin (1963)* proposed a method of determining 
thermal stresses in an infinite rigidly constrained slab with tempera- 
ture viscoelastic properties. Assuming the coefficient of thermal 
expansion, Og, Of the slab to be a constant value the stress equation 


suggested by the authors may be expressed as 
t 


6, * 4 sa [ R[E-&"]dT(t') =) 


O 


where R is a time and temperature dependent modulus of the slab. The 
authors presented a numerical method for evaluating the modulus 
function R(&) from known stiffness-reduced time relationships. In this 
appendix, equations used for this numerical evaluation are presented 


and a listing of a computer program developed for the solution of these 


equations is given. 


D.2. Numerical Method for Evaluating R(€) 


Assuming the material is elastic in dilatation, Humphreys and 


Martin presented the following expressions for the numerical evaluation 


of RCE). 
a en ee 


i "Determination of Transient 

“ revs. J.S. and C.J. Martin, 1963. : : 

eae in a Slab with Temperature Dependent Viscoelastic 
Properties”, Transactions of the Society of Rheology, Vol. VII. 
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£(Enei) +5 ¥ (En) [£(0)-f(Ensy-En)] 
See ee eee 


y(En+)) = i 
ae 5 [£(0)-f (Ens )-En)J 
ies! 
i S yy Ly (Eiai) + ¥ (ER) VE CEns 1-244) - EU SnGiee tA 
NEE ane ETE (D=2) 


i= = [f(0) * £(En+7-En) ] 


which is made dimensionless by the substitutions 


ple) = =e (tay) (tae) 

P(e o> P2EC/SE a 

Ho = -B(1-v9)/4(1-v9) 

where 

R(&) = a modulus function, 
Vo = Poisson's ratio, 
ES = initial elastic Young's modulus, 
E(&) = time dependent Young's modulus. 


Using these relationships R(&) can be evaluated provided that the 


latter three listed parameters are defined. 


Das. Program Operations 


The listing of the computer program developed for calculating 
R(E) contains a description of the necessary input parameters. In the 
program a linear relationship is assumed between logarithm of the 
stiffness (E) and the logarithm of reduced time (XI). Therefore, in 


order to obtain accurate moduli, E and XI values should be specified 


using small increments of reduced time. 


; { ¥ Se . 7 {¥ 


isu 1) ~ {0} 2} GF) « $ P Coad)? 


eo ee oe men 


iW yo a 4 Hi 
he Saal jana) (o)2} : ad ne on 


3 2 
* ; : > { “1 : 7 
ar 


[pS quo Geek Pp 31 EC; D+ Card ce 


~ —_— — tex) 
(0) aa 


a . 4 at 
7 , 
anoitutitedse edt yd eeolmotenomib aban ak = 


py rly 


(oA\Q)A DS. = ax 
7 ote\(3)aS = (3)? ye 


~ Fy s 


. Geet = a 


we J pe 

foitsm? esivboa 8 = .)a 4 a 
“ - nee o>} - 

.oltsr e'noeziot =" ov S 

a . 


eutabom e*ynmvoY oitests Isktin ps * ee ' 


» 
/ 


.eiiiubom e'gauoY Inobmeqeb emit - ‘Ma 
' = ; | 

ort ted bobiverq betsulave od 82 (3)4 eqiilenoksaior 220 + x 

$5 unenkh oe 

.boniteb ers sinimiand boteil 29 mip a 


Pita ae 
A, Pé 


Tage: hee » mae ae 

a ‘ i ie > ms Ot 

a - - ei ae era gee meri 

| gaits luotso 3 ite i glove nea. xosucneo aur to niger « s ar 
i 


a a 7 i, ; ont pale =A et: 


‘i ee qi sasezooen ot to oan £ ants. s3n02 2 Gy 


alt ot 


HL treewsod bomuezs zi ail oO] 


Le : in sy, — 
me seat , i) ie it eet: 20 aiid 3% ek et ba 


ole A> 2A ed 
Wiad baittogga adi f ite 298 fi 


ene ao . 


D3 


Operations within the program are outlined in the listing. 
The output consists of stiffness and reduced time values, together with 
corresponding modulus function values. From this output modulus 


function-reduced time plots can be readily established. 
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THIS PROGRAM SOLVES FOR THE MODULUS FUNCTION (R) USING 
A NUMERCTAL PROCEDURE PROPOSED BY HUMPHREYS AND MARTIN; 1963.6 


DIMENSION TITLE LT (20) 6X1140)sE( 40) oF (40) »R( 40) 0 ¥(40) sA1 (40) sB1 (40) 5 


1€(40) »SUMY (40) d 


DOUBLE PRECISION FeRsV¥oCsSUMYsDOXIsEl 


THE FOLLOWING VALUES ARE SPECIFIED. 


NUMMAT — THE NUMBER OF SETS OF DATA SPECIFIED. 


NUMEXI —- THE NUMBER OF STIFFNESS AND REOUCEO TIME VALUES 
IN EACH SET OF DATAe 
TITLE — A TITLE DESCRIBING THE MATERIAL es 


U-—- INITIAL PCISSON*®S RATIOe 
EO — INITIAL STIFFNESS MODULUS. 
XI -—- REOUCED TIME VALUES FROM STIFFNESS—-REDUCED TIME 
. PLOT. 
, E — CORRESPONDING STIFFNESS VALUES FROM STIFFNESS— 
REOUCED TIME PLOT. 


READ( Ses 2)NUMMAT »NUMEXI 
FORMAT (215) 

DO 36 I1f=1 »>NUMMAT 
READ(5s4)(TITLEL( I )ol=15,20) 
FORMAT (20A4) 
READ(5+s6)VeEO 
FORMAT (1FS 02s 1E6 02) 
00 10 L=2.NUMEXE 
READ(5s8)XI(L)eE(L) 
FORMAT (1E601,1E602) 
CONTINUE 


INITIAL VALUES OF Y AND F IN EQUATION(D-2) ARE CALCULATED. 


E(1)=E0 

X1I¢1)=0-0 

U=( 30 O¥ (1 00-2  O¥V) (400 (1 00-V)) 
F(1)=2-0/3-0 

¥CL)=2e 081 60-VIF(RILIZECL)D 

NEX 1 1=NUMEXI-1 

NEXT 2=NUMEX [-2 


00 12 J=2eNEXI1 
A1( J) =—1LeO*ALOGIOCE(JIZECJ+I)) 
B1 (J) =ALOGIOCE( J) I—-AL CJD ¥ALOGLO(XEOS)) 


CONTINUE 
THE MODULLS FUNCTION R IS CALCULATEDe 
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14 
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20 
22 


24 


26 


28 


30 


DS 


DO 28 KK=1.NEXI2 


THE CHANGE IN REDUCED TIME IS CALCULATEO ANO THE PARAMETER 
F IS DETERMINED. 


DO 18 KL=1eKK 

OXT=XICKK +1 )—-XI (KL) 

OO 14 L=2sNEXI1L 

TF(OXI «GE. XE(L) sANOs OXI oLE~ XI(L#1)) GO TO 16 
CONT INUVE 

E1=DEXP((B1(L) +A1(L)*OLOGIO(OXI)) 402434) 

FCKL41)=( 2207320) *({EL/E(1)) 

CONTINUE 


THE DENOMINATOR OF EQUATION(D—-2) IS CALCULATED. 


C(1)=(U/220) €( FC 1 FC KK4+1)) 
IFCKK eEQe 1) GO TO 22 
IL=KK+1 

DO 20 KM=2sKK 
C(KM)=(U/20e0) (FC IL )—-FCIL~-1)) 
[L=Itt-1 

CONTINUE 


TF (KK eGTe 1) GO TO 24 
VCKK4+1)=€02 0/3 200) CEC KKELD/SECKK) DD +0CC1 VY 01999701 2-0—-C(1)) 


RCKK +1) =CY (KK 41) *#E(1)9)70220-(120-V)) 
SUMYCLI=YCKKt1) YC 1) 

GO TO 28 

LA=KK~-1 

LB=KK 

SUMYF=0-0 


THE NUMERATOR OF THE SECOND TERM OF EQUATION(O-2) [IS CALCULATED. 


DO 26 MW=1eLA 
PRODYF=C(LB)*SUMY ( MW) 
SUMYF =SUMYF + PRODYF 
LB=LB-1 


CONT (NUE 
Y (KKt1)=00 200/320) (EC KK41) SEC 1) FSUMYF4+ (COLI *VOKK) DD 70120-CO1)D) 


RUKK +1) = CY (KK 41) EC 1) 70200401 00-V)) 

SUMY (KK) =YCKK41 )+YCKK) 

CONTINUE 
TITLELsREOUCED TIMEsSTIFFNESS ANO THE MODULUS FUNCTION 
(R) ARE WRITTENe 


WRITE(6e30)(TITLEI (I) of=1 520) 
FORMAT (6X0 20A4) 
OO 34 I[=1+eNEXI? 
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32 
34 


36 


WRITE (6232) XI(T) eEC I) RI) 
FORMAT (5X9 1E7 02 25X 0 LEVe 405Xe LED 2 4) 
CONT INUE 


OPERATIONS RETURN TO BEGINNING OF PROGRAM WHERE A SECOND 
SET OF DATA MAY BE SPECIFIEDe 


CONT INUE 
STOP 
END 
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APPENDIX E 


LISTING OF A COMPUTER PROGRAM FOR 


THERMAL STRESS PREDICTIONS USING VISCOELASTIC 


ANALYSES WITH CONSTANT TIME INCREMENTS 
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APPENDIX E 


LISTING OF A COMPUTER PROGRAM FOR 
THERMAL STRESS PREDICTIONS USING VISCOELASTIC 


ANALYSES WITH CONSTANT TIME INCREMENTS 


ee ae Introductory Remarks 


This appendix contains the listing of a computer program used to 
perform thermal stress calculations in which viscoelastic analyses are 
employed. Stress equations for the analyses are those given in 
APPENDIX C. Thermal stresses are computed using a numerical inte- 
gration scheme which incorporates a constant time increment. In the 
listing this time increment is defined by the term DELT. With the 
exception of differences in numerical integration schemes, operations 


within the program parallel those described in APPENDIX C. 
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THIS PROGRAM USES VISCOELASTIC ANALYSES TO CALCULATE 
THERMAL STRESSES IN AN ASPHALTIC CONCRETE SUBJECTED TO A 
KNOWN TEMPERATURE FIELDe 


STRESSES ARE CALCULATED AT EQUAL TIME INTERVALS 
OVER THE TOTAL TIME PERIOD IN WHICH THE TEMPERATURE FIELD 
IS OEFINED. 


THE PROGRAM CAN BE USED TO DETERMINE THE THERMAL STRESSES 
ASSUMING ONE OF THE FOLLOWING CONDITIONS$ 
1) RESTRAINED BEAMsE=FUNCTION(REOUCED TIME) 
.2) RESTRAINED SLABsE=FUNCTION(REOUCED TIME ) 
3) RESTRAINED SLABeR=FUNCTION( REDUCED TIME) 


. MAIN PROGRAM 
; EREKKEKKAKKKEKER 


COMMON TITLEL (20) ¢T(15) eAT(IS) eECISIVsXLCI1S)2A1(15) eA2(15)»sBiC IS)» 
182(15) «6 TEMP( 1300) »TT( 3000) eA( 3000) »X11( 3000) »€1( 3000) »-ST( 3000) 6 
2NUMDAY » NUMEXI sNUMATT eNUMINC op ENCTIMsDELT » ALPHAsNNeoNATTIsNEX IL oNTo 
3CODE eV 

OOUBLE PRECISION OX!sXIt 

CALL READIN 

CALL LINEAR 

CALL XIS 

CALL STRESS 

WRITE (Gol IC TITLEICI ) of=1520) 

FORMAT (SX 20A4/ ) 


TEMPERATURES AND STRESSES ARE WRITTEN AT TIME INCREMENTS 
CORRESPONOING TO THOSE AT WHICH TEMPERATURES ARE SPECIFIED. 


00 3 LT=1leNNoNT 

WRITE (6Ge2) TT(LT) eo STI(LT) 
FORMAT(12X%e2F100e1) 

CONT INUE 

STOP 
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ES 


SUBROUTINE READIN 


SUBROUTINE REAOIN 
PRREKEEEKKKEK ERE RE 


THIS SUBROUTINE READS IN ALL DATA NECESSARY FOR THE 
GENERATION OF THERMAL STRESSESe THIS DATA CAN BE DIVIOEO INTO 
THE FOLLOWING THREE CATEGORIES$ 
1) THE ASSUMPTION OF A RESTRAINED BEAM OR SLAB» 
CONDITIONS 1¢62 OR 3 MENTIONED PREVIOUSLY. 

2) MATERIAL FROPERTIES. 

3) TEMPERATURE AT A GIVEN DEPTH IN THE ASPHALTIC 
CONCRETE AS A FUNCTION OF TIME. 


COMMON TITLE1(20)6T( 15) sAT( 15) .EC 15) oXT (15) sAL (15) sA2(15).B81015) 9 
182015) »5TEMP(1306)»5TT( 3000) »A( 3000) eX11( 3000) £14 3000) »ST( 3000) >, 
2NUMDAY sNUMEXI s NUMATT »NUMINCs INCTIMsDELT sALPHAeNNe NATIT Le NEXLLoNTs 
3CODE»U 

DOUBL'E PRECISION DXIeXTI1 


THE INPUT PARAMETERS ARE DEFINED AS FOLLOWS: 
CODE - CODE=1.0 OENOTES CONDITION le 
CODE=2-0 DENOTES CONDITIGN 26 
COOE=3-0 DENOTES CONDITION 3-6 
U - POISSONS RATIO. ONLY IF CONDITION 2 {IS SPECIFIED 
A VALUE OF POISSONS RATIO MUST BE READ INe 
NUMDAY — THE NUMBER OF DAYS IN WHICH THE TEMPERATURE 
REGIME IS SPECIFIED. 
NUMEXI — THE NUMBER OF VALUES SPECIFIED FROM THE STIFFNESS MOOULUS-— 
RECUCED TIME RELATIONSHIP. 
NUMATT — THE NUMBER OF VALUES SPECIFIED FROM THE SHIFT FACTOR- 
TEMPERATURE RELATIONSHIP. 
INCTIM — THE TIME INCREMENT IN WHICH THE TEMPERATURE 
REGIME IS SPECIFIED.(HRS) 
DELT -— THE TIME INCREMENT IN WHICH THE THERMAL STRESSES 
ARE TO BE CALCULATEDe(MINSe) 
THE COEFFICIENT OF THERMAL EXPANSION AND CONTRACTION 


ALPHA —- 
OF THE ASPHALTIC CONCRETEe({IN/IN/DEG F) 
E - THE STIFFNESS MODULUS OF THE ASPHALTIC CONCRETE.(PSI) 
THESE VALUES ARE SPECIFIED IN DESCENDING ORDER 
OF MAGNITUDE. 
XL —- REDUCED TIME«(SECS) THESE VALUES ARE SPECIFIED IN 


ASCENDING ORDER OF MAGNITUDE AND MUST BE COMPATIBLE 


WITH THE E VALUES. 


AT - SHIFT FACTORS» 
T — TEMPERATURES COMPATIBLE WITH THE SHIFT FACTORS.«(0EG F) 


TEMP — THE TEMPERATURES AT A GIVEN DEPTH WITHIN THE 
ASPHALTIC CONCRETEe(DEG F) 


READ(Se 1 (TITLEL( 1)s1=1 020) 
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FORMAT (20A4) 

READ(5»22) CODE 

FORMAT(1F56e0) 

IF (CODE-22 3982498 

READ(5»6)U 

FORMAT(1F5e2) 

GO TO 10 

U=0.0 

READ (5012) ( NUMDAY sNUMEXI »>NUMATT »INCTIM »DELT »ALPHA) 
FORMAT (411061 FI0e2,1E8e1) 
READ(5014)¢€E( 1) ol=1+.NUMEXI ) 
FORMAT(C1L1LE7 e2) 
READ(S5e616) (XI (1) T=i1.+NUMEXI ) 
FORMATCIIE71) 
READ(5218) (ATC I) eo T=1.eNUMATT) 
FORMAT(9E7e2) 
READ(5s20)(T( 1) sL=1 »NUMATT) 
FORMAT (9F8.0) 
NUMINC=24/TUNCTIM 
NUMTEM=NUMINC#NUMDAY #1 

OO 24 L=1+NUMTEM 
READ(S5e22)TEMP(L) 

FORMAT (1F10-0) 

CONT INUE 

RETURN 

END 
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SUBROUTINE LINEAR 


SUBROUTINE LINEAR 
EREKEEEKEKKEAE EEK EEE 


A LINEAR RELATIONSHIP IS ASSUMED BETWEEN LOG(E) AND 
LOG(XI)»sAND LOG(AT) AND T VALUESe THIS SUBROUTINE CALCULATES 
THE A AND 8B COEFFICIENTS OF THE RESULTING LINEAR EQUATIONS. 


COMMON TITLE1 (20) o T(15) eAT(LS) oE(1LS) oX1(15)sA1(15)sA2(15) eBiCi5)» 
182(015) 5 TEMP(1300)5TT( 3000) eA( 3000) eX11( 3000) s£1( 3000) »ST( 3000)» 
2NUMDAY » NUMEX Is» NUMATT eNUMINCs INCTIMsDELT s ALPHAsNNoNATT1L»sNEXIL ONT»? 
3CODE eU 

DOUBLE PRECISION OXI.XI1 

NEXI1=NUMEXI-1 

NATTI=NUMATT—-1 


THE A AND B COEFFICIENTS FOR THE LOG(E)—-LOG(XI) 
‘EQUATIONS ARE EVALUATED. 


00 1 J=1eNEXI1 
A1L(J)=—-1eO*ALOGIO(E(CIIJSE( Jt1)) 
B1(J)=ALOGIO(E( J) ALCS) FALOGIO(XE(I)) 
CONTINUE 


THE A ANO 68 COEFFICIENTS FOR THE LOG(AT)—-—T EQUATIONS 
ARE EVALUTEDe 


DO 3 J=1sNATTI 
A20J)=CALOGIO (ATC J) )—ALOGIOCAT( S41) I70TCII—-TC S41 DD 
B2€I)=ALOG1IO(CAT( J) I—A2C ID #TCU) 

CONT INUE 

RETURN 

END 
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SUBROUTINE XIS 


SUBROUTINE XIS 
SEEKER EEE KKE EK EK 


THIS SUBROUTINE CALCULATES THE REDUCED TIME VECTOR(XI1) 
ASSOCIATED WITH THE INPUT TEMPERATURE FIELD.» THIS VECTOR IS 
DETERMENEO AT TIME INTERVALS GF OELT OVER THE CCOMFLETE 
TIME PERIOO IN WHICH THE TEMPERATURE FIELD IS DEFINED. 


COMMON TITLE (20) oT (15) sATCIS) FEC 15) oXT(15)eA1 (15) sA2(15) oBIC15)>s 
1B2(15) 5 TEMP( 1300) 6 TT( 3000) sA( 3000) »XI1{ 3000) .£1(3000) »ST( 3000)» 
2NUMDAY » NUMEX IT »NUMATT eNUMINCo INCTIMseDELT eo ALPHAsNNeNATT I oNEX IL NT» 
3CODE »U 

DOUBLE PRECISION OXI,XIl 

TIME=INCTIM 

NT=(TIME*606)/DELT 

NTL=NT+4+1 

N=NUMD AY * NUMI NC 


TEMPERATURES (TT) AT INTERVALS OF DELT ARE CALCULATED. 


TT(1)=TEMP( 1) 

L=2 

OO 3 M=1.N 

00 1 LM=2-NTI 

TT(L)=TEMP (M) #( 0 TEMP( M41 )—TEMP(M) D/NT) &0LM—1) 
L=L+1 

CONTINUE 

CONTINUE 


SHIFT FACTORS (A) AT INTERVALS OF DELT ARE CALCULATED. 


NN=N¥NT#1 

O00 9 L=1.NN 

DO 5 J=1eNATTI 

IFCTT(L) eb Te TOI) eANDo TT(L) eGEe TK J+1)) GO To 7 
CONTINUE 

A(LI=ATC(1) 

GO TO 9 

ACL)J=EXP( (A209) *#TT(L) +8205) 9702434) 


CONTINUE 
THE REDUCEO TIME VECTOR (XI1) IS CALCULATED. 


X12101)=0.0 

DO 11 L=2eNN 

X11 (L) =(DELT#60 0/20) #CL SACL) HL eZA(L—-1) DX ELCL-1) 
CONTINUE 

RETURN 


END 


AY Ln) 
Pa ry, Uae a 
7 7 ’ / 
a ie 7 
1 t 7 A 
ay a 
f R 
‘= ‘ ‘1. 
4 
q > 
J I ro 
7 i 
t) 
A 
: — 
1 i] : i : 
z 
4 / tl 
i ; 
i — 
aT 
‘ " i 
~ “ 


4 4 


r ye) 4 . Paws Li 
2X Sul voMeVve 


4 ~ 
~— WAS TOOMAVE 
“$9 08 FO SH Se HERE OE 
tepayaety AT (PPA AU DAAD On PPUOROUE RANT 1) 
eh Woroay \< “es jad ULAR IBAIT Tyee att HTje DATATIOSeA } 

a,s9Ma9 | DHVO. 4438 Ao SAVE Be27 TA OSHIMASTSO, 9° | 

HPAsn 2) ae LW?RADAMSYT SHT BOTHW HF OC RAG BMLT mh 

- rar 

eC2RIFMCLSEL PSA Loy eS 2s Kw FSR PO eS Aet@4¥? 409) 232087 WoO? 7 

. (200 } +¢ , t.3 ’ i ‘£95 i Ke - (000k 1AsAOUDE ITT LODEL AAMAT FRU. 
ete i}: t etitAs 1 tho Tbe BM 1 IAT OHTA F TA™ la Mercian the YAR - 

a * : vs 3009 

TRI KO MOFZIIIAG, 3jBUOG 
_ “wE¥oniesarr™~ 
” “Aaanaw cowamtta arn, 
3 ia seve kT Ww 
f au tian eae 
1 Li u i 7 + ™ 
SOStAMUI IAS SP i320 4) SsRVASIMS VA OTT F SRP EAMEMMAT 9. 5) iis 
} ; ,, ~ 
: ie Caenarely yet 
/ ca ; a ar & go y 
§ 7: £FHeSe44 8 oe 


C8 CTRL I OSE 19 28MB TD De OO Sa 
yr m ra 

Sie 4 . ) Seeonraa; 

| ofl then es 

| F is a ahee we a. 
SABTARUDIAD BRA T4130 3 BIAVMATKT FA CA) 200785 TaIR2 en ' 


ay ie ; ps , em 
a) [A Ahaha @ OG 
ay i | T.| _ETTAMs Pet. 8 98 
' ¥ OT GO UCL 9 O6 ae OAs ter ae 
. Ae? | rh | 7 i a ong r aa eon 


vd 7 ai . (ER ONL CL ISEH CITT t ; eee . 
vo - Ara) j y mae 5  BEMITME 
os 7 \ Ca P 9 
wre hit - @OBTA UD gAD 21 (21K? ROTOSV away a2uo38 a 
y 


_ 
ae | a7 a 


Hie . ; y 
ace itnne ete svANede tan 


(oviaiay (avtaltay (a) 


(ol aia) fa 


Anan 


10 
11 


SUBROUTINE STRESS 


SUBROUTINE STRESS 
SOR IOI $8 ORR Fok toto 


SUBROUTINE STRESS CALCULATES THERMAL STRESSES AT TIME 
INTERVALS GF DOELT. 


COMMON TITLE1 (20) oT( 15) sATCI15)¢E0C15)>XI( 15) sA1(15) sA2(15) »BIL15) 6 
182015) » TEMP( L300) oTT( 3000) 6A(3000) oXI143000) £1 (3000) »ST(3000)> 
2NUMDAY ¢ NUMEX T>NUMATT »eNUMINCo ENC TIM sDELT sALPHA sNNeoNATT I »sNEXIL ONT 
3CODEsU 

DOUBLE PRECISION OXIeXIt 

STR1=020 

00 11 KK=1.NN 

KT=KK-1 


THE CHANGE IN REDUCED TIME (OXI) AS A FUNCTION OF TIME 
"AND TEMPERATURE IS EVALUATEDe 


DO 5S [=15KK 

OXIT=XIL(KK)-XTICT) 

IF(OXI elEe XI(1)) GO TO 3 

00 1 L=1s,NEXI1 

IF( OXI eGEe XI(L) eANDe OXI eLEe XI(L#1)) GO TO 4 
CONT [NUE 

E1(1)=ECNUMEXT) 

GO TO 5 

E1( I )=E(1) 

GO TO §$ 

E101 =OEXP((B1(L)FAL(L) #OLOGIO(0XI))/00434) 
CONTINUE 

IF(KK eEQe. 1) GO TO 9 


THE THERMAL STRESS INCREMENT (STR2)sTOGETHER WITH THE 


TOTAL STRESS (ST)0IS COMPUTED. 


DO 7 KL=1.KT 

STR2= (CODE/( beo—-U)) #ALPHAX(ELCKLIFEL(KL ALD) * (CTT CKLY-TTCKL +1972 60) 
STRIL=STR2+tSTRI 

ST(KK)=STRI 

GO TO 10 

ST(KK}=020 

STR1=02-0 

CONT TINUE 

RETURN 

END 
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